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ПАМЯТИ УЧЕНОГО 





27 января 2021 года на 84-м году жизни скоропостижно скончался замечательный 
ученый, инженер и организатор, доктор технических наук, профессор Анатолий 
Андреевич Рыжкин. 

Анатолий Андреевич Рыжкин — заслуженный деятель науки и техники РФ, 
доктор технических наук, профессор. 

Анатолий Андреевич родился 21| января 1938г. Окончил Ростовский-на-Дону 
институт сельскохозяйственного машиностроения в 1960 г. Вся его более чем 60- 
летняя трудовая деятельность связана с родным вузом. 

В переходное для РИСХМА и всего российского образования время возглавил вуз 
— в феврале 1988 г. был избран ректором и занимал эту должность до 2007 г. 
Анатолий Андреевич Рыжкин внес большой вклад в преобразование отраслевого 





института сельхозмашиностроения в технический университет: была создана 
современная база науки и учебного процесса, увеличено число специальностей, расширена аспирантура и 
открыта докторантура, развивались объекты социальной инфраструктуры. В результате РИСХМ превратился в 
региональный центр образования, науки и культуры. В 1992 году по его инициативе и под личным 
руководством решением Правительства Российской Федерации РИСХМ был преобразован в Донской 
государственный технический университет (ДГТУ). 

Анатолий Андреевич Рыжкин — известный ученый в области трения и износа в условиях резания 
материалов. Он основал и возглавил научное направление «Физико-химические и термодинамические основы 
управления износостойкостью — инструментальных режущих — материалов». Прикладные —научно- 
исследовательские работы А. А. Рыжкина были направлены на создание технологий упрочнения 
инструментальных материалов. Результаты НИОКР внедрены на различных предприятиях Ростовской области 
и широко используются при обучении студентов по направлению «Конструкторско-технологическое 
обеспечение машиностроительных производств». 

А. А. Рыжкин — автор более 450 научных и методических работ. Более 50 из них изданы за рубежом 
на английском, немецком, польском, испанском и других языках. На протяжении многих лет профессор 
А. А. Рыжкин руководил аспирантурой и докторантурой, являлся председателем совета по защите докторских 
диссертаций при Донском государственном техническом университете. Анатолий Андреевич лично подготовил 
16 кандидатов и 5 докторов технических наук. 

Более 15 лет Анатолий Андреевич Рыжкин возглавлял Совет ректоров вузов Ростовской области, 
ставший под его руководством действенным общественно-государственным органом управления сферой 
образования Ростовской области. 

В последние годы Анатолий Андреевич заведовал кафедрой «Металлорежущие станки и инструменты» 
и являлся вице-президентом Ассоциации инженерного образования России, членом проблемного Совета по 
машиностроению при Министерстве образования и науки РФ. 

За многолетний и добросовестный труд, большой вклад в решение проблем высшего образования, 
масштабную научную и общественную деятельность Анатолий Андреевич Рыжкин был награжден различными 
государственными наградами и удостоен почетных званий. Ветеран труда, почетный профессор ДГТУ, 
заслуженный деятель науки и техники РФ, почетный работник высшей школы РФ. Он лауреат премий 
Президента РФ и Правительства РФ в области образования; обладатель ордена «За заслуги перед Отечеством 
ГУ степени», ордена Дружбы, «Знака почета», медали «За доблестный труд», ордена «За заслуги перед 
Ростовской областью». 
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Введение. В настоящее время остро стоит проблема получения облегченных бетонных и железобетонных 
изделий и конструкций с улучшенными структурой и характеристиками. Этого можно добиться 
центрифугированием или более совершенным способом — виброцентрифугированием. В то же время влияние 
центробежных и центростремительных сил инерции в данных видах технологий приводит к различиям свойств 
по сечению бетонных изделий и конструкций. Чтобы отразить это в расчетах, необходимо экспериментально и 
аналитически исследовать качественную и количественную картины таких различий в характеристиках 
бетонов, полученных по разным технологиям. 

Материалы и методы. В исследовании применялись усредненные по сечению характеристики бетона — 
«интегральные характеристики бетона». Применяемые сырьевые материалы: портландцемент марки 500, 
щебень фракции 5-20, песок средний. Было изготовлено и испытано девять контрольных образцов кольцевого 
сечения, полученных вибрированием, центрифугированием и виброцентрифугированием. Суть методики 
заключалась в том, что каждый изготовленный экспериментальный контрольный образец использовался сразу в 
нескольких видах испытаний. Из общего кольцевого сечения каждого образца выделялось по три условных 
квадранта из которых выпиливались стандартные образцы малого размера. Впоследствии их испытывали на 
осевое сжатие, растяжение, а также растяжение при изгибе. Применялось следующее испытательное 
оборудование: механический пресс с электронным управлением ИПС-10 — для испытания призм на сжатие и 
разрывная машина Р-10 — для испытания образцов на осевое растяжение. Для измерения деформаций бетона 
применялись тензодатчики и индикаторы часового типа. Для получения деформативных и прочностных 
свойств бетона, в том числе полных диаграмм деформирования с нисходящими ветвями, применялись также 
осциллографы. 

Результаты исследования. Проанализированы результаты расчета интегральных конструктивных 
характеристик бетонов, полученных вибрированием, центрифугированием и виброцентрифугированием. 
Установлено, что вследствие влияния центробежных и центростремительных сил инерции при 
центрифугировании и виброцентрифугировании характеристики бетона по сечению становятся различными. В 
ряде случаев различия эти могут быть весьма значительными. Разработаны и апробированы: новая методика 
оценки зависимости интегральных (усредненных по сечению) конструктивных характеристик бетона 
(плотность, кубиковая и призменная прочности на осевое сжатие); предельные деформации при осевом сжатии; 
прочность на осевое растяжение и растяжение при изгибе; предельные деформации при осевом растяжении; 
модуль упругости; диаграмма «напряжения б,_деформации =» при сжатии; диаграмма «напряжения бь— 
деформации ши» при растяжении от технологии изготовления (вибрирование, центрифугирование, 
виброцентрифугирование). 

Обсуждение и заключения. По результатам исследований сформулированы выводы о положительном эффекте 
предложенной технологии совместного вибрирования и центрифугирования, заключающемся в улучшении 
интегральных конструктивных характеристик и структуры бетона от вибрирования к центрифугированию и от 
центрифугирования к виброцентрифугированию. 


Ключевые слова: вибрирование, центрифугирование виброцентрифугирование, расчет колонн, вариатропная 
структура, интегральные характеристики бетона, предельные деформации, прочность на сжатие, модуль 
упругости. 
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Введение. Известно, что получить бетоны с улучшенными структурой и характеристиками, при 
сравнении с теми, которые получаются центрифугированием и вибрированием, можно путем 
виброцентрифугирования ”” [1-11]. 

Однако в немногочисленных работах по бетонным и железобетонным конструкциям кольцевого 
сечения исследовалось влияние вида технологии на усредненные (общие) по сечению характеристики бетона 
[12—15]. В то же время очевидно, что вследствие влияния центробежных и центростремительных сил инерции 
при центрифугировании и виброцентрифугировании характеристики бетона по сечению становятся 
различными [16]. 

В настоящей работе экспериментально и аналитически исследуются качественная и количественная 
картины таких различий в характеристиках бетонов, полученных по разным технологиям. Очевидно, что для 
расчетного учета этих различий понадобится ввести некие усредненные характеристики по сечению элементов. 
Введем с этой целью термин «интегральные характеристики бетона». 

Материалы и методы. Всего произведено и испытано девять контрольных образцов кольцевого 
сечения, изготовленных вибрированием, центрифугированием и виброцентрифугированием. Размеры данных 
образцов: внешний диаметр Р — 450 мм, внутренний диаметр & — 150 мм; общая высота Н — 1200 мм. 

Применяемое оборудование и методы испытаний описаны в [8—11]. 

В качестве заполнителя применялся щебень фракции 5—20, что приближает свойства получаемого 
бетона к свойствам мелкозернистого бетона. 

В экспериментах варьировался вид технологии изготовления, что зафиксировано в обозначениях 
образцов: вибрирование — В, центрифугирование — Ц, виброцентрифугирование — ВЦ. 

Рассматривалась задача оценки зависимости интегральных (усредненных по сечению) конструктивных 
характеристик бетона (плотность, прочность на осевое сжатие (кубиковая и призменная); предельные 
деформации при осевом сжатии; прочность на осевое растяжение и растяжение при изгибе; предельные 
деформации при осевом растяжении; модуль упругости; диаграмма «напряжения 6, деформации &5» при 
сжатии; диаграмма «напряжения бы деформации 5» при растяжении) от технологии изготовления 
(вибрирование, центрифугирование, виброцентрифугирование). 

Результаты исследования. Методика испытаний отличалась тем, что каждый изготовленный 
экспериментальный контрольный образец использовался сразу в нескольких видах испытаний. Контрольные 
образцы в единичном экземпляре были отобраны и испытаны на 7, 28 и 180 сутки. 

Из общего кольцевого сечения каждого образца выделялось по 3 условных квадранта 4, В и С, из 
которых выпиливались образцы малого размера. Впоследствии их испытывали на осевое сжатие, растяжение, а 
также растяжение при изгибе (рис. 1, 2). 

Четыре образца-куба с ребром 15см выпиливались из квадранта А для проведения испытаний на 
сжатие и растяжение (уровни № |-4), одна призма (15х15х60 см) — для проведения испытаний на растяжение 
при изгибе (уровень №5). Для проведения испытаний на осевое сжатие две призмы (15х15х60 см) 
выпиливались из квадранта В (уровни № 1—2). Далее три призмы (15%х15х60 см) выпиливали из квадранта С для 
проведения испытаний на осевое растяжение (уровни № 1-2). 

После испытания кубиков на осевое сжатие получены значения А». призм на осевое сжатие — 
значения А», &ю. Вы, &ыю, Еь=Еы и диаграммы деформирования «6,—&»», а после испытаний призм на осевое 
растяжение — значения А», и диаграммы деформирования «6,—&ы» и после испытания призм на растяжение при 


изгибе — значения А». 


' Аксомитас, Г. А. Прочность коротких центрифугированных колонн кольцевого сечения с продольной арматурой класса Ат-У при 
кратковременном сжатии : дис. ... канд. техн. наук / Г. А. Аксомитас. Вильнюс : ВИСИ, 1984. 261 с. 

* Петров, В. П. Технология и свойства центрифугированного бетона с комбинированным заполнителем для стоек опор контактной сети: 
дис. ... канд. техн. наук / В. П. Петров. Ростов-на-Дону : РИСИ, 1983. 175 с. 

°Раджан Сувал. Свойства центрифугированного бетона и совершенствование проектирования центрифугированных железобетонных стоек 
опор ЛЭП: дис. ... канд. техн. наук / Сувал Раджан. Ростов-на-Дону : РГАС, 1997. 267 с. 
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Рис. 1. Схема изготовления образцов бетона малого размера из квадрантов по высоте экспериментальных контрольных 
полноразмерных образцов кольцевого сечения для расчета интегральных характеристик: 
а — квадрант 4; 6 — квадрант В; в — квадрант С 
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Рис. 2. Опытные образцы бетона малого размера из квадрантов А, В, С экспериментальных контрольных 
полноразмерных образцов кольцевого сечения для анализа интегральных конструктивных 
характеристик а — квадрант А; 6 —квадрант В; в —квадрант С 


Использовалась методика испытаний по ГОСТ 10 180. Применялось следующее испытательное 
оборудование: ИПС-10 — для испытания призм на сжатие и Р-10 — для испытания образцов на осевое 
растяжение. 

Для измерения деформаций бетона применялись тензодатчики с базой измерения 50 мм и индикаторы 
часового типа с ценой деления 0,001 мм. 

Для получения деформативных и прочностных свойств бетона, в том числе полных диаграмм 
деформирования с нисходящими ветвями, эксперименты проводили, используя постоянную скорость 
деформирования. 

С этой целью, кроме тензометрических датчиков, применялись также осциллографы. 

Для испытаний с одной и той же скоростью подачи нагрузки был выбран шаг нагружения 0,1А, а 
деформации призм повышались с шагом 0,1» (рис. 3). 
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Рис. 3. Режимы испытаний опытных малоразмерных образцов: 
а — поэтапно, с постоянной скоростью нагружения (Ло=соп$$); 
б — поэтапно, с постоянной скоростью деформирования (Л==соп$1) 


Заключительный режим испытаний состоял в возрастании нагрузки до максимума и последующем ее 
снижении в процессе увеличения деформации. Таким образом во время испытаний фиксировалась нисходящая 
ветвь диаграмм «6 — =» бетона, имеющая довольно четкое очертание примерно до значения с = 0,8 А как при 
сжатии, так и при растяжении. Впоследствии зависимость приобретала весьма неустойчивый характер. 

Проанализировав полученные результаты, можно сделать выводы о влиянии технологии изготовления 
образцов. Результаты экспериментальных исследований изменений интегральных характеристик опытных 
образцов бетона, в зависимости от технологии изготовления, представлены на рис. 4—10. 
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Рис. 4. Влияние технологии изготовления на прирост кубиковой прочности бетона при сжатии 
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Рис. 5. Влияние технологии изготовления на прирост призменной прочности бетона при осевом сжатии 
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Рис. 6. Зависимость прироста прочности на растяжение при изгибе бетона от технологии изготовления 
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Рис. 7. Зависимость прироста прочности на осевое растяжение бетона от технологии изготовления 
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Рис. 8. Зависимость снижения предельных деформаций при осевом сжатии бетона от технологии изготовления 
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Рис. 9. Зависимость снижения предельных деформаций при осевом растяжении бетона от технологии изготовления 
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Рис. 10. Влияние технологии изготовления на прирост модуля упругости бетона 


Обсуждение и заключения. Влияние технологии изготовления образцов на плотность во всех 
возрастах было минимальным (около 2 %), поэтому при расчетах показателем «плотность» можно пренебречь. 


Виброцентрифугированные бетоны по показателям «прочность при сжатии» и «прочность при 
растяжении» продемонстрировали более высокие значения, нежели вибрированные и центрифугированные 
бетоны, а именно: призменная и кубиковая прочности на сжатие — до 22,0 %, прочность на осевое растяжение 
и растяжение при изгибе — до 27,0 %. 

Вследствие упорядочения протекающих процессов гидратации цементного камня увеличение 
прочностей при сжатии и растяжении бетонов с одновременным вибрированием и центрифугированием, по 
сравнению с центрифугированными и вибрированными образцами, с возрастом имеет место в исследуемом 
диапазоне возраста бетона (7 — 180 суток) и находится практически в одинаковых пределах во всех возрастах. 

У виброцентрифугированных образцов происходит незначительное снижение (до 6 %) предельных 
деформаций при осевом сжатии и растяжении. Это соответствует максимальной прочности исследуемых 
бетонов. 

Виброцентрифугированные бетоны, в сравнении с вибрированными и центрифугированными, 


продемонстрировали наименьшие предельные деформации в любом возрасте. 
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Показатели «модуль упругости при сжатии Еь» и «модуль упругости при растяжении Евы» во всех 
возрастах бетона были до 8,0 % больше у виброцентрифугированных бетонов, чем у вибрированных и 
центрифугированных. 

Увеличение прочности, с параллельным снижением предельных деформаций, являлось причиной 
повышения Еь и Еыу бетонов с одновременным вибрированием и центрифугированием, нежели у бетонов с 
одним из видов уплотнения. Данный факт влиял на диаграмму «напряжения-деформации» — максимум 
смещался вверх и влево. 

Отличия диаграмм «напряжения-деформации» при сжатии и растяжении, характерные для бетонов с 
одновременным центрифугированием и вибрированием, во всех возрастах следующие: 

— возрастание прочности и уменьшение предельных деформаций (максимум смещается вверх и влево); 

— увеличение начального модуля упругости (возрастание в начале координат угла подъема). 

Во всех возрастах бетонов была характерна следующая тенденция: увеличение подъемистости 
восходящей ветви диаграмм, снижение нисходящей ветви у центрифугированных и виброцентрифугированных 
бетонов по сравнению с вибрированными. 

По данному исследованию целесообразно сделать следующие заключения. 

1. Исследования интегральных конструктивных характеристик при сжатии и растяжении 
рассмотренных видов бетонов различных технологий изготовления в возрастах 7-28 и 180 суток показали: 

— характеристики бетонов улучшаются от вибрирования к центрифугированию и от 
центрифугирования к виброцентрифугированию; 

— увеличение во всех возрастах прочностей при сжатии и при растяжении (до 23 %); 

— снижение всех предельных деформаций (до 8 %); 

— возрастание модулей упругости при различных видах нагружения (до 10 %). 

2. Для всех интегральных диаграмм деформирования «напряжения-деформации» бетонов с 
одновременным вибрированием и центрифугированием свойственно: 

— перемещение максимума вверх и влево; 

— возрастание в начале координат угла подъема; 

— увеличение подъемистости диаграммы в восходящей ветви. 

3. Посредством проведенных многочисленных экспериментальных исследований установлено, что 
одновременное вибрирование и центрифугирование способствует получению бетонов с улучшенными 
структурой и характеристиками, нежели бетоны, полученные только одним видом воздействия — 
центрифугированием либо вибрированием. 
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Трехмерная интегральная модель сухого трения для движения 
прямоугольного корпуса 





М. С. Салимов, И. В. Меркурьев 
ФГБОУ ВО «Национальный исследовательский университет «МЭИ» (г. Москва, Российская Федерация) 


Введение. Исследуется трехмерная модель сухого трения при взаимодействии прямоугольного тела и горизон- 
тальной шероховатой поверхности. Предполагается, что отрыв корпуса от горизонтальной поверхности отсут- 
ствует. Движение тела происходит в условиях комбинированной динамики, когда помимо продольного движе- 
ния тело участвует в верчении. 

Материалы и методы. Предложены дробно-линейные аппроксимации Паде, которые заменили громоздкие 
аналитические выражения, наиболее точно описывающие движение тел по шероховатым поверхностям. Предло- 
жены новые математические модели, описывающие скольжение и верчение тел с прямоугольным основанием. 
Результаты исследования. Разработаны и научно обоснованы аналитические выражения главного вектора и 
момента сил трения для прямоугольных площадок контакта. Разработана модель трения, которая учитывает 
взаимосвязь между скоростями скольжения и верчения, позволяющая находить решения для зависимостей Па- 
де. После численного решения уравнений движений, получены и построены зависимости скорости скольжения 
и угловой скорости от времени. Построены графики зависимостей сил трения и их момента от угловой скоро- 
сти и скорости проскальзывания, которые позволили сравнить интегральную и нормированную модели трения. 
Результаты сравнения показали хорошее соответствие интегральной модели и модели на основе аппроксимаций 
Паде. 

Обсуждение и заключения. Полученные результаты позволяет учесть динамическую связь компонентов, кото- 
рая определяет силовое взаимодействие прямоугольного корпуса и горизонтальной поверхности. Эти результа- 
ты могут быть использованы в мобильной робототехнической сфере. Анализируемое движение корпуса проис- 
ходит за счет управления движением материальной точки внутри корпуса. Такие мобильные роботы могут ис- 
пользоваться при решении широкого класса задач: при создании автономных роботов для исследования косми- 
ческого пространства и планет; при диагностике и лечении в части прохождения по сложным структурам вен и 
артерий; при исследованиях под водой, в местах больших перепадов температур; при подземных работах. 


Ключевые слова: сухое трение, прямоугольное тело, твердое тело, динамика, скольжение, верчение, сила тре- 
ния, аппроксимации Паде. 
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Введение. Исследование движения прямоугольного корпуса является актуальной задачей в мобильной 
робототехнической сфере [1]. Это движение происходит за счет управления материальной точкой внутри 
корпуса. Такие мобильные роботы могут использоваться при решении задач широкого класса. Например, при 
создании автономных роботов для исследования космического пространства и планет. В медицинских целях 
для диагностики и лечения, например, при прохождении по сложным структурам вен и артерий. А также при 
подземных работах и проведении исследований в тяжелых условиях, например, под водой и в местах больших 
перепадов температур [1,2]. 

Таким образом, всё новые задачи ставятся перед робототехникой, для решения которых требуются 
теоретические исследования, в т. ч. исследование моделей трения между корпусом и поверхностью в условиях 
комбинированной динамики [3, 4]. Так как движение мобильного робота происходит в разных направлениях, 
необходимо учитывать продольное движение и верчение. Таким образом, в структуре модели трения нужно 
обеспечить взаимосвязь между скоростями скольжения и верчения [5]. Важное развитие в описании данной 
взаимосвязи было сделано в работе [6]. Её автору удалось разрешить уравнения для главного момента и вектора 
сил трения, где в качестве площадки контакта рассматривался прямоугольник. Такие аналитические выражения 
позволяют наиболее точно описывать движение тел по шероховатым поверхностям, но они являются 
громоздкими и сложными, так как содержат в себе интегральные выражения. Поэтому авторы работы [7] 
построили дробно-линейные аппроксимации Паде, которые позволили находить решения для получаемых 
зависимостей. 

С помощью аппроксимации Паде можно объяснять действия комбинированного сухого трения для 
линейных и угловых скоростей. На основе приближений Паде появилась возможность создавать новые модели 
трения [8, 9], которые впоследствии начали классифицировать для лучшей интерпретации [10]. Классификация 
происходит в зависимости от числа параметров. Так, в работе [11] авторы ввели понятия размерности и порядка 
модели сухого трения в зависимости от порядка используемых Паде-аппроксимаций. 

Модель трения скольжения и верчения, которая предлагается в статье [12] дает возможность учитывать 
динамическую связь компонентов, которые определяют силовое взаимодействие прямоугольного корпуса и 
горизонтальной поверхности [13]. 

Постановка задачи. Рассматривается твердое тело массой т’, которое представляет собой 
прямоугольный корпус с однородными гранями длиной а, шириной В и высотой 2й. Вводится неподвижная 
система координат Оху2, связанная с корпусом (рис 1). Точка О находится на горизонтальной плоскости. 


Система Оху2 имеет начало в точке О, которая соответствует геометрическому центру корпуса. Ось О2, 
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параллельна оси О2. Ось Ох параллельна длинному ребру корпуса. Введем единичные векторы е,,е, осей 


О 


у и Ох соответственно. 





Рис. 1. Системы координат 


Рассмотрим безотрывное движение корпуса по поверхности (рис. 2), которое состоит из 
поступательного перемещения и вращения относительно оси О2. Три координаты определяют положение 
корпуса. Координаты х,, у, и Й задают начало системы координат Оху2 в координатах Оху2. Поворот 


корпуса относительно его исходного положения по оси Ох задается углом ф. В настоящей статье 


рассматривается случай, когда центр масс корпуса С’ и центр масс системы О совпадают (рис. 2) [14]. 





Рис. 2. Движение корпуса системы 


Материалы и методы. Область контакта представляет собой прямоугольник со сторонами а и, в 


котором нормальное напряжение зависит от расстояния от точки Р до граней прямоугольника (рис. 3). 





Рис. 3. Скорости точек О и Р 


Рассмотрим бесконечно малую площадь @5 в произвольно выбранной точке М на поверхности 


контакта. Введем угол ф между относительной скоростью скольжения и осью Ох. Проведем радиус-вектор 


’ ИЗ точки Р к точке М. Вектор скорости точки М обозначим чи, для его нахождения используем формулу 


Эйлера, описывающую распределение скоростей в абсолютно твердом теле: 
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м = ЭР + @.Тур. 
Скорость скольжения в точке М раскладывается на две составляющие вдоль осей 01х/ и О1 Ул: 
Фмх = 9х — У, 
Эму = у + ХФ. 
Используя закон Кулона, найдем малое приращение силы трения, направленной против относительной 
скорости в точке М [15]: 
аЁ = —Го(х,у) т, 


где / — коэффициент трения; о(х,у) — функция распределения контактных напряжений, зависящая от 
координат х иу; 45 = ахАу — малое приращение площади [15]. 
Перепишем дифференциал силы трения и момент этой силы в проекциях на рассматриваемые оси: 
= омх 
АЕ; = —Хо(х,у)  тахау; 


АЕ, = —Хо(х, у) МУ хау; 
|ом| 





р К 
ам, =|х У 09| = хаЕ, - уаЕ.. 
а, аЕ, 0 
В качестве частного случая рассмотрим равномерное распределение напряжения при отсутствии 
внутренних масс в корпусе, тогда эти напряжения будут равны: д = =. но далее запись продолжим в общем 
виде: о(х,у). 
Проинтегрировав выражения для сил трения, получим: 
а/2 гЬ/2 1; — Ум 
Е и о. Рю с(х,у)- т ахау; (1) 
а/2 гЬ/2 0,+х@ 
Е, — —/ Ра р 0(х, у) ти ахау. (2) 


Модуль относительного скольжения [ом | вычисляется по формуле: 


[бя | = о + уму = |5 +0 + 02 (х2 + у?) + 2%(5ух — охУ). (3) 


Представим взаиморасположение векторов переменной скорости скольжения \ и составляющих силы 
трения: Е — составляющая, направленная в противоположном направлении относительно скорости 
скольжения 1; ЁР, — составляющая, перпендикулярная мгновенной скорости проскальзывания. Одновременно 
представим координацию этой системы относительно осей О1х. и О1 у. (рис. 4). 





Рис. 4. Составляющие силы трения и скорости 


Осуществим переход от проекций скорости скольжения: 
ух = 0С05Ф, 
оу = о5ШФ, (4) 
к модулю скорости и углу скольжения: 
Е = Р;со$ф + Ё,5шФ, (5) 
Е, = Е, (—-$шф) + Е, 5шф 
Проинтегрируем момент силы трения на площадке контакта: 


М обе усозф)ни (2+2) 
М, = 78 Га У о(х, у) 02+62(х2+у2)+200(5шф-усо5Ф) ахау. 


Подставим в систему (5) выражения (1)—(3), а также перепишем выражение для момента силы. В 
результате получим трехмерную модель трения скольжения и верчения: 


5(с0$2 ф+зт?ф)-@(усо$ф+х$ф) 


а/2 Пы 


Ви= У Ла ь/2 9 У) хау; (6) 


02 +02 (х2+у2)+20(хзтф-усо$Фф) 
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—_ а/2 гЬ/2 , @(узшф+хсо$Ф) | 
7 Г а/2 Гы 0(х, у) 92 +02 (х2+у2)+2050(хзшф-усо$ф) оо (7 
р [4/2 [5/2 , о(хзшф-усо5ф)+(х? +2) 
— Г- а/2 Ь/2 (ху) 92 +02 (х2+у2)+2%0(хзтф-усо$ф) (8) 


Для того, чтобы не решать громоздкие интегралы, используют замену соответствующими 
разложениями Паде [16, 17]. Таким образом, исходя из теории Паде [18], эти выражения могут 
формулироваться как отношение двух функций нескольких переменных во всей области определения при 
условии, что функции должны иметь одинаковый порядок [7]. Для определения этих функций необходимо 
определить поведение интегральных выражений (6)—(8) при следующих условиях: 





аи  _ _ Л [4/2 [5/2 ‚ УХУ? (х’+у”)+(усозф+хф) (езтф-9с0$Ф) И 
0% =0 а -а/2 Гр ( ) (х2+у2)2 ахау — =; 
| 
дМ»2 —_ _ Л [а/2 (5/2 ‚ (2У?хзтф-2у3с05ф) а 
д\ 1№0=0 — «*-а/2 Тр о( ) (х2+у2)? ахау — =; 
а/2 гЬ/2 
М езвы = ГЛ а/2 Г ь/2 9 о(х ‚У). хе захау = — Те, 


а/2 гЬ/2 у$шф+хсо$ф 
ВА — ет ХУ = 1; 








0% ш=0о д =о 90 ю=о Пос 0% ш=0 0% =0 
с, дом, дом _ ь 
о до 1%=0 — до 1%=0 — до 10=0 — Мару — 





Значения выражений —— и 


поэтому их написание опущено в силу громоздкости. Тождественное равенство нулю реализуется при условии, 


не участвуют в нахождении последующих аппроксимаций Паде, 


когда напряжение о симметрично относительно центра прямоугольного пятна контакта, т. е. точки Р. 

Точная трехмерная интегральная модель [13] (6)-—(8) дает логичное описание явлений сухого трения, 
однако для решения задач динамики такая модель является сложно-приемлемой в силу необходимости 
вычисления внушительных интегралов [10]. Во избежание этой процедуры используется [6] замена точной 
интегральной системы на соответствующие выражения с помощью аппроксимаций Паде во всей области 
изменения переменных. Дробно-линейные разложения Паде дают трехмерную модель трения скольжения и 


верчения первого порядка [19]: 
о-+Ь1 4. 











Н =Ю а (9) 
_ -+Ь20. 

М» вы Мо «+420? (10) 
© +Бзо 

а (11) 


Для определения коэффициентов Паде необходимо изучить свойства этой модели в граничных точках 
по аналогии с интегральными выражениями. Для этого продифференцируем параметры Ё|, Ё,, М» и тем самым 
удовлетворим соответствующим интегральным выражениям: 


Н = и —— 
+1 о 


М» 2 = 16 
(09) 


Е = 5 
от 
9 


Система уравнений движения имеет вид: 
40, 


. 1 а 
ох 
(ть + т.) Г = В, + @то + т. )5уФо,; (12) 





— М; 


(то + т!) —— Е, — (то + т: 0х 020. . 


и = 
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Выразим производные по времени от скорости скольжения и угловой скорости с помощью формул (3)— 


(5) 
Фо _ 1 > 9х ро фо, 
а 7 02 + 02 ав У а 
а 1 0 $ 0 С0$ 
= — = (-- = г = В). 
а тт: 02 02 


Перепишем эти уравнения с помощью формул (3)-(5) при «Фо, = @/а и добавим первое уравнение из 
системы (12): 


40 
ар 7 М4 
а 
(тт + т.) в = Рр (13) 


(то + т. оф = Е. 
Результаты исследования. Далее вычислим интегральные выражения параметров К, И, [5, [6, 9 с 
помощью программного комплекса У\/оШгат Мафетайса при следующих значениях: 
т 
= 1 а= 0,5м; Б = 0,2 м; ту = 1 кг, а = ОЗ, = 87 КГ/ > 
и подставим в систему уравнений (13). На основе численных выражений построим графики интегральных и 


нормированных функций, зависимых от параметра К = У/‹› . На рис.5 приведены графики функций 
интегральных моделей трения (11)—(13), а также моделей на основе аппроксимаций Паде (14)—(16). 
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Рис. 5. Графики интегральных (сплошные линии) и нормированных (пунктирные линии) функций касательной (а), 
нормальной (6) составляющих силы трения и момента трения (в) 


Исходя из графиков функций (рис. 5), можно говорить о хорошем соответствии рассмотренных 
моделей. Далее получим графики зависимостей характерных параметров от времени (рис. 6). 
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Рис. 6. Зависимости скорости скольжения У (а) и угловой скорости и_(6) от времени Ё 


Обсуждение и заключения. Выполнено описание движения мобильного робота, начиная с контакта 
его корпуса и опорной плоскости, в условиях комбинированной динамики, когда присутствует скольжение и 
верчение. Получены аналитические интегральные выражения для касательной и нормальной составляющих 
силы трения [19] и момента трения применительно к прямоугольной площадке контакта. Для полученных 
выражений определены соответствующие аппроксимации Паде. Проведено сравнение интегральной и 
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нормированной моделей посредством построения графиков зависимостей сил трения и момента трения от 
угловой скорости и скорости проскальзывания. Результаты сравнения показали хорошее соответствие 
интегральной модели и модели на основе разложений Паде. Графики соответствуют логичному поведению при 
движении прямоугольного корпуса, т.к. скорость скольжения и угловая скорость увеличиваются 
соответственно заданным параметрам. Следовательно, комбинированная модель трения, реализованная с 
помощью аппроксимаций Паде, может быть применена для решения задач, связанных с мобильными 
подвижными роботами с прямоугольным основанием. 
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Задача о толстостенной сферической оболочке 


М.А. Артемов', Е.С. Барановский', А. А. Верлин', Э. В. Сёмка? 
' ФГБОУ ВО «Воронежский государственный университет» (г. Воронеж, Российская Федерация) 





Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора 
Н. Е. Жуковского и Ю. А. Гагарина» (г. Воронеж, Российская Федерация) 


Введение. Цилиндрические и сферические оболочки широко используются в технике. Они подвергаются 
внутреннему и/или внешнему давлению и тепловому воздействию. Распределение напряжений и деформаций в 
упругопластических оболочках изучалось многими исследователями. Большинство работ связано с 
использованием условий пластичности Мизеса, максимального касательного, максимального приведённого 
напряжения. Эти условия не учитывают зависимость от первого инварианта тензора напряжений и знака 
третьего инварианта девиатора напряжений. Для тел со сферической и цилиндрической симметрией при 
осесимметричном тепловом и силовом воздействии в ряде случаев удается получить численно-аналитические 
решения для напряжений, перемещений и деформаций. 

Материалы и методы. Решение задачи о состоянии толстостенной упругопластической оболочки проводится в 
рамках теории малых деформаций. Предложено условие пластичности, учитывающее зависимость от трех 
инвариантов тензора напряжений, а также знак третьего инварианта девиатора напряжений и трансляционное 
упрочнение материала. Решается несвязная термоупругопластическая задача. Для оценки напряжений в 
области упругого состояния сферической оболочки вводится эквивалентное напряжение, равное выбираемой 
функции пластичности. В качестве метода верификации напряженного состояния используется построение 
годографа вектора напряжений. 

Результаты исследования. Для линейных функций пластичности задача имеет аналитическое решение. 
Получено решение, учитывающее упрочнение материала. Определены аналитические и графические 
зависимости между параметрами внешнего воздействия для упругого и упругопластического состояния шара. В 
случае комбинированной нагрузки возможны варианты, когда пластическая область зарождается на 
внутренней, внешней границах шара или между этими границами. 

Обсуждение и заключения. Результаты вычислений показали, что учет пластической сжимаемости и 
зависимости предела пластичности от температуры может оказать существенное влияние на напряженное и 
деформированное состояние полого шара. При этом учет первого инварианта тензора напряжений в условии 
пластичности приводит к тому, что не только перепад давления между внешней и внутренней границами 
сферической оболочки, но и значения давлений на этих границах могут изменяться в ограниченном диапазоне. 
В данной постановке задачи, когда имеет место только тепловое воздействие, полый шар полностью не 
переходит в пластическое состояние. Результаты исследования позволяют прогнозировать поведение объекта 
(полого шара), испытывающего центрально-симметричные распределенные силовые и тепловые внешние 
воздействия. 


Ключевые слова: полый шар, толстостенная сферическая оболочка, термоупругопластическое состояние, 
эквивалентное напряжение, ассоциированный закон пластического деформирования, годограф напряжений, 
параметры управления поведением модели. 
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Гас{ фай по{ ошу Фе ргеззиге Агор Бебмееп фе ощег апа шпег Боцпдатез ое зрВешса1 зВей, Би Ве ргеззиге уа[чез 
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деГогтайоп 1а\у, згезз ВодоэтарВ, то4е] Бевахтог сопёго| рагатщег. 


Еог сйапоп: А. М. Ацетоху, Е. 5. ВагапоузКи, А. А. УегИа, её а1. Тыск-уаПеа зрБемса| зВе| рго ет. Адуапсеа 
Епошеенио КезеагсВ, 2021, уо1. 21, по. 1,р. 22—31. В рз://401.0г2/10.23947/2687-1653-2021-21-1-22-31 


Введение. Решение задачи о толстостенной сферической оболочке, испытывающей разные внешние 
воздействия, приводится в монографиях [1,2] и ряде научных статей по теории упругости, пластичности, 
термоупругопластичности [3—9]. Обычно рассматривается случай, когда процесс нагружения является простым. 

Задача о толстостенной сферической оболочке является одной из простейших упругопластических задач, 
когда поля внешних воздействий и внутренних параметров обладают центральной симметрией. В силу 
центральной симметрии в пластической области выполняется режим полной пластичности. Для идеального 
пластического тела задача статически определимая, что позволяет ее решить для любых условий пластичности. 
В монографии [1] приведено наиболее полное решение задачи о сферической оболочке, когда не учитывается 
зависимость предела пластичности от температуры, а условие пластичности не зависит от первого инварианта 
тензора напряжений и знака третьего инварианта девиатора напряжений. Рассматриваются случаи только 
теплового и комбинированного нагружения, когда задана температура на стенках сферы, давление на 
внутренние стенки, а на внешней стенке давление отсутствует. В работах [10—13] выполнено исследование 
термоупругопластического состояния разных объектов, а в [3-9] рассмотрен процесс теплового нагружения и 
разгрузки свободного от внешних усилий шара и полого шара для условия Треска с учетом зависимости 
предела пластичности от температуры. Решение этой и аналогичных задач представляет интерес, поскольку 
можно получить аналитическое или частично аналитическое решение для разных математических моделей. 
Аналитическое решение можно получить при выборе кусочно-линейных функций пластичности [11, 13]. 
Математическое моделирование объектов позволяет прогнозировать их состояние и поведение в зависимости 
от значений исходных параметров [ 14, 15]. 
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Материалы и методы. Постановка задачи. Рассматривается задача о толстостенной сферической 
оболочке (полом шаре), испытывающей центрально симметричные внешние воздействия: давление р, на 


внешнюю стенку при р=Ь и давление р, на внутреннюю стенку при р=а. Также рассматривается тепловое 
воздействие на шар: на границе р =а поддерживается температура Г, , на границе р=Б — температура Т,. 
Предполагается, что шар проявляет упругие и пластические свойства. Искомыми параметрами состояния в 
каждой точке шара являются компоненты тензора напряжений, компоненты тензоров деформаций и векторы 
перемещений. В области упругого состояния упругие деформации являются полными (остаточные деформации 
отсутствуют). 

Основные соотношения. Все соотношения приведены к безразмерному виду. В качестве масштаба 
длины выбирается внешний радиус шара Б. Все величины, имеющие размерность напряжений, отнесены к 
пределу пластичности при одноосном растяжении К. Масштабная единица для температуры — 1 С. 

В силу указанной симметрии внешних воздействий, в сферической системе координат р,0,ф матрицы 


компонент тензора напряжений и деформаций будут иметь вид: 


о 0 0 = 0 0 
(0)= 0 « 0, (=)= 0 в 0 
О 0 о ОО = 


ф ф 


В данном случае выполняются равенства 0, =0,, & =& 


ф? ф* 


Если функции пластичности не зависят от первого инварианта тензора напряжений и знака третьего 
инварианта девиатора напряжений, то при решении задачи о шаре функции пластичности будут приводиться к 
виду: 

Л=°,-©, ЕК. (1) 


Рассмотрим условие пластичности: 


(а, — 67)" +2(0, —85:)")" | 


о ,о,, 57,5 )= 
Л( р 0 р —. с+1(1+ о)" 


(2) 
. 1((, -с,-5(= — =’ |" +а(с, -в,-5(=’ —=)")" 
с+п + а)" 


— компоненты тензора пластических деформаций; ТГ — температура. 


= АСТ), 


Р 
0 


где 5, , & 
Когда параметры имеют значения: с =0, ц=1, 6 =0, а = 0, т=1 и=0,А=А, из условия (2) следует 
условие (1). На рис. | в плоскости о,,0, представлены кривые пластичности, определяемые по формуле (2) для 


разных значений числовых коэффициентов в функции пластичности. 


<] 
о, =—,ед 


| 


‚ед 


= | а 





Рис. 1. Кривые пластичности: а — для параметров: с=0,2; т=2;ч= 0,5; д=0; т=3; А =1; сплошная линия 


© = 0,5; пунктирная линия © = 0; 6 — для параметров: = 0; с =0;0=0; А =1 
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Результаты, представленные на рис. |, показывают, что при учете первого инварианта тензора 
напряжений радиальное и окружное напряжения, когда точка шара находится в упругом состоянии, могут 
изменяться в ограниченном диапазоне. Соответственно давление на границы шара также должно быть 
ограниченным. Когда первый инвариант в условии пластичности не учитывается, упругое состояние возможно 
для любого значения давления на границах шара, но ограничен перепад давления Др=р, -р,. Учет знака 
третьего инварианта девиатора напряжений, как отмечалось выше, влияет на значения пределов пластичности. 

Если значения параметров состояния ©,,0, определяют точку области, ограниченную кривой 


пластичности, принимается, что определяющими уравнениями, связывающими напряжения и деформации, 
являются соотношения закона Дюамеля-Неймана [1, 2]: 


Ев, = (1-у)0,-уб,+ЕаТ, Её, =в, -2ув, + ЕОТ, (3) 
где модуль Юнга Ё и коэффициент Пуассона У — константы. 
Если параметры состояния о,,с, определяют точки на кривой пластичности, то принимается аддитивное 
представление полных деформаций через обратимые и необратимые деформации: 
де Р зе р 
5 = +5, &=5 +5, . (4) 


Полные деформации определяются через перемещения по формулам: 





и Чи 
# ==, & =—. (5) 
`а 
р р 
Полные деформации связаны условием совместности деформаций: 
ЧЕ 
го +&-& =0. (6) 
р 
Чи 
Приращения необратимых деформаций связаны с напряжениями нормальным законом: 
ЧЕ? Ч? 
ое (7) 


д//до, 0//до,` 


Соотношение (7) при выборе нелинейных функций пластичности в общем случае неинтегрируемое [16]. 
В квазистатическом приближении напряжения должны удовлетворять уравнению равновесия: 





Чо’, 
р +2(6,-6,) =0. (8) 
ар 
Эквивалентное напряжение. Эквивалентное напряжение — это выпуклые изотропные скалярные 


функции тензора напряжений. В частных случаях термин «эквивалентное напряжение» является синонимом 
других терминов, например, «интенсивность напряжений» [17]. В настоящей работе эквивалентное напряжение 
совпадает с функцией пластичности. В этом случае на упругопластической границе эквивалентное напряжение 
не будет иметь разрыва. 

Поле температур. Поле температур в шаре находится из решения краевой задачи [ | |: 


а О 
ар ар (9) 


и, 





Р 


Решение задачи (9) представим в виде: 





г-т+-4АР В 1 АГАТ т. (10) 
(Б-а)\р 
Упругая область. В области упругого состояния полого шара формулы для напряжений имеют вид: 
В Ест 
с. а, в Е Е 
р р 2р 2р И-у)(6-а) 


Пластическая область. Выберем условия (1). Рассмотрим случай только теплового воздействия (10). 


Тогда пластическая область будет зарождаться на внутренней границе оболочки при условии [1]: $ 
в а дв _ 249 чай г 

(1—у) Ба+2Ь) з 

Обозначим через с — радиус упругопластической границы р=с. В процессе нагружения, когда 2 


В >В, пластическая область а <р<с, увеличивается. При выборе условия (1) напряжения в области а<р<с 


вычисляются по формулам: 25 
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о, =2кАШш(р/ а), ©,’ =’ +кА, к, =511(0,-0.), 
где к, =51(0, -в,)|.. .Если АТ >0,то к, =-1, если АТ <0,то к, =+1. 
Если область с <р<Ь остается упругой, то величины А, В и радиус упругопластической границы 
с, определяются из условий непрерывности напряжений на упругопластической границе и граничного условия 
о |= О. Так, если А и В определяются только из условий непрерывности напряжений на 


упругопластической границе, то имеют место выражения: 


. 
А=2к к ь [2+1 а ве ки: (11) 
Нл эл 56 33 





Уравнение для вычисления с, будет иметь вид: 


3 2 
2к^ ь (2+5 + а вв Л=0. (12) 
а) 3 3’ 3< 3’ Ь 


Если А и В определяются из условий непрерывности напряжений на упругопластической границе и условия 








с |, =0, то имеют место выражения: 


А= 











Е т м и, (13) 
3] 36 33 


Выбор формул (11) или (13) влияет на шаги алгоритма решения задачи, но не влияет на окончательные 
результаты. 
На границе р =ь будет зарождаться вторая пластическая область, если выполняется условие: 
(9,-6,) | ,=к,К, к, =—к,. (14) 
Для определения значения ДТ =ДТ, когда выполняется условие (14), необходимо решить совместно систему 
уравнений (12), (14). Поскольку параметр В в уравнения (13) и (1) входит линейно, то можно получить 
отдельно уравнение для определения радиуса упругопластической границы: 


дк А ОХе 9 —0 
а З(Б+с,)Бс, 


2 


(15) 


а также формулу для вычисления параметра В: 
в о 
° (5? —с”)аБ 


В процессе дальнейшего нагружения, когда выполняется неравенство В > В,, область шара с, <р<Ь 


(16) 


переходит в пластическое состояние на границе р=Ь. 
Полый шар при тепловом воздействии. Учет упрочнения материала. Рассмотрим случай, когда 
условие пластичности имеет вид: 
АО О Е, в Ле, = 56, =. ЛЕК, (17) 
Если до нагружения в шаре остаточные деформации отсутствуют, то в результате теплового нагрева 
пластическая зона будет зарождаться на внутренней границе р =а ‚ когда выполняются условия (15), (16). При 


дальнейшем нагружении образуется пластическая область а<р<с. Для нахождения напряжений в этой 


области надо из системы уравнений (3—6), (8), (17) получить соответствующие уравнения: 








‚ 4`в, 4, 6 баБЕоАТ 
р + 4р |= 0, 
Ар ар 1+30(1-уУ) р(Б-а) (18) 
рас, 
о, ==—+0.. 
2 ар 
Решение системы (18) запишем в виде: 
и ыы Уре сРвивАт |5 т 
+30(1-у —а 
р р (19) 





ИЕ И +2 Ш А + ы ВС 
1+361-у) 2р(Ь-а) 2р’ 
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Величины С,, С, ‚, входящие в формулы (19), определяются из граничного условия ©, | ,=0 и условий 
отсутствия пластических деформаций на упругопластической границе р=с: 


С =2к (К-К,)с' -ЗЕ а М, с», 


Я) — ° Е . 
С = кА па №), В№ 3“ р (20) 
: | За` а 


— барАТ 
° (+35а-у)6-а). 
В результате подстановки (20) в (19) получим: 
ЗбарЕалТ | 11 | _ 


1 
о ——_—_—_—_—__—_—_ НН ВНЕ 
` 1+360@-\) р-а а р 


2к Ас’ бабс,ЕаАТ 1 1 кет 1 
_ и |, 
3 р-а а р’ 3 Цар 


Ка ] | 


[кор 
а 





с, Е ВЕ кА+2кАшШР+ 
1+360(1-\) а р-а 


к Ас’ бабе ЕаАТ || 2 1 кА [2 1 
_ | ее 
3 2(Ь-а) а р’ за р 


Из решения упругой задачи следует: если ДТ > 0,то к =-1[; если ДТ <0,то к =+1. 


а 2р 





Учег пластической сжимаемости. Рассмотрим случай, когда функция пластичности является 
линейной относительно компонентов тензора напряжений: 
6(26. +6,)+к(5,-с,)+1(0, -о,)=&(-ХГ), 1 
к =51и(6, —С,). 
Условие пластичности (21) можно представить в виде: 
ас, +Во, =А(1-ХГ), 
&=2с+к-—\1, 
В=с-к-+1. 
С учетом введенных обозначений, для определения напряжений в пластической области получаем 
задачу: 
ос, +В, = (1-Х), 


Чо, 
Ро 620, (22) 


и 7 


Решение задачи (22) имеет вид: 


| (М -П 2МЕ |< ] г 1 2м | 
о =| р | № " + - К, 
Г " а-+В а(&-+2В) Др о+В р(а+28) 


В (М. -П  2М® а)“ [М- М, р 
= р | — + А, 
" (9) — о-в а(&-+2В) Др а+В  р(а+28))‘ 

адТ м аБАТ 
Ь-а’ р-а 


Для получения правильного результата из (20), когда, например, @+В=0, надо выполнить 





(23) 











где введены обозначения: М = Т, — ‚› М, =ХМ, М =ХМ. 


предельный переход при решении (22). Более просто получить правильный результат можно непосредственно в 
(23), при этом учесть, чтобы выполнялось условие @+[р}=0. В этом случае с = 0, поэтому получим: 


2 
с, 52% д-р [1-1 |2 | 
(9 а (ра 


с, А а-моанор+х, [1-2 =. 
(9) а р а 


В процессе нагружения пластическая зона зарождается на границе р =а ‚ когда выполняется условие: 
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_ @в+к-п(а- 25)БЕ а 
21-—у/(а? +аБ +?) } 
Результаты исследования. На рис.2 приведены графики напряжений и годографов вектора 


напряжений, когда область шара, соответствующая условию а <р < с ‚ находится в пластическом состоянии, а 


Ага -я в 


области шара, соответствующая условию с <р<ьЬ, — в упругом состоянии. 





ВЕТ 5 
а) 6) 
Рис. 2. Графики напряжений (а) и годографов вектора напряжений (6) 


при значениях параметров: А =; у = 0,3; а = 0,5; Б =ВАТ =170; с, =0,57 


На рис. 3 приведены графики напряжений и годографов вектора напряжений, когда области шара, 


соответствующие условиям а<р<с и с, <р<Ь, находятся в пластическом состоянии, а область шара, 


соответствующая условию с <р< с, , — в упругом состоянии. 
о {е] 
га , ед 9 = ра ь ед 





а) 0) 


Рис. 3. Графики напряжений (а) и годографов вектора напряжений (6) 


при значениях параметров: У = 0,3; а = 0,5; ДГ = 270; с, = 0,62; с, = 0,88 


На рис. 4 приведены графики напряжений и годографов вектора напряжений, когда область шара, 


соответствующая условию а<р<с ‚, находится в пластическом состоянии, а область шара, соответствующая 


условию с <р<Ь, — в упругом состоянии. 


Артемов М. А. и др. Задача о толстостенной сферической оболочке 





а) 0) 
Рис. 4. Графики напряжений (а) и годографов вектора напряжений (6) 
при значениях параметров: У = 0,3; а = 0,5; АТ = 215; с =0,58 
На рис. 5 приведены графики напряжений и годографов вектора напряжений, когда области шара, 
соответствующие условиям а<р<с и с, <р<Ь, находятся в пластическом состоянии, а область шара, 
соответствующая условию с <р<с, , — в упругом состоянии. 


7. ед о =—,ед 





а) 0) 


Рис. 5. Графики напряжений (а) и годографов вектора напряжений (6) при значениях параметров: 


у= 0,3; а = 0,5; АГ = 270; с, =0,61; с, =0,86 


На рис. 6 приведены графики напряжений и годографов вектора напряжений, когда области шара, 


соответствующие условиям а<р<с и с, <р<Ь, находятся в пластическом состоянии, а область шара, 


соответствующая условию с <р<с, — в упругом состоянии. 
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Рис. 6. Графики напряжений (а) и годографов вектора напряжений (6) 


при значениях параметров: У = 0,3; а = 0,5; АТ = 79; с = 0,1; ч = 0,1; Еа = 0,012; ух = 0,0017 


Обсуждение и заключения. Результаты вычислений показывают, что в данной постановке задачи, 
когда имеет место только тепловое воздействие, полый шар полностью не переходит в пластическое состояние 
(рис. 2—6). Упрочнение приводит к увеличению эквивалентного напряжения в пластической области и 
уменьшению радиуса упругопластической границы (рис. 4, 5). Упругая область в процессе нагружения не 
может полностью исчезнуть. Пластическая сжимаемость и зависимость предела пластичности от температуры 
оказывают существенное влияние на напряженное состояние полого шара (рис. 6). 
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Обоснование использования физической модели материальной точки 
при исследовании кинематических характеристик автомобильного 
транспортного средства при косом наезде на боковые тросовые ограждения 





Г. П. Колесникова 
ФГКВОУ ВО «Военная академия Ракетных войск стратегического назначения им. Петра Великого» (г. Балашиха, 
Российская Федерация) 


Введение. Выполнен обзор применения методов теоретической механики для разработки алгоритмов 
приближенного аналитического моделирования движения автомобильного транспортного средства (АТС) при косом 
наезде на боковые тросовые ограждения. Обосновано представление АТС как материальной точки при таком типе 
наезда. Цель работы — обосновать использование физической модели материальной точки для описания динамики 
АТС при косом наезде на боковые тросовые ограждения. 

Материалы и методы. Предложена новая физическая модель, описывающая противодействие движению АТС со 
стороны тросового ограждения при косом наезде; представлены новые методы приближенного аналитического 
построения перемещений АТС при косом наезде на боковые тросовые ограждения. Результаты аналитического 
расчета подтверждены данными конечно-элементного (КЭ) моделирования наезда по данным натурных испытаний. 
КЭ-моделирование проводилось с помощью многоцелевого конечно-элементного комплекса [$5 -Вупа. 

Результаты исследования. Разработаны новые алгоритмы для аналитического расчета перемещений АТС при косом 
наезде на боковые тросовые ограждения, а также новая физическая модель, описывающая противодействие 
движению АТС со стороны тросовых ограждений. Научно обосновано использование физической модели 
материальной точки для исследования динамики АТС при косом наезде на боковые тросовые ограждения, в том 
числе на основе сравнительного анализа кинематических результатов виртуального испытания с кинематическими 
расчетами, полученными на основе алгоритмов аналитического построения перемещений АТС. 

Обсуждение и заключение. Анализ результатов виртуального испытания в сравнении с аналитическими 
кинематическими расчетами показал, что представление АТС в качестве материальной точки при косом наезде на 
боковые тросовые ограждения является обоснованным, так как его движение при таком типе наезда близко к 
поступательному. Полученные результаты могут быть использованы при разработке и анализе корректности КЭ- 
моделирования бокового наезда АТС на тросовые ограждения. 


Ключевые слова: тросовое ограждение, аналитическое моделирование, траектория движения, перемещения, 
материальная точка, автомобильное транспортное средство, косой наезд. 
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Введение. Системы безопасности, предназначенные для организации транспортных потоков на дорожных 
трассах, способствуют минимизации неблагоприятных последствий дорожно-транспортных происшествий (ДТП). В 
настоящее время приобретают популярность тросовые дорожные ограждения, которые устанавливаются на 
разделительной полосе автомобильных дорог для предотвращения преднамеренных и непреднамеренных переездов 
транспортных средств через дорогу, а также на обочинах — для предотвращения выездов транспортных средств за 
пределы автомобильной дороги‘. 

Различные конструкции тросовых ограждений состоят из общих элементов: тросов, стоек, анкерных блоков, 
стяжных устройств, но разнятся по типу их монтажа” [1-4], в частности, по способу крепления тросов (рис. 1), что 
существенно влияет на формирование математических моделей. Также различно поведение стоек при наезде АТС. 
Стойки могут сминаться либо, выходя из гильз, компоноваться одна с другой и обеспечивать таким образом 
дополнительное торможение АТС [5]. 

Столкновение автомобиля с тросовым ограждением характеризуются сложной механикой взаимодействия, 
т.к. тросы обладают высокой степенью геометрической нелинейности, а стойки и грунт — высокой степенью 
физической нелинейности, при этом почти все процессы являются переходными [5]. 
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Рис. 1. Некоторые виды стоек тросовых ограждений [6] 


' Отраслевой дорожный методический документ ОДМ 218.6.017—2015. Методические рекомендации по применению дорожных ограждений 
различного типа на автомобильных дорогах федерального значения / рекомендован распоряжением Федерального дорожного агентства от 
23.12.2015 № 2489-р. ЧВГ: Вйрз://ез.зкоушЁЕ.ги/Оаа2/1/4293757/4293757596.рдЕ (дата обращения: 08.01.2021). 

`Ногпе РА. Верог 350 Ассерапсе оЁ Ме\и УогК ТЬгее-бгапа СаЫе Тегиатпа]. ОЁЯсе оЁ Н1оВууау баеу шйазгасваге, ЕНУГА, 9.3. Оерагитеп 
о Тгапзроцаноп — Уазтеюп, О.С.; Еебгаату 14, 2000. АуаЙаЫе Нот: Врз://61е\’ауз.401.20у/ (ассеззеа: 21.01.2019). 
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Для определения траектории АТС можно использовать уравнения колебаний струны, используемые в 
аналитических расчетах аэрофинишеров"” и вантовых конструкций” “. Но в связи со сложностью конструкции 
ограждений возникают сложности формулировки краевых условий. Энергетические методы также дают набор 
уравнений, которые не решаются аналитически ° [7]. Все это приводит к необходимости исследования 
ограждений с помощью инженерных пакетов прикладных программ как для моделирования самого тросового 
ограждения, так и системы краш-теста. При этом часто используется численный метод конечных элементов в 
явной постановке, реализуемый в программных комплексах [,5-2У МА, МАКС, МАЗТКАМ и др. [5, 8,9]. 

Недостатками построения КЭ-моделей являются: сложность разработки; длительность периода расчета 
(120 ч. и более); необходимость проверки правильности построения. Таким образом, период подготовки основы 
для исследования является довольно длительным. В связи с этим предложен метод приближенного 
аналитического расчета движения АТС по данным натурных испытаний, основанный на методах теоретической 
механики. Аналитический расчет позволяет проверить правильность разработанной КЭ-модели, а также 
сократить время исследования пассивной безопасности в 2-3 раза, поскольку для анализа пассивной 
безопасности уже нет необходимости строить КЭ-модель самого тросового ограждения, а использовать 
аналитически полученные перемещения. 

Материалы и методы. В результате наблюдения за АТС в процессе натурных краш-тестов была 
выдвинута гипотеза: при наезде на боковые тросовые ограждения АТС можно представить как твердое тело, 
совершающее поступательное движение (рис. 2). 


траектории точек АТС 


тросовое ограждение 





АТС 


Рис. 2. Модель АТС как твердого тела, совершающего поступательное движение 


Кроме того, делается ряд фундаментальных предположений о динамике АТС во время ДТП, 
обосновывающих использование предложенных теоретико-механических моделей явления: 

® боковые и продольные замедления транспортного средства постоянны в течение интервала времени, 
необходимого для ориентации АТС параллельно недеформированному барьеру; 

® вертикальные и вращательные ускорения автомобиля игнорируются; 


' Нуралиева А. Б. О динамике троса космического лифта : автореф. дис. ... к-та физ.-мат. наук. М., 2012. 20 с. 

` Михайлюк Д. С. Конечно-элементное моделирование и исследование динамики палубного аэрофинишера : автореф. дис. ... к-та техн. 
наук. СПб., 2009. 19 с. 

°Дядькин С. Н. Обоснование, технология навесной сборки и мониторинг байтовых пролетных строений мостов с учетом климатических 
факторов (на примере моста через реку Обь у г. Сургута) : автореф. дис. ... к-та техн. наук. Волгоград, 2005. 20 с. 

* Ле Тху Хыонг Оптимизация параметров пролетных строений висячих мостов, усиленных и не усиленных наклонными вантами, при их 
проектировании с применением ПК : автореф. дис. ... к-та техн. наук. М., 1999. 25 с. 

> Лось М. В. Численное моделирование поведения системы «тело-трос» с учетом изгибной жесткости троса и механизм петлеобразования : 
автореф. дисс. ... к-та физ.-мат. наук. М., 2000. 19 с. 
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® боковая составляющая скорости равна нулю после того, как транспортное средство перенаправлено 
параллельно перилам ограждения; 

® поскольку транспортное средство перенаправляется, оно не входит в зацепление с перилами 
ограждения; 

е деформация транспортного средства происходит в зоне удара, но положение его центра масс 
существенно не изменяется; 

® центр масс автомобиля движется так, как будто вся его масса сосредоточена в этой точке; 

® барьер может быть жестким или гибким; 

® силы трения шин автомобиля о поверхность проезжей части игнорируются; 

® система ограждений не содержит разрывов, которые могут привести к резким вертикальным 
движениям автомобиля. 

При таком представлении движения автомобильное транспортное средство (рис.3) возможно 
рассматривать как материальную точку некоторой массы т [10]. Начало отсчета декартовой системы координат 
соответствует точке начала контакта транспортного средства и ограждения. Движение материальной точки 
начинается со скоростью %, направленной под углом а к плоскости ограждения (оси х) [6]. 


у 





Рис. 3. Расчетная схема АТС [6] 


Противодействие движению транспортного средства со стороны  тросового ограждения 
моделировалось силами [6]: 


»е при боковом отклонение вдоль оси у — силой упругости Р, =-К и диссипативным 
д - . 
сопротивлением Р/" =-Бу, то есть Е, = —Ку-Бу; 
» сопротивление движению вдоль оси х — силой трения Г. =-—ЛМ\,, где № — сила инерционного 


давления, то есть №, = ту [6]. 
В результате движение АТС описывается системой дифференциальных уравнений [6, 11]: 
#=-/} 
9+257+р’у=0. 
Решением системы являются перемещения в функции времени #: 


х(1) = -№ (+, (с0з9+ Узша)!; 


У(г) = ТЯ зи ре. 


В [6] доказано, что сила трения Кулона со стороны дорожного ограждения не оказывает существенного 
влияния на моделируемое движение АТС. Физическое моделирование АТС как материальной точки не является 
новым при исследовании наезда на боковые ограждения. Такой подход был применен Р. М. Олсоном, 
Э.Р. Постом и Ф. Ф. Макфарландом при описании наезда автомобиля на жёсткие ограждения мостов [12]. 
Здесь расчет сопротивления бокового ограждения базировался только на классической модели трения 
Амонтона-Кулона. 

В отличие от этой задачи характер взаимодействия АТС и ограждения в настоящей работе описан 
принципиально новой моделью. При внешней похожести постановки задачи физическая сущность 
взаимодействия АТС и ограждения качественно отличается от сущности взаимодействия с более жесткими 
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ограждениями. В настоящей работе это было учтено посредством силовых моделей Амонтона-Кулона с учетом 
силы инерционного нормального давления и модели сопротивления Кельвина-Фойхта. 

Также, в контексте исследования характера наезда на тросовые ограждения следует отметить работу 
М. В. Вщетап и др. [9]. Данные натурных испытаний, представленных в этой работе, явственно 
демонстрируют, что при косом наезде на тросовые ограждения движение АТС близко к поступательному, когда 
угол рыскания не превышает 10°. Здесь процесс наезда описан двумя моделями: 

1. Простая динамическая модель транспортного средства, где использовался метод Рунге-Кутта 
четвертого порядка для численного решения дифференциального уравнения движения АТС. При этом 
учитываются силы воздействия со стороны ограждения (натяжение троса) и дороги (сила трения по 
классической модели Амонтона-Кулона). Во время удара полагается, что колеса не требуют рулевого 
управления. 

2. Квазистатическая модель ограждения, которая предназначена для расчета изменения 
напряженности тросов при изменении их геометрии в результате движения автомобиля. 

В результате расчетная модель получилась достаточно громоздкой и требует численного решения, 
тогда как модель, предложенная в настоящей работе, предусматривает алгоритмы расчета, легко проводимые 
аналитически без применений численных методов. 

Приближенный аналитический расчет кинематических характеристик АТС. В результате 
сделанных предположений получен алгоритм построения перемещения АТС, названный гармоническим [6]. 
Далее предложены квадратичный и кубический алгоритмы, основанные на построении полиномных функций 
перемещений [13]. Здесь противодействие движению АТС со стороны тросового ограждения вдоль оси у 
определяется некоторой функцией Р,(Е), в результате интегрирования которой функция перемещения АТС 
вдоль оси у имеет характер полинома второго, либо третьего порядка. Моделирование сопротивления 
движению вдоль оси х осталось прежним. 

Проведем сравнительный анализ результатов приближенного аналитического расчета кинематических 
характеристик АТС. На рис. 4—6 представлены расчетные данные для автобуса Мерседес-Бенц-0345. При 
построении перемещений АТС использовались данные натурных испытаний Государственного научного 
центра РФ ФГУП «НАМИ». 


у, м 





7.2 8.4 95 — 10.8 12 18 192 


Рис. 4. Взаимозависимость отклонений х, м и у, м, рассчитанных методами: кубического приближения — черная сплошная 
линия; квадратичного приближения — пунктирная линия; гармонического приближения — серая линия [13] 


Колесникова Г. П. Обоснование использования физической модели материальной точки при исследовании 
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Рис. 5. Зависимости отклонений х, миу, м, от времени [, с, рассчитанные методами квадратичного и гармонического 
приближения вдоль оси х — сплошная черная линия; вдоль оси у: пунктирная линия — методом квадратичного 
приближения, серая линия — методом гармонического приближения [13] 
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Рис. 6. Зависимости отклонений х, миу, м от времени [, с, рассчитанные методами кубического и гармонического 
приближения вдоль оси х — сплошная черная линия; вдоль оси у: пунктирная линия — методом кубического приближения, 
серая линия — методом гармонического приближения [13] 


Как видно из приведенных расчетных данных, расходимость аналитических методов между собой не 
превышает 20 %. Такие же вычисления были проведены и для автомобиля ГАЗ-3102 [12]. 

Построение конечно-элементной модели контакта тросового ограждения и автобуса Мерседес- 
Бенц-0345. Для анализа эффективности приближенного метода построения траектории движения АТС 
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разработана КЭ-модель взаимодействия тросового ограждения и автобуса Мерседес-Бенц-0345 в соответствии 
со стандартом организации СТО 521000-005-10690827-2015 - согласованным с Росавтодором в 2017 г. Согласно 
этому стандарту проводился монтаж тросовых ограждений для испытаний, результаты которых также были 
использованы при аналитическом моделировании траекторий и отклонений АТС. Объектом виртуальных 
испытаний были тросовые ограждения марки 14ДД/У4(300)-П-1, 1-3, 0-ГБ [11, 13] с параметрами: 

® гильзы восьмигранные марки ГЗ 500/У; 

® стойки СТД-2 состоят из двух элементов в виде квадратной трубы длиной 1500 м с размерами 
поперечного сечения 50х50Х3 мм; 

® шаг стоек с забетонированными гильзами на рабочем участке — 3 м; 

® шаг на начальном и конечном участках — 2 м; 

® высота стоек над поверхностью дорожного полотна — 1,1 м. 

Параметры троса: 

® трехпрядный, семижильный; 

® диаметр — 19 мм; 

® расчетное разрывное усилие = 173 кН; 

® количество ветвей — 4. 

Параметры автобуса Мерседес-Бенц-0345 [12, 13]: 

® полная масса — 14050 кг; 

® габаритная длина — 12000 мм; 

® габаритная ширина — 2500 мм. 

Скорость вхождения АТС в контакт с ограждением — 69 км/ч. Угол между плоскостью ограждения и 
направлением движения АТС — 20° в соответствии с методикой ГОСТ Р 52721—2007 (пп. 6.1-6.6) [11, 13]. 
Покадровые изображения АТС по данным натурного” и виртуального испытаний представлены в таблице 1. 


Таблица 1 


Сравнительные покадровые изображения АТС 


Момент Вид испытаний 
времени &, с натурные виртуальные 
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Результаты натурных испытаний практически совпали с результатами виртуальных испытаний, т. к. в 
обоих случаях динамический прогиб ограждения составил 1,5 м; длина пути взаимодействия АТС с 
ограждением по расчету — 23,5 м, по результатам натурных испытаний — 25 м. Относительная погрешность 
по длине пути взаимодействия АТС с ограждением не превысила 6 %. 

Результаты исследования. В результате виртуального испытания с использованием разработанной 
КЭ-модели тросового ограждения определены траектории характерных точек автобуса Мерседес-Бенц-0345. На 
рис. 7 отмечены точки АТС, для которых снимались показания, точка С соответствует центру тяжести. 


тросовое ограждение 








АТС 


Рис. 7. Выбранные для анализа точки АТС: А, В — на лобовой поверхности; 0, ЕЁ — на задней поверхности 
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Рис. 8. Зависимость перемещений характерных точек АТС от времени. Пунктирная линия — для точек А и В вдоль оси у; 
сплошные линии — вдоль оси х: черная — для точки 4, серая — для точки В 
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Рис. 9. Зависимость перемещений характерных точек АТС от времени. Пунктирная линия — для точек Д и Е вдоль оси у; 
сплошные линии — вдоль оси х: черная — для точки /), серая — для точки Ё 


Как видно из представленных данных, перемещения точек 4, Ви р), Е вдоль оси у попарно совпадают. 
Тоже можно наблюдать и для скоростей. Таким образом, при косом наезде АТС на боковые тросовые 
ограждения движение АТС действительно близко к поступательному. 

Обсуждение и заключения. Проведенный анализ кинематических результатов виртуального 
испытания в сравнении с аналитическими кинематическими расчетами показывает, что представление АТС в 
качестве материальной точки при его косом наезде на боковые тросовые ограждения является обоснованным, 
так как движение АТС близко к поступательному. Однако, при необходимости для более точного 
аналитического расчета кинематических характеристик движения точек, находящихся в областях, близких к 
точкам С, ОР, Е необходим дополнительный аналитический расчет, например, с использованием уравнений 
плоскопараллельного движения твердого тела. 
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Метод обратного анализа при математическом моделировании 
гидродинамического балласта в буровой установке 





И.Р. Антибас,, А. Г. Дьяченко, Б. И. Саед? 
'ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 
* Университет «Алеппо» (г. Алеппо, Сирийская Арабская Республика) 


Введение. При организации бурильных работ одна из основных проблем — точность и плавность опускания 
связок труб в шахту бурильной установки. Это зависит от многих факторов, в том числе от работы 
гидротормоза подъемного устройства. Цели данной работы: создание и исследование математической модели 
гидродинамического балласта в буровой установке. С помощью метода обратного анализа изучено влияние 
некоторых эксплуатационных показателей на тормозной момент гидравлического тормоза. 

Материалы и методы. Эксперименты проводились с использованием лабораторной установки, которая 
представляет собой модель гидродинамического тормоза. Его клапан закрывали при различных условиях, 
чтобы получить несколько значений давления с расчетом тормозного момента при подвешенном грузе 
определенной массы. Моделировались реальные (полевые) условия работы гидродинамического тормоза и 
сравнивались полученные результаты. При создании математической модели использован метод обратного 
анализа. Он основан на результатах экспериментальных измерений и позволяет выразить совокупность влияний 
отдельных переменных на тормозной момент. 

Результаты исследования. Создана и протестирована математическая модель гидравлического тормоза. 
Определена зависимость тормозного момента от давления, плотности и вязкости балластной жидкости. 
Влияние каждой переменной определено экспериментально, т. к. рассматриваемую зависимость нельзя 
представить в виде прямой связи. Методом обратного анализа получен набор значений констант, которые дают 
оптимальное решение. Принимая во внимание стандартный массив ошибок и минимальную стандартную 
ошибку, учтены статистические ошибки, допущенные во время экспериментальных измерений. 
Визуализирован физически приемлемый диапазон значений предложенной математической модели. С 
помощью базовой (нелинейной) математической модели рассчитан вспомогательный тормозной момент 
гидродинамического тормоза в зависимости от давления, плотности и вязкости. Установлена адекватность 
предложенной модели. В качестве критерия корректности использовались расчетные показатели тормозного 
момента. Ошибочное расхождение не превысило 6 %. Для дополнительного тестирования модели проведен 
расчетный эксперимент, имитирующий полевые условия. 

Обсуждение и заключения. Для математического моделирования гидродинамического балласта в буровой 
установке целесообразно использовать метод обратного анализа. Модель, предложенная в данной работе, 
связывает тормозной момент гидродинамического тормоза с рабочими параметрами жидкости внутри балласта: 
давление, вязкость и плотность. Доказана объективность модели. Предложена поправка к ней для симуляции 
работы тормоза в полевых условиях. Опираясь на полученные результаты, в будущих исследованиях 
целесообразно испытать созданную модель в полевых условиях с реальной полезной нагрузкой. 


Ключевые слова: гидромеханический балласт, математическое моделирование, метод обратного анализа. 
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Введение. При обустройстве буровых установок операции подъема и опускания бурильных труб 
обеспечиваются подъемными устройствами (катушками). Они снабжены тормозами двух типов, основной из 
которых — механический. В процессе бурения вес труб превышает полезную нагрузку. Для его компенсации 
вводят вспомогательный регулировочный балласт, предназначенный для снижения скорости при установке 
группы бурильных труб [1-4]. 

Балласты используют для контроля нагрузки и скорости движения крюка, а также для поглощения 
кинетической энергии группы скважинных бурильных труб. Помимо этого балласты: 

— снижают усилия на буровой установке, особенно при больших нагрузках; 

— уменьшают износ основных механических балластных элементов; 

— помогают медленно и плавно остановить груз, прикрепленный к крюку. 

Один из типов вспомогательных тормозов — гидродинамический. В этом случае находящаяся в 
балласте вода преобразует часть поглощенной кинетической энергии подъемной оси в тепловую, и за счет 
этого опускаются и поднимаются трубы [5-7]. 

При манипуляциях с трубами ось тормоза соединяется с осью подъемных устройств. Во время 
операции подвижная часть установки вытесняет воду в направлении наклонных лопастей внутри статора и 
вращается со скоростью, равной скорости вращения оси подъемных устройств. Вода попадает на ребра статора, 
а затем на его детали. Процесс повторяется, и увеличиваются силы, препятствующие движению ротора. Это 
создает тормозной момент, который отражает вращательное движение осей подъемных устройств и уменьшает 
скорость вращения. В результате скорость спуска группы бурильных труб уменьшается [8—10]. 


Антибас И. Р. и др. Метод обратного анализа при математическом моделировании гидродинамического балласта 





Коэффициент безопасности торможения (если его значение не превышает допустимого) координирует 
вспомогательную гидродинамическую работу и эффекты от работы двух основных тормозных систем, а также 
обеспечивает более длительное время обслуживания элементов главного тормоза. Это происходит за счет 
уменьшения износа фрикционных дисков и плоскостей тормозных колес. Увеличение тормозного момента 
балласта позволяет обеспечить правильное тормозное движение цилиндрического подъемника. Это 
свидетельствует о важности исследования гидродинамического тормоза. 

Гидравлическая мощность тормоза № определяется из соотношения: 
где р — плотность рабочей жидкости (пресной воды); г — ускорение свободного падения; О — количество 
израсходованной рабочей жидкости, равное объему рабочей жидкости, которая проходит через систему ребер 
за один рабочий цикл; Н — высота (уровень) рабочей жидкости в балласте. 

Тормозной момент Мв принудительного торможения двигателя определяется соотношением 

М, =р:=-О-Н/®, (2) 
где © — угловая скорость вращения движущегося участка балласта. 


Показатели кинематического торможения оребренных гидравлических машин определяются по 
известным теоретическим соотношениям тормозного момента: 


М, =^м -р-(2° -4°)-в>, (3) 
2 
МХ -в О в. 4 
вЕ=Аы °Р.( йо (4) 
Здесь р — наружный диаметр «кольца» рабочей жидкости, образованного при вращении колеса ротора и 


принятого равным диаметру колеса подвижной секции; 4 — внутренний диаметр кольца рабочей жидкости, 
который зависит от уровня балласта; п — число оборотов движущегося участка балласта (число вращающихся 
осей), об/мин: 
2пи , . 
о=———>® =—— 
60 100 

Коэффициент гидравлического тормозного момента Ам — безразмерная величина, учитывающая форму 
рабочих полостей в балласте, параметры и количество тормозных ребер. В практических расчетах среднее 
значение коэффициента гидравлического тормозного момента принимается 0,3. 

Цель исследования — изучить влияние некоторых эксплуатационных показателей на тормозной 
момент гидравлического тормоза, который представляет собой давление внутри тормозной камеры, а также на 
плотность и вязкость рабочей жидкости внутри балласта. Этих показателей нет в соотношении (4), поэтому 
примем, что тормоз работает при атмосферном давлении в | атмосферу, жидкость представляет собой пресную 
воду с плотностью 1 г/см” и вязкостью 1 стокс. 

Материалы и методы. Лабораторные эксперименты проводились в университете Алеппо (Сирийская 
Арабская Республика) на установке, представляющей собой модель ингибирования гидродинамического 
тормоза (рис. 1). 


- 























Рис. 1. Лабораторная установка — модель ингибирования гидродинамического тормоза: 1 — шкив, 2 — подлокотник 
шкива, 3 — манометр, 4 — выпускной клапан, 5 — термометр, 6 — отверстие в верхней части бака, 7 — регулирующие 
клапаны уровня жидкости, 8 — выходной канал для жидкости, 9 — выход жидкости из резервуара, 10 — клапан слива 
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жидкости, || — впускной клапан, 12 — входная линия, 13 — гидравлический тормоз, 14 — катушка, 15 — измеритель 
количества циклов, 16 — трос, 17 — полезная нагрузка, 18 — концентратор связи, 19 — линия управления, 20 — бак для 
жидкости 
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Этапы экспериментов 

1. Резервуар и балласт заполняют жидкостью до требуемого уровня (0,106 м до уровня первого 
клапана). 

2. С помощью ручного рычага, связанного с катушкой, на 0,317 м поднимают подвешенный на крюке 
груз. Он весит 8 кг и соединен с тросом. 

3. Позволяют связке труб упасть под собственным весом. 

4. Снимают показания с индикатора оборотов оси катушки. 

5. Затем фиксируют влияние показателей, полученных в ходе лабораторных экспериментов по 
изучению гидравлического торможения: 

— давление, создаваемое в тормозной камере (Р); 
— плотность рабочей жидкости (р); 
— вязкость рабочей жидкости (и). 

Сначала рассчитывают тормозной момент вспомогательного гидродинамического тормоза по 
соотношению (4). С этой целью определяются значения и, /), 4. 

Наружный диаметр «кольца» рабочей жидкости РО образуется при вращении колеса подвижной секции 
(ротора) и принимается равным его диаметру. Диаметр О движителя лабораторной установки в практических 
экспериментах равен 0,33 м. 

Внутренний диаметр кольца рабочей жидкости 4 зависит от уровня балласта. Для определения этого 
диаметра выравниваются объемы: 

— рабочей жидкости внутри балласта; 
— жидкости в «кольце», которое образуется при вращении подвижной секции внутри балласта и 
ограничивается высотой охлаждающей жидкости в баке (0,106 м). 

При ДР = 0,33 м рабочий объем балластной жидкости — 2,85 литра. Расчеты показали, что внутренний 
диаметр 4 = 0,283 м. 

Число оборотов и связано с подвешенным грузом. Если тормоз не работает, то это число оборотов 
намоточного устройства. Если же оно отделено от тормоза, то и берется из лабораторных экспериментов: при 
весе груза 8 кги = 240 об/мин. 

Рассчитаем момент гидравлического торможения для веса 8 кг, прикрепленного к крюку. Во время 
эксперимента бак был заполнен пресной водой до уровня первого клапана (0,106 м) при следующих 
показателях: давление Р = 1 атм, плотность жидкости р = 1 г/см° и вязкость жидкости ц = 1 Ст. Поэтому момент 
гидравлического торможения: 


2 
М, =0,3-1000-(0,33° - 0,283"). 27 =355 Н‘м. 
100 

Рассмотрим соотношение, обратное (2). В первую очередь речь идет об обратной зависимости между 
крутящим моментом и угловой скоростью (6). Уменьшение числа оборотов п на 20% при увеличении 
плотности жидкости означает увеличение крутящего момента на 20%. Учтем допустимое значение тормозного 
момента 355 Н`м. В итоге получим искомую величину тормозного момента при следующих показателях: 
давление Р = 1 атм, плотность р = | г/см” и вязкость жидкости ци = 1 Ст. 

Применительно к данной работе отметим следующее. При определении тормозного момента 
гидравлического тормоза учитывают число оборотов катушки, связанной с осью тормоза (то есть во время 
торможения): уменьшение числа оборотов намотки на определенный процент означает такое же увеличение 
крутящего момента. 

Влияние давления внутри балласта на тормозной момент гидравлического тормоза. Внутри 
гидравлического тормоза создали давление жидкости, закрыв клапан на линии выхода жидкости. Провели 
лабораторный эксперимент и расчет тормозного момента при подвешенном грузе массой 8 кг. Результаты 
представлены в таблице 1. 

Таблица 1 
Изменение тормозного момента гидравлического тормоза при изменении давления внутри балласта 


Г 
Левое граничное значение крутящего момента, Нм |105 | 14 | 103 | 100 | 9% | 9 





Правое граничное значение крутящего момента, Н*м 355 359 363 372 


На рис. 2 показано изменение тормозного момента при изменении давления внутри балласта. 
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Рис. 2. Изменение тормозного момента при разном давлении внутри балласта 


Влияние плотности рабочей жидкости на тормозной момент. В экспериментах использованы 
химические вещества, которые увеличивают плотность воды, не подвергая коррозии элементы лабораторного 
оборудования. 

Рассмотрены различные концентрации растворенных в воде веществ для получения балластной 
жидкости разной плотности. Условия те же: исследуется тормозной момент при подвешенном грузе 
8 кг (таблица 2). 

Таблица 2 
Зависимость тормозного момента от плотности жидкости, числа витков катушки 


и добавок в балластную жидкость 


[Чаю [| мо 
КЕ т ООО ООО 


Соли хлорида натрия 


Пищевой сахар 


И 





Отмечено очень слабое изменение вязкости при увеличении концентрации соли хлорида натрия. 
Можно считать, что вязкость жидкости приблизительно равна | Ст. 
На рис.2 показана зависимость тормозного момента от плотности рабочей жидкости при 


использовании соли хлорида натрия и пищевого сахара. 
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6) 
Рис. 3. Зависимость тормозного момента от плотности рабочей жидкости с введением 
соли хлорида натрия (а) и пищевого сахара (6) 


Влияние вязкости рабочей жидкости на торможение гидросистемы. Использовались различные 
комбинации концентраций растворенного в воде глицерина и силиката для получения балластных жидкостей, 
отличающихся по вязкости. Условия те же: исследуется тормозной момент при подвешенном грузе 
8 кг (таблица 3). 

Таблица 3 
Зависимость тормозного момента от вязкости жидкости и добавок в балластную жидкость 


ООО ПОЗ ПОС О ОИСИ 


Глицерин 


Силикат натрия 





Ь&р://уезбиК-доп$ва.га 


На рис. 4 показаны изменения тормозного момента с изменением вязкости рабочей жидкости при 
использовании глицерина и силиката натрия. 
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6) 
Рис. 4. Зависимость тормозного момента от вязкости рабочей жидкости внутри 
балласта с введением глицерина (а) и силиката натрия (6) 


Полученную таким образом почти линейную зависимость нужно проверить с помощью 
математической модели. 
В ходе исследования индекса вязкости установлено, что параметры плотности и вязкости мало зависят 
друг от друга и плотность не равна 1 г/см” (таблица 4). 
Таблица 4 
Соотношение плотности и вязкости жидкости 





Коэффициент влияния плотности не превышал 6 %. Его рассчитывали как разницу между плотностями 


2 < 2 
вязкой жидкости и пресной воды. Например, если вязкость была 30 Ст, то плотность — 1,025 г/см“. 
Следовательно, коэффициент воздействия плотности: 

1,025-—1 


В 


Результаты исследования 

Создание математической модели гидравлического тормоза. Зависимость тормозного момента от 
давления, плотности и вязкости нельзя представить как прямую связь. Влияние каждой переменной на 
тормозной момент определяется экспериментально (рис. 2-4). 
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Метод обратного анализа. В научной и справочной литературе нет зафиксированных показателей 
взаимосвязи тормозного момента с давлением, плотностью и вязкостью. Поэтому был задействован метод 
обратного анализа [4]. Он эффективен для создания математических моделей на основе экспериментальных 
измерений. С его помощью можно показать, как совокупность названных выше переменных влияет на 
тормозной момент. Это можно выразить соотношением: 


Р-Е(Р,С), (5) 


где А — функция, которая связывает изучаемое явление Ш) с набором значений переменных и набором значений 
для констант математической модели явления. 

Обратный анализ позволяет найти множество значений предполагаемых констант модели С через 
обратную зависимость: 


Р=Е'(2,,С), (6) 


где )„ — множество экспериментальных значений для изучаемого явления [4, 5]. 

Считается, что прямое решение для отношения (6) невозможно. Поэтому следует задействовать 
итерационную систему, чтобы получить набор значений констант, которые дают оптимальное решение. В 
обобщении указывается набор значений Р для соотношений: 


р. =Е(Р,С), (7) 
р, -Б.| 5, (8) 
где В. — набор рассчитанных значений для изучаемого явления; = — требуемая точность в соответствии с 


расчетом изучаемого явления. 

Предложенный метод учитывает статистические ошибки, допущенные во время экспериментальных 
измерений, ориентируясь на стандартный массив ошибок (5,) и минимальную стандартную ошибку. Таким 
образом, согласно статистическому распределению Гаусса зависимость плотности статистических данных: 


Л(Р)=Р = сопя1.ехр -5[(р. -5,) С’ (Б.-Б, )| 


(9) 
Здесь С, — стандартный массив ошибок. Символ { означает матрицу: 
У (1) он Оооо оьеоее 0 
| юны 5. (2) оны 0 а, (10) 
ооо ры» (п) 


где и — количество контрольных точек. 
О. является функцией Р, поэтому проблема связана с определением такого значения Р, чтобы О. дало 
максимальное значение Р.. 


5, =(Р.-Б,) С’ (Б.-Б,). (11) 
Минимальное значение 5. можно получить при более чем одном наборе значений Р. Поэтому 


правильный набор параметров формируется с помощью фактического физического моделирования значения 
этих параметров (рис. 5) [4, 5]. 





Рис. 5. Определение физически приемлемого диапазона значений (зеленая область): 
5 — физически приемлемое математическое решение, 55, 53 — математические решения 
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Поэтому мы используем набор первичных значений для параметров Ро, которые выводятся из средней 
нагрузки физического поля с учетом матрицы стандартных отклонений 50. 


При использовании нового распределения Гаусса получаем следующий статистический параметр 


интенсивности: 
Л, (Р)=Р =сопя1.ехр -[(Р-в) с; (Р-Р) | (12) 
Здесь С, — матрица стандартного отклонения: 
5, (1) Аааа. ее новы 0 
С, =| 0... ини, 5, (2) ааа О |. (13) 
О оао оная 5, (*) 


Проблема, выраженная в отношении (6), решается путем обратного анализа. При этом набор 
значений Р находится в пределах, обозначенных в (7) и (8). Таким образом, можно определить общую область 


параметров Р! и Р› по распределению Гаусса: 


Х(Р)=Р.В = сопз1.ехр(—5), (14) 
1 С —1 —1 
5=>[(Р.-Р.} ‘С/(Р.-Р.)+(Р-Р)-С7(Р-Р)|. (15) 
Итак, для решения нужно найти максимальные или минимальные значения 5 функции Г (Р). 
Единственный способ найти минимальное значение 5 — использовать цифровые методы, например 
Ньютона — Гаусса. Этот подход основан на преобразовании аналитических отношений в цифровые 


итеративные с учетом ошибки, допущенной из-за игнорирования некоторых ограничений в аналитических 


отношениях (рис. 6). 


Метод обратного анализа 


помсЕ множества значений 
констант и переменных 
на шаге частоты: 1+1 
Ва=В 








Расчет значений 
тормозного 
момента 
при наборе 
ен параметров Р.-1 


Мо 





Показ 


результатов 





Рис. 6. Алгоритм метода обратного анализа 


На следующем этапе работы с помощью базовой (нелинейной) математической модели вычислили 
вспомогательный тормозной момент гидродинамического тормоза. Он зависит от давления, плотности и 
вязкости: 


М» =Х(Р, р.и)-+С, (16) 
М, =а-р" +Ь-р"-+а-н"+С. (17) 
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Математическая модель (17) является общей нелинейной моделью, если значения п1, и2, Из отличны от 
единицы. Множество констант для поиска — это а, Б, а, а также и1, п>, Из. 
После применения итеративного цифрового метода с точностью, равной 0,001, найдены параметры, 
представленные в таблице 5. 
Таблица 5 
Значения параметров предлагаемой модели 


ПО ООО ООС ООС ООО ООС ООО ООО ООО 





Согласно таблице 5, исследуемое физическое явление можно представить как линейную модель (п1, И>, 


пз) = 1: 
М, =58,5Р +284, 6р +129 +10,62. (18) 

Математическая модель основана на результатах лабораторных экспериментов и подразумевает 
гомогенизацию единиц измерения в соответствии со значением констант. 

Адекватность математической модели. Первый этап определения адекватности модели — это 
испытание на торможение при: 

— давлении Р = 1 атм, 
— плотности жидкости р = 1 г/см, 
— вязкости жидкости и = 1 Ст. 

В данном случае торможение было равно начальному тормозному моменту Мв = 355 Н\м. Это логично, 
учитывая экспериментальные значения. 

Второй этап: три лабораторных эксперимента с подвешенным на крюке 8-килограммовым грузом. Они 
кратко описаны ниже. 

Первый. Взяли максимальные значения переменных и тормозного момента, а затем выполнили расчет с 
использованием математической модели по соотношению (18). Определили коэффициент ошибок, сравнив 
экспериментальный и математический результат. 

Второй. Взяли случайные значения переменных, относительно далекие от экспериментальных 
значений и от расчетных для тормозного момента (по математической модели). Вычислили коэффициент 
ошибок, сравнив экспериментальный и математический результат. 

Третий. Взяли случайные значения для переменных, относительно далекие от экспериментальных 
значений и от расчетных для тормозного момента (по математической модели). Вычислили процент ошибок, 
сравнив экспериментальный и математический результат. 

В таблице 6 приведены результаты подтверждения адекватности математической модели. 

Таблица 6 


Результаты проверки адекватности модели 


м д Вязкость Плотность | Тормозной момент Тормозной момент Ошибки, 
о авление 
жидкости жидкости (расчет) (математическая модель) % 
1 





Итак, результат подтверждает адекватность модели. Во-первых, уровень ее ошибочности оказался 
приемлемым. Во-вторых, выявленные ошибки имеют научное объяснение. Дело в том, что параметр плотности 
не зависит от вязкости, и это учтено в модели. Но эксперименты, проведенные ранее для подтверждения этой 
точки зрения, выявили между этими показателями корреляцию на уровне 6 % — и это примерно эквивалентно 
проценту ошибок, допущенных при использовании математической модели (в сравнении с расчетными 
данными). 

Для выражения совокупного влияния плотности и вязкости жидкости на работу гидравлического 
тормоза можно ввести в модель параметр и`р. Однако это необязательно, учитывая сравнительно низкий 
процент ошибок. В противном случае предложенная модель значительно усложнится. 

Имитация полевых условий. Предлагаемая математическая модель выражает физическое явление, 
выявленное и исследованное в лабораторных условиях при начальном тормозном моменте 355 Н\м. 

В полевых условиях балласты характеризуются начальным значением тормозного момента Мьо. 
Предлагается ввести в математическую модель параметр, который не изменит ее форму, но, предположительно, 
отразит «полевой» балласт: 


Антибас И. Р. и др. Метод обратного анализа при математическом моделировании гидродинамического балласта 





М, =58,5Р+284,6р+1,29и+10,62+(М„-М»„), (19) 


В 
где Мзв — тормозной момент «полевого» балласта после применения условий (давление, плотность и вязкость); 
Мво — начальный тормозной момент для полевого балласта; М»з„ — начальный тормозной момент 
лабораторного балласта, равный 355 Н\м. 

Это основная гипотеза. Она обосновывается следующим образом. Размеры лабораторного ингибитора 
выбраны не случайно, а так, чтобы соответствовать реальным меньшим размерам тормозов, производимых 
компанией Рагтас Г..Г.С (модель 112—500)". 

Обсуждение и заключения 

1. Математическая модель создана с помощью метода обратного анализа, который связывает 
тормозной момент гидродинамического тормоза с рабочими параметрами (давление, вязкость и плотность) 
жидкости внутри балласта. 

2. Эксперименты при различных случайных значениях переменных доказали адекватность модели. 
Сравнивались значения тормозного момента, определенные экспериментально и с помощью созданной модели. 
Ошибочность не превысила 6 % (таблица 6). 

3. Предложена поправка к модели для симуляции работы гидродинамического тормоза в полевых 
условиях. 

4. Опираясь на полученные результаты, в будущих исследованиях целесообразно испытать созданную 
модель в полевых условиях с реальной полезной нагрузкой. 


Библиографический список 

1. АБдочаеу, а. 5. Орйса| фотостарВу аз а РОЕ-сопзаште орипитаноп ргоет/ С. 5. АБдощаеу, 
К. Кеп, А. Н. Нле5сВег // Тпуегзе Ртоетз. — 2005. — Уо/. 21 (5). — Р. 1507-1530. РОГ: 10.1088/0266- 
5611/21/5/002 

2. Апа[уз1$ ОР Фе 0166$ затшр ег Юг Шегагсшса! шуегзе ргоетз/ 5. Асарлючц, /. М. Вагазеу, 
О. РаразрИороч[о$ [е{ а1.] // ЗТАМ/АЗА Лоштпа оп Опсецалу ОцапиЯсайоп. — 2013. — Уо1. 2 (1). — Р. 511- 
544. РОТ: 10.1137/130944229 

3. А теме\ оЁ зшгогае то4е|$ ап Фет аррПсайоп ю отоиид\ужег тодеПпо / М. 7. АзВег, В. Е. У. СтоКе, 
А. 7. лаКетап [© а|.]// \У/мег Кезомгсез Кезеатсв. — 2015. —[ Уо.. 518). —[ Р. 5957-5973. 
Ьрз://4о1.ого/10.1002/2015\/В 016967 

4. ЧатК, М. Р. Рагзишиае фе тео4 ог шире ууотКше Буроезез Гог Вудгоос1са1 шодепис / М. Р. Сагк, 
О. Кауе5к, Е. Еешсла//  \Умег Кезогсез Кезеасв. — 201. — У... 47 (9). — Р.1-6. 
6рз://401.0г2/10.1029/2010\В009827 

5. 5ей@ш, Е. ЕЕМ АррПе4 ю Нудгодупапис Веагие Оез1еп/ Е. Меат // Ме\м Тгфо|ое1са| У’ауз$. — 
ОптуегзИу Сатриз$ ЭТеР В: : шТесв, 2011. — Р. 451-476. БОТ: 10.5772/15525 

6. \Гапе, У. Соире4 Рупапис Апа[уз1$ оР Оееружег ОгШипе К1зег ипаег СотЫшеа Еогсше апа Рагатей1с 
ЕхсЦцаноп / У. Уапо, О. Сао, У. Еап® // Лоигиа| ог Мага баз Зслепсе апа Епошеегие. — 2015. — \Уо|. 27 (3). — 
Р. 1739—1747. 

7. Мао, Г. Оеер Умег ОптШие В15ег Месвап1са! Вебамюог Апа1у$15 Соп$14еги» Асба! К1зег Зите 
СопйЯсиганоп / Г.. Мао, О. Гла, 5. ХВоц [ей а1.] // Гоигпа[ ог Мага! аз эс1епсе апа Епошеегте. — 2016. — Уо1. 33. 
— Р. 240—254. 

8. \/апе, У. Орипитайоп Апа[уз1$ оЁ фе К15ег Тор Тепзюоп Еогсе ш Оеермужег ОгИп2: Апише аё Фе 
Мшипит Уапапсе ог Го\жег Нехе Лопе ОейЙесйоп Ап / У. \УМапе, О. Сао, Г. Еапе // Тоигпа| оГ Рего]еит 
эстепсе апа Епошеетие. — 2016. — У\о|. 146. — Р. 149—157. 

9. Свапе, У. шИиепна| Еас®юг$ Юг Фе Пез1еп оЁ ОШа-Оеер\ужмег ОиШие ЁК1зегз / У. Свапе, С. СВеп, 
Г. Хи [её а1.] // Рего]еит ЕхрогаНоп апа Оеуеортеп:. — 2009. — Уч]. 36 (4). — Р. 523—528. 

10. Сопе, \. Ап адариуе зитосже то4деппо-Базе затрПие зиажееу Юг рагатеег орипитаНоп ап 
ЧзитТочноп езитайой (АЗМО-РОПЕ) / \.. Сопе, О. Эиап // Епугойтета! МодеШие & ЗоЙ\ате. — 2017. — 
\Уо1. 95. — Р. 61—75. № рз://401.0го/10.1016/.епузой.2017.05.005 


Сдана в редакцию 04.12.2020 
Запланирована в номер 21.01.2021 


' Опоша! НУРВОМАТС® Втаке 620 251-5000 / Рагтлас Г..Г..С // ум. раттасбгаке.сот. 


Механика 


53 


ЬИр://уезК-допза.га 


А4уапсеа Епатеетиз Кезеагсй 2021. Т. 21, № 1. С. 43—54. 155М 2687—1653 





Об авторах: 

Антибас Имад Ризакалла, доцент кафедры «Основы конструирования машин» ФГБОУ ВО «Донской 
государственный технический университет» (344003, РФ, Ростов-на-Дону, пл. Гагарина, 1), кандидат 
технических наук, доцент, КезеагсвеШ: 0О-4789-2018, ОКСШ:  Б@р://огс14.0ге/0000-0002-8141-9529, 


ппа4.атураз(@та!.га. 


Дьяченко Алексей Геннадьевич, доцент кафедры «Основы конструирования машин» ФГБОУ ВО 
«Донской государственный технический университет» (344003, РФ, г. Ростов-на-Дону, пл. Гагарина, 1), 


кандидат технических наук, ОВСТО: В@р://огс14.0г2/0000-0001-9934-4193, @ехеу-а2(@лтай.га. 


Саед Бакир Имад, доцент кафедры «Сельскохозяйственное машиностроение» Университета Алеппо 
(САР, г. Алеппо, Университет Алеппо, механический факультет), доктор технических наук, доцент, ОВСТО: 


В рз://огс14.0г2/0000-0003-3855-7691, БаКи-111(@лтай.га. 


Заявленный вклад соавторов: 

И. Р. Антибас — научное руководство, постановка задачи, определение методологии исследования; 
сбор и анализ аналитических и практических материалов по теме исследования, критический анализ и 
доработка решения; компьютерная реализация решения задачи. А. Г. Дьяченко — анализ научных источников 
по теме исследования, критический анализ и доработка текста. Б. И. Саед — постановка задачи, определение 
методологии исследования; сбор и анализ аналитических и практических материалов по теме исследования. 


Все авторы прочитали и одобрили окончательный вариант рукописи. 


Айуапсей Епатеети? Кезеагсй 2021. Т. 21, № 1. С. 55-61. 155М 2687—1653 
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Моделирование гидравлической системы устройства с 
самоадаптацией по силовым и кинематическим параметрам 
на рабочем органе 





Т.А. Хиникадзе', А. Т. Рыбак’, П. И. Попиков” 

'ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 

`ФГБОУ ВО «Воронежский государственный лесотехнический университет имени Г. Ф. Морозова» (г. Воронеж, 
Российская Федерация) 

Введение. В настоящее время в России принят курс на создание интеллектуальных машин и оборудования. Это 
относится и к мобильным технологическим машинам дорожно-строительного и коммунального назначения. 
Поэтому проектирование и создание такого вида исполнительных устройств с функцией самоадаптации 
является актуальной задачей. 

Материалы и методы. Представлено устройство, оснащённое гидроприводом с самоадаптацией по нагрузке и 
согласованием кинематических и силовых параметров главного движения и движения подачи рабочего органа 
установки для бурения породы. Для исследования и проектирования устройства на основе методов 
математического моделирования гидропривода и адаптивных систем предложена математическая модель, 
разработанная с использованием основ теории объёмной жёсткости гидравлических систем. Это позволяет 
максимально точно описать влияние динамических свойств гидравлической системы (сжимаемость рабочей 
жидкости, упругие свойства трубопроводов, рукавов высокого давления, гидравлических аппаратов) на 
динамические свойства системы в целом. 

Результаты исследования. Математическая модель для устройства с самоадаптацией включает подмодели 
адаптивных связей, взаимосвязи силовых, кинематических и технологических параметров бурения породы, а 
также математическое описание движения элементов системы. Решение полученной математической модели 
производилось в среде программного обеспечения динамического моделирования технических систем 
ош тТесВ. В результате получены общие зависимости адаптивной системы от конструктивных параметров и 
технологических условий функционирования. 

Обсуждение и заключение. Математическая модель представленного устройства показывает принципиальную 
возможность реализации принципа самоадаптации по нагрузке при внешних и внутренних возмущающих 
воздействиях в процессе работы. Полученные результаты могут использоваться при проектировании 
адаптивных систем другого технологического оборудования, например, для глубокого сверления заготовок с 


переменными свойствами по глубине. 


Ключевые слова: гидравлическое устройство с самоадаптацией, функция гидропривода, обобщенная 


математическая модель, адаптивные связи, согласованные движения, рабочий орган, нагрузка, стабилизация. 
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Т. А. КыиКаадхе', А. Т. ВуБак', Р. 1. РорщЖоу” 
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Введение. Принятый в России курс на создание интеллектуальных машин является основой для 
проектирования и создания исполнительных устройств с функцией самоадаптации. Решение этой задачи 
является своевременным и актуальным. 

Свойством самоадаптации обладают приводы с дифференциальными связями элементов внутренней 
структуры [1]. К такому классу технических систем относится устройство для осуществления бурения породы с 
переменными свойствами (патент РФ № 2582691). Свойство самоадаптации реализуется устройством с 
обратной отрицательной и прямой положительной связями [2]. 

На качество процесса самоадаптации оказывают влияние внешние и внутренние воздействия — 
непостоянство нагрузки, сопротивление среды, сухое и вязкое трение, объемная жесткость жидкости и 
трубопроводов, адаптивные связи. 

Материалы и методы. При современном уровне развития вычислительной техники, сложности и 
высокой стоимости проектируемого оборудования широкое распространение получили методы 
математического моделирования на стадии разработки [3-13]. Разработка специальной модели для расчетов и 
вычислительного эксперимента по определению параметров на этапе проектирования дает возможность 
обоснованного выбора типоразмера устройства с самоадаптацией. 

Авторами разработана математическая модель для решения задачи проектирования устройства в среде 
динамического моделирования технических систем ЗииШшТесь (Зипчайоп Ш Тесбис) [14, 15]. Этот 
программный продукт позволяет проводить моделирование технологических процессов, протекающих в 
различных отраслях с одновременным моделированием систем управления. 
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На рис. представлено устройство для осуществления функции самоадаптации в условиях силового 
сопротивления на рабочем органе, которое содержит: насос постоянной производительности (Н), 
предохранительный клапан (КП), фильтр (Ф), регулируемые дроссели (Др1) и (Др2), регулятор расхода (РР), 
гидроцилиндр подачи (ГЦ) и гидромотор главного движения (ГМ), гидрораспределители (РТ), (Р2) и (РЗ), 
манометры (Мн1), (Мн2) и (Мн3), гидравлический бак (Б) и трубопроводы. 
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Рис. 1. Принципиальная схема устройства для осуществления функции самоадаптации в условиях силового сопротивления 
на рабочем органе 


Устройство предназначено для реализации функции самоадаптации по нагрузке (с целью её 
стабилизации) и согласования выходных движений рабочего органа технологических машин и оборудования, а 
также функций собственно силового привода. Поэтому в устройстве предусмотрены многопозиционные 
распределители для формирования потоков и направлений рабочей жидкости. К таким дополнительным 
функциям устройства относятся: разогрев масла в гидросистеме в начале её работы; «взвешивание» подвижных 
частей (холостой нагрузки) механизма подачи при наладке на технологический режим; ускоренный подъём 
(отвод) инструмента с вращением, но без регулирования его скоростей; подачу инструмента «на объект 
воздействия» или «от объекта воздействия» с вращением и регулированием скоростей инструмента. 

Математическая зависимость для внутренней отрицательной связи устройства устанавливается путем 
совместного решения уравнений, вытекающих из уравнений баланса давления по системе и уравнения 
баланса расхода рабочей жидкости в системе обратной связи. 

Из баланса давлений в системе вытекает 


1 1 Лошт 
АРлр5 = с Мм + р. вы (1 — пит) АРлрло, (1) 


где ДрдрБ и Друр1о — перепады давлений на дросселях 5 и 10; 4; м и М‚„ — приведенный объём гидромотора 
ГМ и создаваемый им крутящий момент; 1 и иг — площади поршня гидроцилиндра в поршневой и 
штоковой полостях; № — сила сопротивления со стороны грунта, преодолеваемая при перемещении штока 
гидроцилиндра ГЦ. 
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Уравнение баланса расходов в цепи обратной связи имеет вид 
Фдрло И др о а (2) 
где Олр5 и Одрло — расходы рабочей жидкости через дроссели Др! и Др2; и, — скорость перемещения 
поршня гидроцилиндра ГЦ. 

Имея в виду, что гидронасос Н, питающий систему, имеет постоянную производительность, расход 
рабочей жидкости через регулятор расхода РР в процессе функционирования системы привода остаётся 
постоянным. Становится очевидным, что расход рабочей жидкости через дроссель Др2 также будет 
постоянным. В этом случае зависимость ускорения движения поршня гидроцилиндра ГЦ от суммарной 
нагрузки на гидродвигателях главного движения и движения подачи примет вид 


Чт _ _ (1 @Мы , 14Е 
а — А (= а "> я (3) 





где и, — скорость перемещения поршня гидроцилиндра ГЦ; & — время; А — коэффициент обратной связи, 
зависящий от параметров номинального режима работы системы, конструктивных параметров гидромашин, 
дросселей и его настройки. 

Из уравнения (3) видно, что положительным приращениям крутящего момента на валу гидромотора 
и перемещению поршня подачи гидроцилиндра соответствует отрицательное приращение скорости 
перемещения подачи инструмента. Другими словами, при увеличении сопротивления со стороны грунта 
вращению режущего инструмента или перемещению его подачи, скорость подачи инструмента уменьшается. 

Прямая связь между скоростями и ускорениями вала гидромотора главного движения и штока 
гидроцилиндра движения подачи имеет вид 








1 
и = т (46мм = В 
п 
ИЛИ 
п Догмы 
= 4 4 
Е гм ЧЕ ( ) 
где ® == ловая скорость вращения вала гидромотора, — расход абочей жидкости через регулято 
ГМ рр У 
расхода. 


Из уравнения (4) видно, что при увеличении частоты вращения гидромотора главного движения 
(происходит при уменьшении сопротивления со стороны обрабатываемой поверхности) скорость перемещения 
поршня гидроцилиндра подачи возрастает. 

Математическое моделирование движений рабочих элементов устройства. Уравнения движения 
ротора гидромотора основного движения и поршня гидроцилиндра подачи инструмента имеют вид: 








атм __ 
] м ак ДР; тм — Мсопр› (5) 
Ап 
Та АЕ _ Г. Ри о: Ришт — Еопр, (6) 
где м — суммарный момент инерции всех вращающихся элементов системы, приведённый к валу 
гидромотора; т, — масса всех подвижных частей системы, приведённая к поршню гидроцилиндра; Метр — 


суммарный момент сопротивления вращению рабочего органа, приведённый к валу гидромотора; Копр — 


суммарная сила сопротивления перемещению инструмента со стороны обрабатываемой поверхности, 
приведённая к поршню гидроцилиндра; Ар, — перепад давления на гидромоторе; рии ри ит— давления в 
поршневой и в штоковой полостях гидроцилиндра соответственно. 

Моделирование свойств гидравлической системы 

Математическая модель гидравлической системы устройства для осуществления функции 
самоадаптации в условиях силового сопротивления на рабочем органе разработана с использованием объемной 
жесткости, которая даёт возможность осуществлять моделирование максимально приближенно к реальным 
характеристикам [16-20]. При моделировании особое внимание уделено определению приведенного 
коэффициента объёмной жёсткости рукавов высокого давления [16]. Полученная математическая модель 
включает следующие уравнения: 

—щ уравнение приращения давления в различных точках (на рис. 1 обозначены точками с 1 по 23) 
гидравлической системы имеет вид: 


ар — ба: о. вх: > › Оисхе ЧЕ ) (7) 
где » Оки У, Ос — суммы всех расходов рабочей жидкости, входящих и исходящих из рассматриваемого 
(1-го) объёма системы за время 4; Си»; — приведенный коэффициент объёмной жесткости выделенного участка 


гидравлической системы; 
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— уравнение для определения расходов рабочей жидкости через различные элементы гидравлической 


2 ь 
0; = ВУ | [р; — Р:+1| $181 (ф; — Рича), (8) 


где р; ир;:.1 — давление на входе и на выходе гидравлических сопротивлений; / — площадь живых сечений 


системы имеет вид: 


сопротивлений; р — плотность рабочей жидкости; 
— формула расчёта приведенного коэффициента расхода линейных сопротивлений следующая: 


ее (9) 
С 
1 
где Чи [ — диаметры и длины линейного участка трубопровода; А, — коэффициент гидравлического трения 
трубопровода; 
— формула для определения приведенного коэффициента объёмной жёсткости металлических 


участков трубопровода имеет вид: 








НЕ: ИЕ 
С —_ па21 1-Е (10) 
б Е 


где 4 и [ — внутренний диаметр и длина участка трубопровода; д — толщина его стенки; Ед и ЕЁ! — модули 
упругости жидкости и материала стенки гидролинии. 
Приведенные коэффициенты объёмной жёсткости РВД определяются экспериментально. 
Производительность гидронасоса определяется с учётом объёмных потерь по формуле: 


Чрн®н 
Од = “Пон (11) 


2п 
где рн — рабочий объем гидронасоса, ©, — скорость вращения вала гидронасоса, По,н —_ текущая величина 


объемного коэффициента полезного действия гидронасоса. 
Расход рабочей жидкости через гидромотор определяется по формуле: 


Чрм@м 
Оо — к. (12) 


21 0,м 
где (м — рабочий объем гидронасоса; о„ — скорость вращения вала гидронасоса; Пи„ — текущая величина 
объемного КПД гидромотора. 
Текущая величина объемных коэффициентов полезного действия гидронасоса и гидромотора 
определяются по формуле: 





И .._ РВ 
По = Е (1 — Поло ) (13) 
Рр.пот 
ГДе Попо —Щ Номинальная величина объёмного КИД гидронасоса и гидромотора; рп —- величина 


номинального давления гидравлических машин; р, — текущее значение давлений на насосе или моторе. 
Результаты исследования 

Предлагаемая математическая модель рассматриваемого устройства позволяет с высокой точностью произвести 

теоретические исследования эксплуатационных возможностей при бурении породы ещё на стадии 

проектирования. Расчёт системы бурения, произведённый с использованием программного обеспечения 

эштТесй [10, 11], показал правильность этого утверждения. 
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Рис. 2. Изменение параметров функционирования буровой машины при ступенчатом изменении сопротивления вращению 
основного рабочего органа буровой машины со стороны грунта: а — перемещение поршня гидроцилиндра (х„) и его 
скорости (у„); 6 — момент сопротивления вращению вала гидромотора М, 
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На рис. 2 представлены графики изменения параметров движения поршня гидроцилиндра подачи при 
ступенчатом изменении момента сопротивления вращению основного рабочего органа буровой машины со 
стороны грунта. 

Анализ результатов расчётов, представленных на рис. 2, показывает, что с увеличением момента 
сопротивления вращению рабочего органа буровой машины (рис. 2 6) скорость поступательного перемещения 
инструмента уменьшается и может принимать отрицательные значения (рис. 2 а), что говорит о выполнении 
системой функции адаптации. Свойства системы требуют дополнительного исследования. 

Заключение. Предлагаемая математическая модель позволяет ещё на стадии проектирования делать 
предварительную оценку качества работы и выбирать наиболее рациональные параметры системы бурения, 
оснащённой адаптивным гидромеханическим приводом, при различных условиях эксплуатации и варьировании 
её конструктивных свойств. Полученные результаты и их анализ говорит о том, что предложенный способ 
адаптации основного и вспомогательного движений системы привода может с успехом применяться и в другом 
технологическом оборудовании, например, при глубоком сверлении многослойных металлических заготовок. 
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Введение. Рассмотрена дуговая наплавка путем подачи дополнительной присадочной проволоки, подогреваемой до- 
полнительной дугой, горящей между присадочной и электродной проволокой. В условиях такой наплавки исследована 
минимизация поступления металла переплавляемой подложки в наплавляемый металл. Цели работы: изучение усло- 
вий, обеспечивающих саморегулирование данного перспективного процесса дуговой наплавки, а также оценка воз- 
можности управления мощностью теплового воздействия на металл и на поток наплавляемого металла. 

Материалы и методы. При решении широкого круга задач сварки и наплавки целесообразно использовать методы 
инженерного анализа, основанные на физико-математическом моделировании процессов и явлений. К их числу отно- 
сятся: 

— саморегулирование дугового процесса при совместном плавлении электродной и присадочной проволок; 

— оценка возможности управления мощностью теплового воздействия на металл и на поток наплавляемого металла 
при формировании наплавочной ванны. 

В работе рассмотрены особенности дуговой наплавки антикоррозионных хромоникелевых сталей на низколегирован- 
ную сталь. 

Результаты исследования. Предложены новые математические зависимости, описывающие физические явления при 
наплавке с дуговым взаимодействием между электродной и присадочной проволокой. Разработана физико- 
математическая модель совместного плавления электродной и присадочной проволок. Она позволяет определять зна- 
чения регулирующих параметров. Кроме того, можно узнать, с какой мощностью на подложку воздействует тепло от: 
— тепловыделения в основной дуге, 

— потоков капель наплавляемого электродного и присадочного металла, 

— излучения плазмы дуг. 

Обсуждение и заключения. Установлено, каким образом на ток и длины основной и дополнительной дуг воздействуют 
напряжения питания. Определена скорость подачи электродной и присадочной проволок диаметром 1,6 мм и 1,2 мм из 
сплава шсопе] 625. Показано, какое тепловое воздействие в этом случае испытывает подложка. Отмечено, что ввиду 
большего значения тока основной дуги диаметр электродной проволоки должен быть больше, чем у присадочной. Теп- 
ловой поток в подложке создается преимущественно потоком капель наплавляемого металла. 


Ключевые слова: физико-математическая модель, дуговое взаимодействие, наплавка, плавящийся электрод, присадоч- 
ная проволока. 
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Введение. Наплавка материалами, характеристики которых заметно отличаются от основного металла, 
обеспечивает необходимые эксплуатационные свойства контактных поверхностей изделий. При наплавке следует 
уменьшить перемешивание материалов, поэтому нужно регулировать мощность теплового воздействия на 
поверхность основного металла и количество наплавляемого материала. В работе [1] впервые была обоснована 
возможность такого регулирования путем изменения теплового потока при дуговом взаимодействии плавящегося 
электрода с дополнительной присадочной проволокой. Первоначально данный подход использовался в основном 
при сварке [2-4], а в последние годы — и при наплавке [5]. По мнению автора работы [6], этот способ наплавки 
поможет обеспечить высокие эксплуатационные свойства уплотнительных и контактных поверхностей 
трубопроводной арматуры. Однако ряд исследователей отметили фактор, затрудняющий использование данного 
метода при сварке и наплавке. Основной проблемой они считают регулирование подачи присадочной проволоки. 
При недостаточной скорости подачи она может плавиться вне наплавочной ванны, а при высокой — не 
расплавляться. Поэтому для обеспечения требуемого потока наплавляемого материала при ограниченном тепловом 
воздействии дуги на поверхность металла необходимо точно задавать ток дополнительной дуги, горящей между 
проволоками. Условия выполнения работ, нестабильность подачи проволок, неоднородность их теплофизических 
свойств также непредсказуемо меняют условия плавления основного электрода и присадочной проволоки. Поэтому 
стабильность процесса наплавки можно обеспечить, используя эффект саморегулирования основной и 
вспомогательной дуг. Цели работы: исследование условий, обеспечивающих саморегулирование дугового процесса 
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при совместном плавлении электродной и присадочной проволок, а также оценка возможности управления 
мощностью теплового воздействия на металл и на поток наплавляемого металла. 

Материалы и методы. При решении ряда задач перспективно использование методов инженерного 
анализа, основывающихся на физико-математическом моделировании процессов сварки [7, 8]. Такие модели 
представляют собой систему дифференциальных уравнений, граничные условия которых учитывают множество 
технологических факторов. Уравнения решаются в итерационном цикле. Важное преимущество такого метода 
инженерного анализа на основе фундаментальных физических законов —Й универсальность результатов и 
возможность их использования для изучения механизмов физических взаимодействий в процессах сварки и 
наплавки [9]. К таким явлениям относятся: 

— саморегулирование дугового процесса при совместном плавлении электродной и присадочной проволок, 

— оценка возможности управления мощностью теплового воздействия на металл. 

Однако с учетом эффекта саморегулирования дуг необходимо дополнительно исследовать особенности 
совместного плавления электродной и присадочной проволок. Кроме того, нужно узнать, с какой мощностью на 
подложку воздействует тепло от: 

— тепловыделения в основной дуге, 

— потоков капель наплавляемого электродного и присадочного металла, 

— излучения плазмы дуг. 

В рамках данной работы рассматривается двухэлектродная сварка (англ. доц е-е]еско4е газ те агс 
мет», РЕ-СОМАУ/) и наплавка электродами (плавящийся электрод и токоведущая присадочная проволока) в за- 
щитном газе. 

Основная часть. Для изучения механизма взаимодействия основной и дополнительной дуги между 
электродной и присадочной проволоками необходима физико-математическая модель, учитывающая: 

— особенности плавления электродной и присадочной проволок, 

— эффект саморегулирования тепловой мощности дуг при наплавке. 

Численное решение уравнений модели в итерационном цикле позволяет более глубоко изучить 
регулировочные характеристики процесса и их влияние на тепловую мощность наплавки. 

Физическая модель плавления электродной и присадочной проволок. С учетом многих значимых факторов 
при исследовании особенностей саморегулирования плавления электродной и присадочной проволок целесообразно 
использовать подходы, изложенные в [10]. Однако пространство моделирования нуждается в уточнении. 
Рассмотрим уточненное пространство моделирования условий плавления электродной и присадочной проволок под 
действием дуг. Первая из них — а (рис. 1) горит между электродной проволокой | и подложкой 3, а вторая В — 
между электродом и присадочной проволокой 2. При этом, по нашему мнению, проволоку целесообразно подавать 
со стороны передней части наплавочной ванны. Это позволит исключить приварку («замораживание») 
нерасплавившейся присадочной проволоки в застывающем металле хвостовой части сварочной ванны. 





Рис. 1. Схема процесса двухдуговой наплавки плавящимся электродом с дополнительной присадочной проволокой 
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Для формирования общего плазменного столба дуг а и Ь присадочная проволока подается под углом @ 
к оси плавящегося электрода х. Это позволяет осям присадочной и электродной проволок пересекаться на 
заданном расстоянии Х от их токоподводящих мундштуков и задавать расстояние х. от точки пересечения осей 
проволок до поверхности металла. В пространстве моделирования задаются также скорости подачи и диаметры 
проволок: электродной у, 4, и присадочной у, 4, Для основной дуги а задается напряжение ИО.» а для 
дополнительной Б — Ис. Эти напряжения регулируются ключом К, перераспределяющим напряжение И, 
источника питания между дугами. Непрерывность токов дуг /„, [, обеспечивается дросселями Г., [у и 
замыкающими диодами Д., О, Токи [ь, [; дуг самоустанавливаются эффектом саморегулирования длин [ь, /, 
соответствующих дуг. Длины вылетов электродов Х., Х, определяются процессом плавления электродов. 
Скорости плавления определяются токами дуг [., [- и падением напряжения в анодных областях, а также 
подогревом вылетов проходящими по ним токами дуг. К — электрическое сопротивление источника питания. 

Математическая модель саморегулирования процесса наплавки. Основное управляющее воздействие 
на процесс обеспечивает работа ключа К. Его замыкания (с периодом т и длительностью %, 1,) регулируют 
напряжение на электродах: 

0. = (0. — (Векь); Им = (0. — (В+. В), (1) 

где К., К; — электрические сопротивления вылетов электродной и присадочной проволок. 

В установившемся режиме эти напряжения определяют длины дуг [11]: 


= (И — Иа — к); 14 = = (Мы - Иа — И, (2) 
где И., О; — падение напряжения в анодной и катодной областях дуги; Ё — градиент потенциала в столбе дуги. 
Средние значения токов дуг целесообразно определять по скоростям плавления проволок, которые в 
установившемся режиме равны ры их подачи: 
тета 


= 
с (ер(Ть — Ту) + йе [5 40а (срСТь = Тхе) т т. (3) 
где 7» — температура капель а а металла; ср — теплоемкость металла проволок; 41/5 — теплота 
плавления; Гхе, [ху — температуры подогрева вылетов проволок токами дуг. 
Следует учитывать, что при импульсном регулировании ток будет изменяется с периодом т и 
амплитудой: 








Оет 


А/= 1 (4) 
Однако при малом периоде изменения, характерном для современных источников питания, 
периодическим изменением тока можно пренебречь. 
Ток, протекающий через катодное пятно на поверхности металла: 
ШЕИ. (5) 
Чтобы узнать распределение температуры вдоль осей х вылетов электродной Ге (х) и присадочной Т/(х) 
проволок, следует учесть зависимость удельного электрического сопротивления металла проволоки от 
температуры р. (Г) [12]. Тогда распределение температуры: 


т.о = (#8) 2 [ре (т,60) 4 ТУ - (9. [* рь(1,(%)) 4х. (6) 


Температура подогрева 7. определяется по ее значениям на конце вылета х = Х, 





угер` 0 


Как отмечается в работах [13, 14], температура капель Т;, металла, поступающего в сварочную 
(наплавочную) ванну, определяется сочетанием многих факторов. Однако она в любом случае будет в 
интервале между температурами плавления 17! и кипения Ту металла: 


Ту + ТЕ 
Г = 3 (7) 
Мощность тепловыделения рассчитывается по приведенным ниже формулам. 
— Для катодного пятна дуги: 
Рь = Ик(Е — 1) (8) 
— Переносимое электродными каплями: 
па 
Ра — т сРГ „Те =. (9) 
— От потока капель присадочной проволоки: 
Н. 
РЕ = Ир + срТ хуй т (10) 
Мощность излучения столба дуги: 
Ра = Е(1-1+ + 6). (11) 


На поверхность подложки попадает только часть излучения дуги. Учтем, что интенсивность излучения 
обратно пропорциональна квадрату расстояния от центра излучения. Примем, что этот центр удален от 
поверхности подложки на длину основной дуги. Тогда мощность на поверхности подложки: 
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Руза 45 
12 
2 И 72+ и. ( ) 
где А) —_ Площадь поверхности, г — расстояние от поверхности до оси факела дуги. 


Численное решение уравнений модели сводится к итерационному подбору токов дуг, при котором 
достигается равенство во всех соотношениях (рис. 2). 


Исходные данные: марка металла и наплавляемого материала, 

размеры А\1, Х›, х., параметры режима наплавки у», 4», Ув а Ць, К, т, Ш, 1... 

Определенные из баз данных свойства материалов Т,, Ти, ср, 415, ЦИ, Чь Е, ре(Т). 

Начальное приближение 

Итерационный цикл уточнения токов /ь, // 
Вычисление распределения температур в вылетах: х = 0, Г(х) = То, К. = 0, К,=0 
Цикл х =х + ах 

ре(т;@) 

У рср 


Чх 


то = теб + (№) 2, т.) = пуб + (#+) 


Электросопротивление вылетов: 


врЁ 
№ р Ре (Те(%)) ‚В.о = ва — (10) 


Чх 
4 ла 


В = ЕО 


пока х < Хь илих < Х; 
Параметры вылетов 7%. = Ге (№е), Ги = ТАЁ», Ве = В, (Ёе), Ку = К (Ё) 
Уточнение токов 


1 па? 
р 31; В 4 (т, — Ту’) Е и) = о е (срСТ, И Тхе) + ‘.) 
(#1 


4 т 
пока /. соп${ или /[+с01п$ 
Определение параметров процесса: О = Аь (6 + 1, ОЙ. 
РИ + ерТть Е; р = бу Неро" 
В. = (В+) Дент: |: = (и, 2-1); = =(Иу- И. - 0%). 


Рис. 2. Алгоритм расчета параметров совместного процесса плавления электродных проволок 





Представленный на рис. 2 алгоритм расчета параметров совместного процесса плавления электродных 
проволок позволяет детально исследовать условия саморегулирования основной и дополнительной дуг в 
процессе наплавки. 

Результаты исследования 

Регулировочные характеристики процесса. Регулирующие воздействия: 

—щ скорость подачи электродной и присадочной проволок, 

— напряжение холостого хода источника питания, 

— относительная длительность замыканий ключа, управляющая соотношением напряжений основной 
и дополнительной дуг. 

Эти воздействия определяют токи дуг, их длины и распределение мощности процесса. 
Дополнительными параметрами, определяющими результаты регулирования, являются диаметры проволок, их 
вылеты, углы между ними. 

Диапазоны значений параметров регулирования имеют ограничения. Важнейшее из них — условие 
существования дуги между подложкой и электродом, т. е. тока /№ = /. — [+> 0 через катодное пятно на подложке. 

Соответственно, должно выполняться соотношение: 


а2 уг 
Е А 


| (13) 
Па Ок 
Второе важное ограничение — минимальное напряжение на дугах, обеспечивающее их стабильное 
горение при отсутствии коротких замыканий: 
Ге 9 


Дополнительно нужно обеспечить корректное расположение точки схода электродов над поверхностью 


подложки и наплавляемого слоя и отсутствие коротких замыканий каплями электродного металла, т. е. длина 
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основной дуги должна быть: 
[> хе вх. > 4... (15) 

Для этого необходимо напряжение основной дуги не менее: 

И. = > Иа Ин ЕЁ. (16) 

Длина дополнительной дуги [, является технологическим параметром, позволяющим удалять поток 
капель присадки от основной дуги. 

Блияние регулирующих параметров на мощность процесса. Мощность тепловых потоков, 
воздействующих на подложку, складывается из тепловыделения в катодном пятне основной дуги, мощности 
излучения основной и дополнительной дуг, теплопереноса каплями электродного и присадочного металла. На 
рис. 3 показано влияние скорости подачи электродной проволоки диаметром 1,6 мм из сплава шсопе] 625 на: 

— ток /, основной дуги, 

— мощность Ру, выделяемую в катодном пятне, 

— мощность Р, теплового потока капель электродного металла, 

— температуру Г. подогрева вылета электрода, 


— длину [/, основной дуги. 
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Рис. 3. Зависимость параметров основной дуги от скорости подачи у, электрода диаметром &, = 1,6 мм из сплава шсопе 625 
при вылете /[, = 10 мм, напряжении питания (И, = 32 В, скважности {,/т = 0,5 и токе дополнительной дуги /.= 108 А 


С увеличением скорости у, подачи электродной проволоки при фиксированном токе /. дополнительной 
дуги ток /, протекающий по электроду, возрастает почти пропорционально скорости. Отклонение от 
пропорциональности объясняется заметным увеличением температуры Т. вылета электрода. По этой же 
причине мощность Р‚„, нагревающая электрод, увеличивается быстрее, чем токГ, дуги. Мощность Рь, 
выделяемая в катодном пятне на поверхности подложки, практически пропорциональна скорости подачи, но 
много меньше тепловой мощности Р, потока капель. Длина /, основной дуги линейно убывает по мере 
увеличения скорости подачи электрода. 

На рис. 4 показано влияние скорости у; подачи присадочной проволоки диаметром 1,2 мм из сплава 
Гпсопе| 625 на: 

— ток /+ дополнительной дуги, 

— мощность Ру, формируемую в катодном пятне, 

— мощность Р‚теплового потока капель металла присадочной проволоки, 

— температуру /, подогрева вылета проволоки, 


— длину ‚дополнительной дуги. 
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Рис. 4. Зависимость параметров дополнительной дуги от скорости подачи у: присадочной проволоки диаметром 4, = 1,2 мм из 
сплава [псопе! 625 при вылете /, = 10 мм, напряжении питания Ц. = 32 В, скважности („/ т= 0,5 и токе основной дуги /. =185 А 


С увеличением скорости у; подачи присадочной проволоки при фиксированном токе /, основной дуги 
ток /- несколько непропорционален скорости подачи из-за заметного увеличения температуры 1; вылета 
проволоки. Мощность Р., затрачиваемая на нагревание проволоки, увеличивается быстрее тока /; 
дополнительной дуги. Мощность Ри формируемая в катодном пятне на поверхности подложки, сильно 
снижается при увеличении скорости подачи присадочной проволоки из-за перераспределения тока основной 
дуги на присадочную проволоку. Тепловая мощность Р; потока капель с присадочной проволоки меньше 
мощности Р. капель с электрода. Длина /- дополнительной дуги нелинейно убывает по мере увеличения 
скорости подачи. Это объясняется влиянием электрического сопротивления, которое у присадочной проволоки 
больше, чем у вылета электрода, из-за меньшего диаметра проволоки и большей температуры 7‚. Напряжение 
дуг определяет их длину, т. е. расположение над подложкой. 

На рис. 5 показано, как от режима наплавки зависят расстояния между электродной и присадочной 
проволоками и их удаление от подложки. 
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Рис. 5. Влияние напряжения основной дуги Оз, при разном напряжении дополнительной дуги Она: мощность излучения дуги 
Ра ‚ мощность излучения дуги на подложку Р‚„,, длину основной дуги /», расстояния х, от точки схода электродов до подложки 


Полосков С. С. и др. Регулирование дугового взаимодействия между плавящимся электродом и дополнительной присадочной проволокой 





Минимальное напряжение дуг (12 В) определяется суммой анодного и катодного напряжений, а также 
их минимальной длиной. Между электродной и присадочной проволоками задано напряжение Це. При этом 
повышение напряжения О, между электродом и подложкой увеличивает длину [. основной дуги, что 
сопровождается повышением мощности Р„„ тепловыделения в плазменном факеле дуг. Повышение 
напряжения ЦИ, между электродом и присадкой увеличивает длину /, дополнительной дуги, что также повышает 
мощность Ра. Некоторая часть Р„, 
зависимость общей мощности Р„„ факела дуг от напряжения основной дуги О.,. При этом менее значительно 


этой мощности Р,„.4 излучается на поверхность подложки. Заметна 


изменение мощности Р,,, нагревающей подложку. Это объясняется удалением плазмы дуги от подложки при 
увеличении длины /, основной дуги. Напряжение дополнительной дуги Це, влияет на мощность подогрева Р,., 
излучением значительно сильнее, так как расстояние между факелом дуги и подложкой не зависит от этого 
напряжения. Напряжения Ц, и И’ существенно влияют на расположение х. точки схода проволок, что налагает 
ограничения по питанию дуг. 

Обсуждение и заключения. Наплавка с дуговым взаимодействием между электродной и присадочной 
проволокой позволяет регулировать мощности, расходуемые на плавление наплавляемого металла и на 
подогрев подложки. Стабильность дуговой сварки плавящимся электродом обычно обеспечивается эффектом 
саморегулирования, который выравнивает скорости плавления и подачи электродных проволок при их 
изменении в широких пределах. Физико-математическая модель совместного плавления электродной и 
присадочной проволок разработана с учетом саморегулирования, при котором управляющими воздействиями 
являются скорости подачи проволок и напряжения на сварочных горелках, а токи и длины дуг — это 
результаты саморегулирования. Модель позволяет определять значения регулирующих параметров, а также 
мощности теплового воздействия на подложку при наплавке плавящимся электродом с дополнительной 
присадочной проволокой. 

По итогам проведенных работ можно сделать следующие выводы. 

1. При дуговой наплавке плавящимся электродом с дополнительной присадочной проволокой для 
эффективного регулирования соотношения между объемами наплавляемого материала и переплавляемой 
подложки следует обеспечить саморегулирование: 

— основной дуги между электродом и подложкой, 

— дополнительной дуги между электродами. 

2. Саморегулирование дугового процесса плавления электродной и присадочной проволок описано 
системой уравнений. В ней регулирующие воздействия — это скорости подачи проволок и напряжения 
основной и дополнительной дуг, а результатом управления являются токи дуг, их длины, распределение 
мощности тепловыделения, тепловые потоки капель электродного металла и излучения дуг на поверхности 
ПОДЛОЖКИ. 

3. Установлены ограничения на: 

— соотношения параметров подачи электродной и присадочной проволоки, которые обеспечивают 
наличие дуги между электродом и подложкой; 

— минимальное напряжение на дугах, поддерживающее их стабильное горение при отсутствии 
коротких замыканий. 

4. Выполнены расчеты режимов наплавки плавящимся электродом с дополнительной присадочной 
проволокой сплава Шшсопе| 625. Вычисления показали возможность регулирования процесса в широком 
диапазоне скоростей подачи проволок и напряжения на горелках. 
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Прочностной анализ корпуса плунжерного насоса типа ТУ 600 
в системе автоматизированного проектирования эЗоПа У\/огк$ Знишайоп 


Э. В. Коледа', С. О. Киреев', М. В. Корчагина', А. В. Ефимов', Й. Шперлинг” 


'ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 
* АСТ \!е1 Зегу1сез ОшЬН (г. Дайсдорф, Федеративная Республика Германия) 





Введение. Актуальность представленной работы обусловлена широким применением плунжерных насосов в 
производственной практике, в частности в газо- и нефтедобыче. От их надежности во многом зависит качество 
проведения технологических операций и эффективность дальнейшей эксплуатации скважин. 
Совершенствование плунжерных насосов предполагает повышение их надежности, увеличение срока службы, 
коэффициента полезного действия, уменьшение габаритов, массы и трудоемкости монтажных и ремонтных 
работ. Модернизация механизма включает его силовое исследование, так как по найденным силам проводятся 
последующие расчеты на прочность. До появления программ численного анализа твердотельных объектов 
аналитическое решение задачи прочностного расчета станины приводной части насоса высокого давления было 
весьма трудоемкой и затратной процедурой. Ситуация изменилась с развитием компьютерных технологий и 
включением метода конечных элементов в системы автоматизированного проектирования. Цель настоящей 
работы — выполнение прочностного расчета корпуса плунжерного насоса типа ТУЗ 600 из стали 09Г2С. 
Материалы и методы. Разработана методика определения реакций опор кривошипного вала плунжерного 
насоса высокого давления и прочностного расчета корпуса приводной части. Направление и величина 
результирующих сил и реакций опор определены графически по принципу суперпозиции силового воздействия 
на опоры. Прочностные расчеты проводились с помощью метода конечных элементов в системе 
автоматизированного проектирования Зо14 Уогк$ Зпищайоп. При этом использовались твердотельная и 
конечноэлементная модели корпуса с наложенными граничными условиями, выявленными при анализе 
конструкции и расчете усилий, возникающих в процессе работы насоса. 

Результаты исследования. Описаны реакции в опорах коленчатого вала с учетом формируемых плунжером сил 
в зависимости от режима его работы и положения кривошипа. Определены силы, действующие на каждый из 
плунжеров, и результирующие реакции в каждой из опор. Представлены эпюры напряжений и коэффициента 
запаса прочности, позволяющие оценить прочность корпуса и выработать рекомендации по созданию более 
рациональной его конструкции. 

Обсуждение и заключения. В результате проведенных расчетов выявлены области конструкции, имеющие 
минимальные коэффициенты запаса прочности, и области, в несколько раз превышающие рекомендованные 
значения. Это позволяет оптимизировать исследуемую конструкцию, усилив первые и уменьшив толщину 
металла на вторых. С точки зрения массогабаритных характеристик и ремонтопригодности результаты 
выполненного прочностного расчета могут быть использованы для оптимизации конструкции корпуса насоса 
при типовых режимах работы. 


Ключевые слова: плунжерный насос, реакции опор, прочностной расчет, оптимизация конструкции, принцип 
суперпозиции, расчет корпусных деталей. 
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Е. У. Кое4а', 5. О. Кгееу', М. У. Когсваоша', А. У. Ейтоу', У. Зрег што” 


Роп Зе Тесбиса! ОштуегзИу (Воз®юу-оп-Ооп, Визап Еедеганоп) 
"АСТ \еП Зегу1сез ОшЬН (ГуездотЕ, Еедега! ВерчБйс оЁ Септапу) 
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Фелг зегулсе Ше, еЁРйслепсу, до\п171е, тедисйоп ш \ееВЬ [абог ищепзИу ог шуаПайоп ап4 гераш \уотКк. ТВе 
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Введение. Добыча нефти не обходится без использования плунжерных насосов. Они необходимы для 
таких технологических операций, как цементирование и кислотная обработка скважин, гидропескоструйная 
перфорация, гидравлический разрыв пласта и др. [1]. Уменьшение габаритов и веса насосов делает их крайне 
привлекательными для использования в мобильных нефтепромысловых установках. От надежности насосов во 
многом зависит качество проведения технологических операций и эффективность дальнейшей эксплуатации 
нефтяных и газовых скважин [2]. 

Несмотря на общий довольно высокий уровень конструкций плунжерных насосов, их продолжают 
улучшать. Совершенствование плунжерных насосов направлено на повышение надежности, увеличение срока 
службы, коэффициента полезного действия, уменьшение габаритов, массы и трудоемкости монтажных и 
ремонтных работ. С этой целью меняются конструкции узлов и деталей приводной и гидравлической части [3— 
5]. Модернизация механизма всегда включает его силовое исследование, так как по найденным силам 
проводятся последующие расчеты на прочность. В последнее время в связи с развитием компьютерных 
технологий все шире используются численные методы прочностного анализа с применением прикладных 


программ [6, 7]. 
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До появления программ численного анализа твердотельных объектов аналитическое решение задачи 
прочностного расчета станины приводной части насоса высокого давления было весьма трудоемкой и 
затратной процедурой [8-10]. Развитие метода конечных элементов (МКЭ) в области механики 
деформируемого твердого тела и включение его в системы автоматизированного проектирования (САПР, 
например Зойа У\Уо!Кз ЗиищаНоп) открывает новые возможности в решении задач подобного рода. 
Прочностные расчеты, проводимые в САПР, позволяют оптимизировать конструкцию корпусных деталей. Для 
рассматриваемой конструкции корпус является одной из ответственных деталей, воспринимающих нагрузку и 
обеспечивающих правильное взаимное расположение элементов приводной части. 

Для получения достоверных результатов при прочностном анализе корпуса насоса с помощью 
численных методов в САПР необходимо определить все внешние нагрузки, действующие на корпус. 

В качестве объекта модернизации выбран широко применяемый трехплунжерный насос Т\УЪЗ 600. В 
России для изготовления корпусных деталей плунжерных насосов высокого давления применяется 
конструкционная низколегированная сталь марки 09Г2С [11]. Замена материала корпусных деталей 
плунжерного насоса типа ТУЗ 600 на сталь марки 09Г2С позволит снизить себестоимость. В связи с этим цель 
настоящей работы — анализ возможности замены материала корпуса плунжерного насоса типа ТУЗ 600, на 
сталь 09Г2С для оптимизации цены и обеспечения возможности ремонтных работ. 

Материалы и методы. Замену материала следует обосновать прочностным расчетом. Расчет корпуса 
проводился МКЭ в САПР ЗоПа У\Уогк$ зиичайоп. 

Внешние силы, действующие на корпус насоса, — это реакции в опорах коленчатого вала. Они 
возникают от действия сил инерции возвратно-поступательно движущихся деталей крейцкопфной шатунной 
группы и сил давления жидкости на плунжер. Для определения реакций в опорах коленчатого вала необходимо 
выполнить динамический расчет кривошипно-ползунного механизма плунжерного насоса". 

Как отмечается в [12, 13], в задачи динамического анализа входит изучение влияния внешних и 
внутренних сил на звенья и кинематические пары механизма, а также выявление способов снижения 
динамических нагрузок. 

Рассматриваемый коленчатый вал состоит из двух коренных шеек на основной оси вала и трех 
кривошипных шеек, эксцентрично расположенных со смещением 120°. Подшипниковые опоры расположены 


на коренных шейках и на щеках вала между кривошипными шейками (рис. 1). 





Рис. 1. Коленчатый вал плунжерного насоса типа Т\/$ 600: 1 — коренные шейки; 2 — шатунные шейки; 3 — щеки 


Во время работы насоса три плунжера последовательно совершают возвратно-поступательные 
движения, обеспечивающие нагнетание или всасывание (рис. 2). 


' Тимофеев, Г. А. Теория машин и механизмов: курс лекций. М., 2010. 351 с. 
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Рис. 2. Расчетная схема механизма приводной части плунжерного насоса: положение 1-го плунжера (а); 
положение 2-го плунжера (6); положение 3-го плунжера (в) 


Рассмотрим такое положение механизма, при котором первый плунжер находится в крайней точке 
нагнетания (кр.т) и испытывает максимальную нагрузку. Угол поворота кривошипной шейки для этого 
плунжера равен нулю (рис. 2 а). Тогда второй плунжер совершает всасывание (всас.) (рис. 2 6), его шейка 
смещена на 120°. Третий плунжер работает на нагнетание (наг.) (рис. 2 в), шейка смещена на 240°. Таким 
образом, в одной ситуации коленчатый вал нагружен силами различной величины, и реакции опор будут 
различны. Это происходит при положении кривошипа, когда первый плунжер находится в крайней точке 
нагнетания, второй совершает всасывание, а третий работает на нагнетание. Проведем расчет для указанного 
положения кривошипа (механизма) [8—10], а затем для других положений: когда второй и третий плунжеры 
будут достигать крайней точки нагнетания. 

Массы звеньев рассчитаны в программе Аию4дезК шуетщог Рго{ез1опа| 2018 по созданной ранее ЗО 
модели привода плунжерного насоса типа Т\УЪ 600. Была построена расчетная схема механизма приводной 
части плунжерного насоса с указанием сил, действующих на его звенья (рис. 2). 

Обозначим внешние силы для нахождения реакций кривошипного вала: Р„, — сила давления жидкости, 
Р, — сила инерции. 

Силу давления жидкости, действующую на плунжер, находим из выражения: 


— 2 
Рк = фк.п Г", 
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е | 
где р. — давление перекачиваемой жидкости, приложенное к плунжеру; 7’ — радиус сечения плунжера. 
Сила инерции возвратно-поступательно движущихся частей равна: 


Р, = -т-1= —т-г: о? - (с0о5ф + (2: с05*ф — 1), 
т = ть + 0,275 ть, 
где т, — масса крейцкопфа, т. — масса верхней головки шатуна, г — радиус кривошипа, « — частота 
К Ш 


вращения коленчатого вала при 230 об/мин, ф — угол поворота коленчатого вала. 

Сила, действующая на плунжер, определяется по формуле: 

Р; = Ра + Рш, 

где Р„; — сила давления перекачиваемой жидкости в каждом положении плунжера, Р‚; — сила инерции 
возвратно-поступательно движущихся частей в каждом положении плунжера. 

Суммарная сила Р;, приложенная к оси пальца крейцкопфа и направленная по оси цилиндра, может 
быть разложена на: 

— силу №, действующую перпендикулярно оси цилиндра; 

— силу 5, действующую по оси шатуна. 

№ прижимает крейцкопф к стенке цилиндра, что вызывает износ их поверхностей. Эта сила изменяется 
по значению и направлению, поочередно прижимая крейцкопф то к одной, то к другой стороне направляющих. 

Силу о, перенесенную на ось шатунной шейки, можно разложить на: 

— касательную силу Т, действующую перпендикулярно кривошипу коленчатого вала; 

— радиальную силу К, направленную по оси кривошипа. 

Сила К определяется по формуле: 


с05(ф; + В) 
Ко Ро 
со5В 
где ф; — угол поворота кривошипа (рис. 2); В; — угол отклонения шатуна от оси (рис. 2); Ру; — суммарная 


сила, действующая на плунжер. 
Результаты расчета радиальных сил для всех кривошипных шеек в каждом рассматриваемом 
положении механизма сведены в таблице 1. 
Таблица 1 
Результаты расчета радиальных сил для всех кривошипных шеек в каждом рассматриваемом положении 
механизма 


ПО ПО ОЕ ООО ООО ООО О ОО ОЕ О ООО 





Реакции опор направлены противоположно радиальным силам (рис. 3), а их величины определяются по 
формулам: 


— К1 

КВ = № 
К 

! м.) 

К>, Кз — 2? 
К 

! м З 

Кз, Кл О? 


2 
где К\, К>, Кз — радиальные силы, направленные по оси кривошипа. 


'Каталог плунжерных насосов высокого давления производства У\Уен ЗРМ/ ООО «Норд-СПМ»// виерспм.рф: [сайт]. ОВГ: 
Вр://виерспм.рф/саа]озиез/ Л21К_КАТАЛОГ НАСОСОВ.рАЕ (дата обращения: 14.02.2021). 
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Рис. 3. Схема реакций опор 


Направление и величина результирующих сил РЁ. и Е. реакций опор 02, Оз определяются графически 
(рис. 3) по принципу суперпозиции силового воздействия на опоры. 

Результирующие силы ЕЁ. и Ё; по величине и направлению совпадают с силами КЁ. и К- соотвественно. 

Найдены численные значения и направления опорных реакций для положения механизма, при котором 
первый плунжер находился в крайней точке нагнетания (рис. 2 а), второй совершал всасывание (рис. 26), а 
третий — нагнетание (рис. 2 в). Для других положений механизма реакции опор определяются аналогично и 
будут численно равны, поочередно меняясь местами. 

В таблице2 представлены результаты расчета реакций опор коленчатого вала для положений 
механизма, в которых каждый из плунжеров поочередно находился в крайней точке нагнетания. 

Таблица 2 
Численные значения реакций опор для всех положений механизма 


Силы реакций опор 


Положение И 
оложение плунжера 


225000 225000 153343 153343 





153343 329602 225000 1251 
1251 153343 329602 225000 


Корпус плунжерного насоса высокого давления моделировался в САПР Ашо4езК шуешог 


Рго!е5$10опа1 2018. Модель включает стаканы подшипников, т. к. наличие стакана существенно влияло на 
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характер приложения нагрузки. 
Прочностной анализ выполнен МКЭ в системе автоматизированного проектирования ЗоПа Уо1К$ 


опич[аноп. 
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На рис.4а представлена твердотельная модель рассчитываемой конструкции, на рис.46 — 
конечноэлементная. 


Стаканы подшипников 


Корпус насоса 





6) 


Рис. 4. Модели корпуса трехплунжерного насоса высокого давления Т\/5 600: твердотельная (а), конечноэлементная (6) 


Конечноэлементная модель содержит 100747 элементов (175172 узлов). Материал корпуса — сталь 
09Г2С. Модель описывает нагрузки и фиксирует перемещения по всем направлениям поверхностей отверстий 
под болты. На рис. 5 представлены граничные условия для контактирующих поверхностей. По типу контактов 
поверхности связаны с общими узлами. 
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Рис. 5. Граничные условия: 1 — рабочая нагрузка (радиальная сила, распределенная по цилиндрическим поверхностям); 
2 — дистанционная нагрузка (масса редуктора и гидроблока); 3 — ограничение — «зафиксирован», нулевые перемещения 
по всем осям — Х, У, 1; 4 — рабочая нагрузка (сила давления жидкости во время нагнетания и всасывания); 5 — сила 
тяжести; 6 — тип контакта — связанные (зазор отсутствует) поверхности с общими узлами 


Результаты исследования. Эпюры эквивалентных напряжений по Мизесу представлены на рис. 6. 
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Рис. 6. Эпюры, полученные в результате прочностного анализа корпуса плунжерного насоса типа Т\М$ 600: эпюра 
эквивалентных напряжений (а), эпюра коэффициента запаса по пределу текучести (6) 


Коэффициент запаса по пределу текучести находится из соотношения: 





(] 
К =—>, 
Оэкв 
где о, — предел прочности стали марки 09Г2С, с, = 470 Мпа; о, — максимальное эквивалентное 
напряжение, о; = 396,392 МПа. 
К = = 118. 
” 396,392 


Как видно, эквивалентные напряжения’ имеют наибольшую интенсивность в нижней области 
конструкции корпуса насоса. 
Минимальный коэффициент запаса по пределу текучести материала — 0,82 (рис. 6 6). 
Прочностной расчет корпуса плунжерного насоса проводился для третьего положения механизма 
(рис. 2 в). Для других положений механизма (рис. 2 а, 6) прочностной расчет выполняется аналогично. 
Результаты прочностного расчета корпуса для трех положений механизма представлены в таблице 3. 


Таблица 3 
Результаты прочностного расчета корпуса плунжерного насоса 


Положение Эквивалентные Коэффициент запаса по пределу 
Положение плунжера 


механизма напряжения. МПа 
Р : текучести прочности 


231,772 


284,783 
396,392 


Обсуждение и заключения. Описаны реакции в опорах коленчатого вала с учетом сил, формируемых 





плунжером в зависимости от режима его работы и положения кривошипа. Определены силы, действующие на 
каждый из плунжеров, и результирующие реакции в каждой из опор. 

Для прочностного расчета корпуса плунжерного насоса из стали 09Г2С применен МКЭ в САПР ЗоПа 
УМ’отК$ Зиач[аноп. При оценке прочности корпуса отмечено, что в третьем положении механизма нижняя 
область конструкции имеет минимальные коэффициенты запаса по пределу текучести и прочности (таблица 3). 


' Марочник сталей и сплавов / А. С. Зубченко, М. М. Колосков, Ю. В. Каширский [и др.]. М., 2003. 784 с. 
° Фещенко В. Н. Справочник конструктора. Книга 2. Проектирование машин и их деталей. М., 2017. 400 с. 


Машиностроение и машиноведение 


1 
х® 


ЬИр://уезш-допза.га 


А4уапсеа Епатеетиз Кезеагсй 2021. Т. 21, № 1. С. 71-61. 155М№ 2687—1653 





Указанная область требует усиления, поскольку значения коэффициентов недостаточны для заданных условий 
эксплуатации. Также в корпусе есть ненагруженные участки, где значения коэффициентов по запасу прочности 
и текучести в разы превышают рекомендуемые. В дальнейшем следует оптимизировать корпус плунжерного 
насоса: усилить области конструкции с неудовлетворяющими коэффициентами по запасу прочности и 
текучести, а в местах, где коэффициенты превышают рекомендуемые значения, применить металл меньшей 
толщины. 

Таким образом, результаты выполненного прочностного расчета могут быть использованы для 
оптимизации конструкции корпуса насоса при типовых режимах работы. 
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Силовой расчет цепной передачи с зубчатой цепью и эвольвентными 
звездочками 





С. Б. Бережной, Г. В. Курапов 


ФГБОУ ВО «Кубанский государственный технологический университет» (г. Краснодар, Российская Федерация) 


Введение. В силовом анализе цепной передачи с зубчатой цепью и эвольвентными звездочками учтены 
центробежные силы и коэффициенты трения между пластиной и зубом звездочки. Цели работы — определение 
всех сил, действующих в зацеплении, и расчет коэффициента сцепления зубчатой цепи с эвольвентной 
звездочкой в приводных передачах. 

Материалы и методы. При оценке тяговой способности цепной передачи важны такие параметры, как силовой 
расчет и коэффициент сцепления зубчатой цепи со звездочкой (показывает, во сколько раз давление на данный 
зуб больше давления, воспринимаемого впереди идущим зубом). В представленной работе визуализированы 
схемы: расположения пластин в зацеплении с зубьями эвольвентной звездочки и сил, действующих в таком 
зацеплении. Рассмотрены факторы, которые действуют на эвольвентный профиль зуба звездочки. Это вес 
пакета пластин цепи и силы: центробежная, трения, нормального давления и растяжения. При этом учтены 
изменения коэффициента сцепления на последующие зубья, участвующие в тяге. Исследовано равновесие 
звеньев Ги 1+ | зубчатой цепи в системе координат ХОТ с центром на оси вращения эвольвентной звездочки. 
Метод позволяет определить все искомые силы геометрическим расчетом значений угловых параметров 
передачи. С помощью полученных уравнений уточняются: коэффициент сцепления В, натяжение ведущей 
ветви передачи 51 и ведомой ь.. 

Результаты исследования. Представлен патентованный стенд передачи с зубчатой цепью и эвольвентными 
звездочками. Проведенные на нем испытания подтвердили итоги исследований цепной передачи с зубчатой 
цепью и эвольвентными звездочками с заданными параметрами. Была доказана правильность силового анализа 
передачи с учетом центробежных сил и коэффициентов трения пластин и зубьев звездочки. 

Обсуждение и заключения. Отмечено, что центробежные силы и коэффициенты трения при зацеплении 
существенно влияют на тяговую способность передачи с зубчатой цепью и эвольвентными звездочками. 
Полученные данные могут быть использованы для точной оценки тяговой способности таких передач. 


Ключевые слова: цепная передача, зубчатая цепь, эвольвентная звездочка, пластина цепи, тяговая способность, 
коэффициенты сцепления, центробежные силы, силы трения, коэффициенты трения, угловая скорость шарнира, 
ведущая ветвь, ведомая ветвь. 
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асНоп сое слет оЁ Ше р!а{е$ ап Ще зргоске {ее \аз ргоуе4. 

Обси5яоп апа Сопсизоп5. [1$ пе фа Фе сепёРига1 Гогсез апа Фе Ёлсйоп сое слет 4итих епоасетеп{ аНес 
1отиНсап у Фе тасйоп сарасйу оРа гап$1п115510п УИ а 1юое4 сваш ап4 шуо[\ие зргоске$5. Тве жа оМаше4 сап 
Бе изеа фо ассигже]у езитае Пе тасйоп сарасиу оЁ зисВ сеагв. 


Кеунога$: спал ап$1115$101, юофе4 сваш, шуопие зргосКеь сваш ре, Насйоп сарасйу, соирПие сое сев, 
сепаРиса| Гогсез$, Нлсйоп Ююгсез, ИлсНоп соеЁйслен, ] о апошаг у@юоску, дтушо 14е, $асК з14е. 


Еог сйапоп: 5. В. Веге7Ипоу, СО. У. Кигароу. Ро\ег апа[уз$1$ ог сВаш гап$1п71$$10п У сеаг сваш ап шуопие 
5ргосКе5. Адуапсеа Епошеегие ВезеагсВ, 2021, уо1. 21, по. 1, р. 82—88. №рз://401.0г°/10.23947/2687-1653-2021- 
21-1-82-88 


Введение. Цепные передачи с зубчатой цепью недостаточно исследованы. Их выпуск ограничен, 
однако они все шире используются в производстве. Например, такими передачами комплектуют приводные 
механизмы в машинах, в прецизионных станках и другой технике. В представленной работе исследуется 
силовой анализ и расчет тяговой способности передач с зубчатой цепью и эвольвентными звездочками. 

Тяговая способность определяет работоспособность эвольвентной зубчатой цепной передачи, которая 
зависит от ряда факторов’. В их числе: 

— геометрия передачи, 

—щ предельное контактное напряжение наиболее слабого элемента пары «пластина цепи — зуб 
звездочки», 

— коэффициент сцепления зубчатой цепи с зубьями звездочки, 

— наличие пробуксовки зубчатой цепи по зубьям звездочки, 

— допустимое удельное давление внутри шарнира, 

— ударная стойкость элементов зубчатой цепи” [1-3]. 

Материалы и методы. Коэффициент сцепления — это один из основных факторов тяговой 
способности цепной передачи, оснащенной зубчатой цепью. Во время работы такой передачи возникает 
зацепление и действуют силы, которые перемещают пластины цепи по рабочему и затылочному профилям 
зубьев звездочки [4—7]. На дуге обхвата эвольвентной звездочки зубчатой цепью любые ее звенья (звено 1, 1+ 1, 
+2 ит. д.) испытывают: 

— силу трения №; цепи о зуб звездочки, 

— растягивающую силу О; соседних звеньев, 

— силу нормального давления №; профиля зуба звездочки (рис. 1). 


' Курапов Г. В. Исследование работоспособности цепной передачи с зубчатой цепью и эвольвентными звездочками : дис. канд. тех. наук. 
Краснодар, 2019. 173 с. 

* Семенов В. С. Обеспечение работоспособности роликовых цепных передач на стадии их эксплуатации // Сб. мат-лов всерос. науч.-практ. 
конф. аспирантов, докторантов и молодых ученых. Майкоп, 2017. С. 121-125. 
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Рис. 1. Схема расположения пластин в зацеплении с зубьями эвольвентной звездочки 


В литературе коэффициент сцепления рассчитывается отношением растягивающих сил О; и О в 


смежных звеньях цепи: 


[6 
В = —. 1 
| О | 


Данный коэффициент непостоянен, поэтому пластины зубчатой цепи перемещаются по профилям 
зубьев звездочек на всей дуге обхвата. На рис. 2 представлено зацепление первого пакета пластин (звена) цепи 


(№ 
КК 


с зубом звездочки. 
у 








т | 















Рис. 2. Схема сил, действующих в зацеплении пластины зубчатой цепи со звездочкой 


На эвольвентный профиль зуба звездочки действуют перечисленные ниже факторы [8—14]. 

1. Сила трения №; /, препятствующая перемещению пакета пластин к вершине зуба звездочки. Здесь /— 
коэффициент трения. 

2. Вес пакета пластин цепи (4, = т;с. Здесь © — ускорение свободного падения, т; — масса звена цепи 
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в зависимости от шага'. 
3. Сила нормального давления М№; профиля зуба звездочки, сосредоточенная в месте зацепления К}. 
4. Растягивающие силы О; и О; смежных звеньев зубчатой цепи. 


т; о? 





5. Центробежная сила, учитывающая массу звена зубчатой цепи Ё, = . Здесь & — угловая скорость 


вращения звездочки. А — радиус центральной части звена (представлена как окружность, меняющая положение 
при вращении звездочки) на эвольвентном профиле зуба звездочки. 
Однако при переходе цепи с ведущей на ведомую ветвь растягивающие силы 0; и О» изменятся на 
величину каждого значения коэффициента сцепления В,, поэтому его нужно определить точнее. 
Построим систему координат ХОТ. Она проходит через центр оси вращения эвольвентной звездочки. 
Вектор оси Х совпадает с вектором линейной скорости зубчатой цепи в передаче. 
Проанализируем равновесие первого звена Г и второго 1 + 1. 
Равновесие сил на оси Х (УХ = 0) будет равно: 
А АВ РОх =0. (2) 
Равновесие сил на оси У (У. = 0) будет равно: 
СС, М ОЕ 


31 аа — НУ =0. (3) 
При УХ =0и У 7 = 0 пластины зубчатой цепи будут неподвижны. Если же УХ =0и У У2 0, то пластины 
будут двигаться по профилю зуба. 
Определим силы, действующие на звенья цепи. На рис. 2 обозначены: Р — центр масс звена зубчатой 


а, 


с050 





цепи; 4. — диаметр окружности места контакта К пакета пластин с профилем зуба звездочки, 4. =  @— 


ИИ 
180 ` 


ят—— 
ре 


диаметр делительной окружности звездочки, 4, = 


Ниже перечислены силы, действующие вдоль осей координат. 


Вдоль оси Х: Вдоль оси 7: 
№, =№, -зшВ, М№,, = М, -созВ, 
№ „Л =№, /-созВ, М№,Л = М, Л -зшВ, 
а Йа ‘зшф, ее "СО5Ф, 
Е А -зШф, а . "с05Ф, 
Ох =О.. "С05Ф, От =О., ‘5шФф. 
Здесь В — угол профиля зуба эвольвентной звездочки; ф = ^ — половина углового шага зубьев эвольвентной 
2 
звездочки. 
Значения угловых параметров: 
Л=Ф+Уу, (4) 
и=д-—у-—0,, (5) 
У=+1+0,, (6) 
2 
п 
р 9, › (7) 
й 
И 
=“ Ф-Уу. в 
В в. (8) 
Здесь у — половина угла заострения зуба эвольвентной звездочки. 
у = 0 — 9, ь (9) 
© = агссо$—^. 10 
ь 7. (10) 


к 


0, — угол между осью ОО, зуба и линией, соединяющей точку К\ с центром О звездочки (рис. 2). 


"ГОСТ 13552-81. Государственный стандарт союза ССР. Цепи приводные зубчатые. Технические условия / Государственный комитет 
СССР по стандартам. — М.: Изд-во стандартов, 1987. 12 с. 
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ый (11) 


5; — толщина зуба эвольвентной звездочки на месте контакта К с зубчатой цепью (рис. 1). аь = т2сС0$49 — 
диаметр основной окружности эвольвентной звездочки, где т — модуль звездочки для передач с зубчатой 
цепью, 2 — количество зубьев звездочки. 

Подставив все значения, входящие в уравнения (2) и (3), получим: 


24, -Е, -с08ф- №, 7 -зтВ- №, -созВ+О., -зтФ-Ё -с08ф=О0, (12) 
(В; -со5ф- №, -созВ- №, Л -зтВ+А „, -с0$ф-20С, ) 
О РТ С иоажи- (13) 
5шф 
О =М, -созВ- №, Л-со$зВ- Е, -заф+ЁЕ и -зШФ+О.. -с05ф=0. (14) 


Преобразуя уравнения (13) и (14), получим коэффициент сцепления В; зубчатой цепи с заданным зубом 


эвольвентной звездочки: 


О ‚ М, -с08р-М Х-с05р-Е, -зто+ЁЕ, -ЗПО+О.. -с05ф 
В = = ми М). (15) 
О. Е, -с0$ф-+ М, -с05 В+ М, 1 - за Р-Р, -с08ф-20, 


Если известна величина натяжения ведущей ветви цепной передачи 51, то натяжение ведомой ветви 5> 
определяется с учетом общего коэффициента сцепления цепи с зубом звездочки: 


5 
в 16 
75. (16) 

5 
$ =, 17 
В (17) 
В. =В,-В,-В....В,. (18) 


Здесь В. — общий коэффициент сцепления цепи со звездочкой; и — количество зубьев ведущей эвольвентной 
звездочки. 

Результаты исследования. Для проверки расчетов был спроектирован и изготовлен стенд’, 
защищенный патентом Российской Федерации (рис. 3). Проведены стендовые испытания исследуемой цепной 
передачи. Параметры: 

1) мощность электродвигателя М№,, = 12 кВт, 

2) шаг зубчатой цепи [,= 19,05 мм, 

3) количество зубьев ведущей и ведомой звездочки 21 = 2>= 23, 

4) частота вращения ведущей и ведомой звездочки и! = и = 1000 об/мин. 





Рис. 3. Испытательный стенд передачи с зубчатой цепью и эвольвентными звездочками: | — ведущая звездочка, 
2 — ведомая звездочка, 3 — зубчатая цепь, 4 — электродвигатель, 5 — генератор 


' Стенд для испытания цепных передач : патент 147242 Рос. Федерация : ©01М 13/02 / С. Б. Бережной [и др.]. № 2014134006/28 ; заявл. 
19.08.14 ; опубл. 27.10.14, Бюл. № 30. 3 с. 
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23.4 
‚^^. Результаты 


Для стендовых испытаний проводились вычисления с учетом справочных данных 
представлены в таблице 1. 
Таблица 1 


Полученные данные для определения тяговой способности цепной передачи с зубчатой цепью 


Параметр Обозначение и размерность Численное значение 
Натяжение рабочей ветви передачи 4662 


Общий коэффициент сцепления цепи с ведущей В 172 
эвольвентной звездочкой 7 
Центробежная сила цепи Вы Н 29,82 


Предложенный подход к определению коэффициента сцепления подтвердил верность силового расчета 





с учетом центробежных сил и коэффициентов трения пакета пластин и зуба звездочки. 

Обсуждение и заключения. В данной работе определены все силы, действующие в зацеплении пакета 
пластин цепи и зуба звездочки. Представлен порядок силового расчета передачи с зубчатой цепью. Получены 
зависимости для определения влияния центробежных сил и коэффициентов трения на конечный расчет 
коэффициента сцепления. Для подтверждения теории цепная передача исследована на испытательном стенде. 

Итак, центробежные силы и коэффициенты трения при зацеплении существенно влияют на тяговую 
способность передачи с зубчатой цепью и эвольвентными звездочками. 
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О выявлении трещиноподобных дефектов сварки существующими 
методами контроля качества 
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'ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 
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Введение. Обобщены материалы исследований, посвященных обеспечению работоспособности сварных 
конструкций опасных технических устройств — элементов нефтегазодобывающего оборудования. 
Систематизированы и представлены в виде схемы факторы, определяющие эксплуатационную прочность 
сварных соединений и конструкций. Описаны причины снижения их эксплуатационных свойств. 
Определена роль в обеспечении жизненного цикла конструкции объемных и трещиноподобных дефектов 
сварки, не обнаруженных при контроле качества на стадии изготовления и при определении остаточного 
ресурса в процессе эксплуатации. Предложены меры по улучшению практики выявления 
трещиноподобных дефектов при визуально-измерительном и ультразвуковом методах контроля качества 
сварных соединений. 

Материалы и методы. В соответствии со СНиП 2.05.06-85 конструкции магистральных 
нефтегазопроводов, вертикальных цилиндрических резервуаров и другого нефтегазодобывающего 
оборудования (НГДО) изготавливают из спокойных и полуспокойных низкоуглеродистых и 
низколегированных сталей с временным сопротивлением до 686,5 МПа. Независимо от класса и уровня 
прочности стали должны хорошо свариваться всеми способами, предписанными нормативно-техническими 
документами (НТД). При этом независимо от состояния поставки эквивалент углерода С» не должен 
превышать 0,46 %. Методы исследований — расчетно-экспериментальные. Для вычисления напряженно- 
деформированного состояния сварных соединений применялись методы: 

— конечных элементов (МКЭ), 

— механики разрушений с использованием коэффициента интенсивности напряжений К, . 

При оценке геометрических размеров трещиноподобных дефектов сварки задействованы методы 
математической статистики. Дефекты сварки выявлялись стандартными методами контроля качества, 
предписанными НТД и ГОСТами. 

Результаты исследования. Рассмотрены сварные конструкции, работающие в режиме усталости. 
Проанализированы некоторые характерные для сварных соединений факторы. Показано, как они влияют 
на формирование прочностных эксплуатационных свойств. Установлена роль опасных трещиноподобных 
дефектов, которые с высокой вероятностью могут образовываться в сварных соединениях при 
изготовлении конструкций и при эксплуатации. Часто причиной является снижение показателей 
механических свойств из-за старения и потери пластичности вследствие случайных механических 
воздействий. Отмечено, что невозможно выявлять эти дефекты существующими методами контроля 
качества, как в процессе контроля при производстве, так и при диагностировании в процессе эксплуатации. 
Это снижает точность прогнозирования эксплуатационного ресурса работы сварной конструкции. 
Обсуждение и заключение. Предложено внести в НТД требования: 

— к качеству сварных соединений опасных технических устройств НГДО, 

— по выявлению острых трещиноподобных дефектов с радиусом кривизны в пределах 0,1-0,25 мм. 

Однако достоверность выявления таких дефектов применяемыми методами контроля остается низкой из-за 
человеческого фактора. Для повышения достоверности выявления острых трещиноподобных дефектов 
определены направления совершенствования ультразвукового контроля качества сварных соединений. 


Ключевые слова: сварка, сварные швы, конструкционная прочность, дефекты сварки, неразрушающий 
контроль, визуально-измерительный контроль, ультразвуковой контроль, достоверность. 
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Введение. Проблемам определения прочности эксплуатируемых сварных металлоконструкций 
посвящены научные исследования [1-7], а также стандарты ', инструкции’ и нормативные документы” 
профильных организаций. Факторы, определяющие прочностные характеристики сварных соединений, 
многочисленны и разнообразны. Разрушения сварных металлоконструкций в процессе работы зависят от: 

— условий эксплуатации (уровень расчетных нагрузок, схема нагружения, совместное влияние 
различных сил, действие неучтенных нагрузок и т. п.); 


' ГОСТ 14771-76. Дуговая сварка в защитном газе. Соединения сварные. Основные типы, конструктивные элементы и размеры / 
Федеральное агентство по техническому регулированию и метрологии. М.: Изд-во стандартов, 2001. 37 с. 

"РД-25.160.10-КТН-001-12. Инструкция по технологии сварки при строительстве и ремонте стальных вертикальных резервуаров ОАО «АК 
“Транснефть”» / ОАО «АК “Транснефть”», ООО «НИИ “ТНН”». 112 с. ОВГ: В@рз://\у\у м. збадтед.тга/т9-25-160-10-Кт-001-12-шзгаКстуа-ро- 
{ей по]о21-зуайа-ри-5тое[фуе-1-тетоще-5аштуп-уегикашув-ге7егупагоу-сва${-1_5с68180807а.ВНи (дата обращения :10.10.2020). 

"СТО Газпром 2-2.4-083-2006. Документы нормативные для проектирования, строительства и эксплуатации объектов ОАО «Газпром». 
Инструкция по неразрушающим методам контроля качества сварных соединений при строительстве и ремонте промысловых и 
магистральных газопроводов / ВНИИгаз, ЗАО «ВНИИСТ-Диагностика», ИРЦ «Газпром». ОКГ: 
Брз://е5.згоушЁ.та/таех2/1/4293831/4293831880.В п (дата обращения :10.10.2020). 
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— конструктивных и технологических причин. 
Основные такие факторы представлены в виде схемы на рис. 1. 
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Рис. 1. Факторы, определяющие эксплуатационную прочность сварных соединений и конструкций 


Обобщение многочисленных исследований позволяет утверждать, что на несущую способность 
сварных конструкций влияют перечисленные ниже условия. 

Г. Неоднородность структуры металла и механических свойств сварного шва и околошовной зоны. 

П. Дефекты сварки, главным образом трещиноподобные. Это подрезы, трещины наплывы, острые 
включения, непровары, несплавления и т. п. (рис. 2). 





Рис. 2. Дефекты в стыковом сварном соединении 


Ш. Несовершенство геометрии сечений сварных швов. Снижение прочности в наибольшей степени 
обусловлено остротой и резкостью сопряжения основного и наплавленного металла (рис. 2, 3): Чем меньше 
радиус перехода тг и больше угол подхода 09 от основного к наплавленному металлу, тем ниже усталостная 
прочность сварного соединения. Эта закономерность действует даже при отсутствии дефектов сварки. 












бер ь О шах 


Рис. 3. Распределение напряжений в зоне сопряжения основного и наплавленного металла, в вершине острых дефектов 


Конструкцию сварного соединения определяют: 
— форма и размеры стыкуемых кромок деталей; 
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— взаимное расположение, конфигурация, длина сварных швов; 

— тип сварного соединения по ГОСТ. 

На рис. 4 показано распределение напряжений в нахлесточных сварных соединениях с фланговыми и 
комбинированными швами. 
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а) 6) 
Рис. 4. Распределение напряжений в нахлесточных сварных соединениях: с фланговыми швами (а); 
с комбинированным швом (6) 


В таблице 1 приведены значения теоретического коэффициента концентрации напряжений К; в 
зависимости от длины и ширины накладки: 


бтах р |9) 
КНР = 6.6.1, - < (461). (1) 





Таблица 1 


Значения К; в зависимости от длины и ширины накладки 





ГУ. Остаточные сварные напряжения (ОСН). 

У. Абсолютные размеры свариваемых деталей (масштабный фактор). 

Наличие сварных соединений в металлоконструкции даже без дефектов сварки обычно способствует 
снижению эксплуатационных свойств. Наименьшую опасность в этом смысле представляют стыковые сварные 
соединения. Нахлесточные и тавровые соединения создают концентрацию напряжений, обусловленную 
резкими изменениями рабочих сечений [8, 9]. Однако для изготовления многих конструкций недостаточно 
стыковых швов. 

Остаточные сварочные напряжения могут изменять картину распределения общего силового потока. 
ОСН, взаимодействуя и суммируясь с напряжениями от внешних нагрузок, влияют на работоспособность 
сварных соединений. Это влияние может быть и положительным, например, если векторы ОСН и напряжений 
от внешних нагрузок разнонаправлены (растяжение + сжатие ит. п.). 
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Наличие сварных соединений в наибольшей степени сказывается на усталости конструкций. 
Возмущающие факторы (рис. 1) приводят к разрушению при напряжениях, которые могут передаваться 
элементами металлоконструкций без сварных соединений‘. 

В многочисленных работах, посвященных влиянию дефектов сварки на работоспособность 
соединений, даются рекомендации по снижению вредного воздействия. При этом отмечается, что некоторые 
дефекты следует отнести к разряду неопасных. Например, одиночные некрупные поры практически не влияют 
на усталостную прочность сварных соединений. В этом случае регламенты по наличию пор призваны 
обеспечить не прочность, а приемлемый товарный вид изделия. 

Традиционно под дефектами сварки понимаются любые несовершенства геометрии и несплошности 
материала сварных швов и околошовной зоны. Таким образом, геометрия локальной зоны сопряжения 
основного и наплавленного металла также может считаться дефектом. 

Материалы и методы. Изучение нормативных документов, регламентирующих наличие и размеры 
дефектов сварки, показало, что в них практически отсутствуют требования к геометрии участков сварных 
соединений в месте сопряжения шва с основным металлом. Однако именно здесь зарождаются разрушения из- 
за образования и роста трещиноподобных дефектов [2, 4, 10]. При разработке нормативов (например, для 
магистральных трубопроводов, работающих под действием статических и повторно-статических нагрузок) 
важно учитывать не только размер дефекта, ослабляющего рабочее сечение сварного шва, но и его остроту, 
определяющую локальную концентрацию напряжений в устье дефекта. 

В [8] представлены методики прогнозирования устойчивости сварных соединений при работе на 
усталость, учитывающие реальную геометрию сварных швов и вероятность дефектов. Такой подход 
целесообразно использовать для разработки отраслевых нормативов допустимости дефектов сварки с учетом 
нагруженности и ответственности металлоконструкций. Данные методики позволяют определять статистически 
обоснованные характеристики качества сварных соединений. 

В металлоконструкциях при статических нагрузках дефекты сварных швов проявляются, если их 
размеры столь значительны, что реально ослабляют рабочие сечения соединений. Для определения прочности 
соединений и приемлемых размеров дефектов проводят расчет на статическую прочность по допускаемым 
напряжениям. 

К расчетам сварных соединений, работающих в режиме повторно-статических нагрузок с 
трещиноподобными дефектами, целесообразно подходить с позиций механики разрушения. В этом случае 
используется коэффициент интенсивности напряжений К, [3, 10-12]: 

К, = от, (2) 
где о — номинальные напряжения; / — полудлина эквивалентной трещины, описывающей характерный размер 
трещиноподобного дефекта (например, глубина подреза, глубина непровара ит. п.). 

При этом к трещиноподобным дефектам можно условно отнести дефекты с радиусом в вершине А 
менее 0,25 мм. Дефекты, для которых К > 0,25 мм, условно можно считать тупыми. 

Результаты исследования. Одним из факторов, определяющих эксплуатационные свойства сварных 
конструкций, являются трещиноподобные дефекты —Й такие, как: непровары, подрезы, межваликовые 
несплавления, микротрещины. 

При визуально-измерительном контроле весьма сложно обнаружить острые дефекты с малыми 
значениями измеряемых величин (0,1-0,25 мм). Это объясняется низкой чувствительностью измерительного 
инструмента и человеческим фактором, так как измерения выполняются с помощью конроллера в ручном 
режиме [13]. 

Радиография и ультразвук — два основных неразрушающих физических метода контроля для 
выявления внутренних дефектов. При исследовании особо ответственных конструкций указанные методы часто 
дублируют друг друга, так как имеют различную чувствительность и достоверность. Для выявления острых, 
нераскрывающихся дефектов (например, трещин) предпочтение отдают ультразвуковому контролю” [14, 15]. 
Однако он часто выполняется вручную. Соответственно, на результаты может повлиять человеческий фактор, 
поэтому данный метод представляется недостаточно достоверным. 

Основные сложности при оценке качества изделия могут возникать у специалистов уже на стадии 
подготовки, при выборе метода ультразвукового контроля (например, при использовании зеркально-теневого 
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метода, эхометода или теневого метода). Для выбора необходимо владеть методиками проведения работ, а 
также учитывать характеристики преобразователей в соответствии с ГОСТР 55725-2013'. Такие знания не 
всегда демонстрируют даже специалисты 3-го уровня квалификации, что выявляется при аттестации. 

Специалисты по оборудованию и контролю качества должны уметь правильно настраивать 
дефектоскоп по стандартным образцам (СО-1, СО-2, СО-3) и по стандартным образцам предприятия (СОП). 
Такие шаблоны изготавливаются из материала, аналогичного материалу объекта контроля. Образцы должны 
соответствовать объекту по толщине и (если требуется) по радиусу (при контроле труб малого диаметра). 
Однако практика проведения аттестационных мероприятий свидетельствует о том, что в некоторых случаях 
специалисты не способны качественно выполнять эту работу. 

В процессе контроля важны также скорость и шаг сканирования объекта. При слишком высокой 
скорости сканирования оператор рискует пропустить дефект или неточно определить его координаты, глубину 
залегания и границы несплошности. Нормы скорости сканирования оговариваются в нормативной 
документации. Однако в полевых условиях процедура контроля осложняется температурными перепадами, 
которые зависят от времени года, места проведения и др. Так, например, в зависимости от температуры среды 
меняется температура объекта контроля, и это может сказаться на вязкости контактной жидкости. С другой 
стороны, при низкой температуре требуется подогрев контролируемого участка. Это меняет некоторые 
физические свойства объекта контроля и, соответственно, прохождение ультразвука в изделии и т. п. 

Главная проблема при подведении итогов контроля и выдаче заключений — вольная и субъективная 
интерпретация полученных данных о количестве дефектов, глубине и области их залегания. Многое зависит от 
опыта оператора, его умения отличать сигналы дефектов от ложных или поверхностных. 

Обсуждение и заключения. Необходимо ограничить влияние человеческого фактора для повышения 
достоверности результатов визуально-измерительного и ультразвукового контроля качества сварных 
соединений. Для достижения этой цели считаем целесообразным в процессе периодической обязательной 
оценки квалификации и аттестации специалистов по контролю организовать их дополнительное обучение. 
Преподавателями должны быть специалисты сварочного производства. В краткосрочных учебных программах 
следует обобщить информацию о влияния острых трещиноподобных дефектов сварки на прочностные 
характеристики сварных конструкций и их эксплуатационный ресурс. 

При аттестации необходимо использовать специальные тренажеры для ультразвукового контроля и 
образцы сварных соединений элементов конструкций с искусственно созданными дефектами, имитирующими 
их реальные размеры, форму и расположение в испытательных тест-образцах. Это позволит получить 
объективную оценку достоверности и информативности ультразвукового контроля. 
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О модификации декодера Б(-рршо кодов с низкой 
плотностью проверок на четность 





С. С. Гурский, Н. С. Могилевская 


ФГАОУ ВО «Южный федеральный университет» (г. Ростов-на-Дону, Российская Федерация) 


Введение. Во всех видах цифровой связи применяются методы помехоустойчивого кодирования. Во многих 
стандартах цифровой связи, например вай-фай (англ. \1-Е1) и 5@, используются коды с низкой плотностью 
проверок на четность. Эти коды популярны потому, что для них возможно построение кодеров и декодеров с 
невысокой вычислительной сложностью. Цель настоящей работы — повышение корректирующей способности 
известного битфлиппинг-декодера (англ. Б\-Ирриш®, ВЕ) ГОРС-кодов. Для этого строится модификация 
декодера, позволяющая динамически управлять одним из его основных параметров, выбор которого 
существенно влияет на качество декодирования. 

Материалы и методы. Рассмотрен известный декодер БИ-Иррше двоичных ГОРС-кодов. Некоторые его 
параметры не имеют жесткой связи с параметрами кода. С помощью имитационного моделирования 
исследована зависимость качества декодирования от выбора выходных параметров декодера Б\-Н1рришэ. 
Показано, что на результаты декодирования в этом случае существенно влияет входной параметр декодера — 
порог Т. Разработана модификация ВЕ-декодера двоичных ГОРС-кодов, в которой предлагается задавать порог 
динамически во время выполнения алгоритма в зависимости от степени повреждения кодового слова 
ошибками. Проведен сравнительный анализ корректирующей способности декодеров методом имитационного 
моделирования. 

Результаты исследования. Сформулирована и доказана лемма о максимальном значении порога Т декодера. 
Найдены верхние оценки для количества операций оригинального и модифицированного декодеров. Построена 
имитационная модель, реализующая цифровой помехоустойчивый канал связи. В модели исходные данные 
кодируются заданным ГОРС-кодом, зашумляются аддитивными равномерно распределенными ошибками, а 
затем поочередно декодируются алгоритмом Би-Е@ррше с различными параметрами порога Т и 
модифицированным декодером. По входным и выходным данным оценивается корректирующая способность 
использованных декодеров. Эксперименты показали, что в диапазоне реального уровня ошибок 
корректирующая способность модифицированного декодера выше, чем у оригинального, вне зависимости от 
выбора его параметров. 

Обсуждение и заключения. Доказанная в работе лемма устанавливает верхнюю границу значения порога в 
оригинальном декодере, что облегчает его настройку. По сравнению с оригинальным декодером разработанная 
модификация способна лучше исправлять ошибки. При этом сложность модификации увеличена незначительно 
по сравнению с оригинальным алгоритмом. Отмечено, что качество декодирования модифицированным 
декодером растет при увеличении длины кода и уменьшении количества циклов в графе Таннера, 
соответствующего проверочной матрице кода. 


Ключевые слова: ГОРС-коды, корректирующая способность декодера, динамический порог, двоичный 
симметричный канал, экспериментальное исследование. 
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Введение. В 1963 году в работе [1] Р. Галлагер впервые описал класс линейных блочных кодов, 
проверочная матрица которых содержит малое количество ненулевых элементов. Такие коды принято называть 
кодами с низкой плотностью проверок на четность или ГОРС-кодами (от англ. |10о\/-4епзИу рагцу сБесКк со4ез). 
Для них возможно построение кодеров и декодеров с невысокой вычислительной сложностью. Таким образом, 
при использовании ГОРС-кодов скорость передачи данных существенно не ограничивается. Многие 
современные работы посвящены ГОРС-кодам и их декодерам [2—5]. ГОРС-коды активно используются в 
разных стандартах цифровой связи, например вай-фай (англ. \/-Е!), 5С и оптической связи [6, 7]. Однако, 
несмотря на популярность этих кодов, некоторые связанные с ними задачи требуют исследования и решения. 
Одна из них — построение новых и улучшение существующих декодеров. 

Цель данной работы — повышение корректирующей способности известного декодера БИ-Иррше 
ГОРС-кодов (далее ВЕ-декодер). Для этого строится модификация декодера, позволяющая динамически 
управлять одним из его основных параметров, выбор которого существенно влияет на качество декодирования. 

Материалы и методы. Основными параметрами двоичных ГОРС-кодов являются длина М, 
размерность К и минимальное расстояние кода 4. Информационные слова [М№, К, 4]-кода С — это векторы 
т = (т, ть, ....тк) Е ЕХ, где Ё› — поле Галуа мощности 2, а кодовые слова — векторы с = (с1, со, ..., См) 
ЕЕ” [8]. Удобно задавать ГОРС-коды проверочной (М -К)жМ матрицей Н. Большее количество ее 
элементов — нулевые [1], поэтому удобнее хранить ее не целиком, а запоминая только позиции ненулевых 
элементов по строкам. 

Различают регулярные [9] и нерегулярные [10] ГОРС-коды. В регулярных кодах все строки и столбцы 
проверочных матриц содержат фиксированное количество единичных элементов (К и ] соответственно), иначе 
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код называют нерегулярным. Для удобства проверочные матрицы регулярных ГОРС-кодов будем называть 
регулярными матрицами, а нерегулярных ГОРС-кодов — нерегулярными. 

Регулярные ГОРС-коды обладают рядом преимуществ: легко оцениваемые параметры кода, удобство 
хранения матриц, невысокая вычислительная сложность алгоритмов кодирования и декодирования и пр. Кроме 
этого, декодеры регулярных кодов исправляют ошибки равномерно, в отличие от нерегулярных, которые 
исправляют ошибки в некоторых частях кодового слова хуже, чем в других. Однако задача генерации 
регулярных матриц с заданными свойствами является сложной, часто для ее решения используются методы 
перебора. 

Для обсуждения свойств матрицы Н удобно использовать соответствующий ей граф Таннера @ = 
(У, Е), где Е — множество ребер, а И = 5 УВ — множество вершин, 5 — множество строк матрицы Н, а В — 
множество ее столбцов [11]. Каждый ненулевой элемент Н, стоящий в 1-й строке и /-м столбце, задает ребро, 
соединяющее -ю вершину множества 5 и ]-ю вершину множества А. На рис. | представлен пример регулярной 


проверочной матрицы 3х6 с параметрами К = 4 и] = 2 и соответствующий ей граф Таннера. 
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Рис. 1. Цикл в графе Таннера и в проверочной матрице 


Верхний ряд вершин графа соответствует столбцам матрицы Н, а нижний ряд связан со строками Н. 
Важная характеристика проверочной матрицы Н ГОРС-кода — наличие и тип циклов в соответствующем ей 
графе Таннера. Цикл — это такая последовательность смежных вершин графа, в которой первая и последняя 
вершины совпадают. Длина этой последовательности называется длиной цикла. Минимальная длина цикла в 
графе называется обхватом. Если граф не содержит циклов, его обхват полагают бесконечным. Пример цикла 
длины 4 выделен жирными линиями на графе (рис. 1). 

Возможности по исправлению ошибок зависят не только от основных параметров ГОРС-кодов, но и от 
структуры проверочной матрицы Н. С одной стороны, наличие циклов небольших длин (таких как 4 и 6) 
заметно ухудшает корректирующую способность декодера [12]. С другой стороны, код, которому соответствует 
граф Таннера без циклов, не исправляет ошибок, т. к. его минимальное кодовое расстояние равно 2. Таким 
образом, задача построения проверочных матриц регулярных ГОРС-кодов является многопараметрической. 
При ее решении необходимо следить за основными параметрами кода, а также за циклами в графе Таннера, 
соответствующем проверочной матрице. 

Рассмотрим в удобном виде известный ВЕ-декодер ГОРС-кода С [13]. 


Вход: ГОРС-код С с параметрами [М№, К, а], заданный проверочной матрицей 


Й11 Й2 нии ам 
й й Бе й 
И 2 22 " м | (1) 
П(м-ю1 Им-к? + Им-Юм 
Вектор с’ =сС+е, СЕ С(С Е”), е(Е ЕЁ) — вектор ошибок; р — количество итераций алгоритма; Т — 


пороговое значение. 
Выход: кодовый вектор Сб Е С(< ЕЛ). 
Шаг 1. Положим счетчик 7 равным нулю. 
Шаг 2. Вычислим синдром 5 = С'НТ. Если $ = 0 илит = р, то переходим на шаг 5. 
Шаг 3. Выделим из вектора 5 = (51, 52,....5м_к) единичные координаты, т. е. $; = 1, = 1,(М-К). 


Составим множество [, = {1|5$; = 1}. Вычислим Й' = (11, А», ..., Ам), где 
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п; = ие, Ви. (2) 

Величины Ви, [= 1,..., М следует полагать неотрицательными целыми числами. Таким образом, А’ Е МХ, где 
№ = № {0}. 

Шаг 4. В векторе й' = (11,1», ..., Вл) находим все элементы й: > Т. Среди них выбираем случайный й! 
и инвертируем бит с! вектора с. Добавляем к счетчику 7` единицу и переходим на шаг 2. 

Шаг 5. С := С’. 

Исследования, проведенные в рамках данной работы, позволяют сделать ряд замечаний по ВЕ- 
декодеру. 

Замечание 1. Входной параметр р задает максимальное количество итераций алгоритма со 2-го по 4-Й 
шаги, но декодер может восстановить кодовое слово за меньшее число итераций. 

Замечание 2. При выборе параметра Т нужно руководствоваться следующими соображениями. Если 
известен параметр 4 используемого [М№, К, а]-кода С, то по нему можно вычислить Ё — число гарантированно 


исправляемых ошибок, и тогда количество итераций декодера ограничивается этим значением: 


а-1 
р=:= [> (3) 
Здесь |х| — округление числа х до меньшего целого. Если параметр неизвестен, то его можно оценить с 
помощью границы Синглтона [5] 
а<М-К+1 
и, используя (3), получить 
Ш к 
| 


Замечание 3. Структура декодера такова, что восстановление корректного кодового слова не 
гарантируется, даже если в зашумленном слове С’ = с + е возникло не более & ошибок (3). 


Замечание 4. В литературе для регулярных проверочных матриц в ВЕ-декодере рекомендуется 
] 
выбирать порог Т зависящим от веса ] столбца матрицы Н, а именно Т = >. Для нерегулярных матриц такие 


рекомендации в литературе не приведены. Корректирующую способность ВЕ-декодера может ухудшить 
неудачный выбор порога Т. При его большом значении на шаге 4 декодера в векторе А’ может не найтись 
координаты, превосходящей порог Т, следовательно, не будут исправлены ошибочные биты. При выборе 
малого значения Т на шаге 4 ВЕ-декодера в векторе А’ может появиться несколько координат, значение 
которых превышает порог. Среди них могут быть и координаты, не содержащие ошибку. Таким образом, выбор 
параметра Т может в значительной мере повлиять на качество декодирования. 

Результаты исследования. Сформулируем и докажем лемму о максимально возможном значении 
порога Т. Затем модифицируем ВЕ-декодер таким образом, чтобы порог устанавливался в процессе 
декодирования динамически, и проведем сравнительный анализ оригинального и модифицированного 
алгоритмов декодирования. 

Лемма. Пусть двоичный [М№, К, а]-код С задан с помощью проверочной матрицы Н, имеющей 
фиксированное количество / единичных элементов в каждом столбце. Тогда максимальное значение порога Т 
для ВЕ-декодера такого ГОРС-кода С не может быть больше 

Т=]- 1. (4) 

Доказательство. Пусть из канала передачи получен вектор С’ = С -е, где СЕС — верное кодовое 
слово, е Е ЕМ — вектор ошибок с весом Хэмминга и’(2). Если и’(г) = 0, то на шаге 2 вектор-синдром $ = 0. 
Следовательно, алгоритм перейдет на шаг 5 и вернет С’ в качестве ответа. В этом случае значение порога не 
используется. Если \’(е) > 0, то из регулярности Н вытекает справедливость неравенства Й/ <}, где й, — 
элементы вектора ЙА’. Инвертирование бит с, вектора С’ происходит в алгоритме, только если |! > Т. 
Следовательно, 

Т<# <]. 
Таким образом, формула (4) верна. 
Внесем в ВЕ-декодер изменения, которые позволят определять величину порога динамически, в 


зависимости от степени повреждения кодового вектора в канале передачи. 
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Вход: [М№, К, а]|-код С, заданный приведенной выше проверочной матрицей (1). Вектор с’ = С-+е, где 
СЕС(с Е), е (Е Е”) — вектор ошибок; р — количество итераций алгоритма; Т — некоторое пороговое 
значение, выбранное заранее. 

Выход: кодовый вектор сб Е С(с ЕЛ). 

Шаг 1. Положим счетчик 7 равным нулю. 

Шаг 2. Вычислим синдром 5 = С'НТ. Если $ = (0,...,0) илит = р, то переходим на шаг 7. 

Шаг 3. Выделим из вектора $ = (51,52, ....5м_к) единичные координаты, т.е. $; = 1, [= 1,(М-К). 
Составим множество [, = {15; = 1}. Вычислим В" = (№1, №»,...,Вм), где Ву такой же, как и в оригинальном 
декодере (2). При суммировании величины Йи [=1,...,М№М, следует полагать неотрицательными целыми 
числами. Таким образом, В’ Е №, где № = ММЦ {0}. 

Шаг 4. Инициализируем значение порога Т := тах(й1)1-1..м- 1. 

Шаг 5. Если Т > 0 

Выберем произвольный элемент Йа вектора р’ — такой, что №. >Т. 
Инвертируем бит са. 


Шаг 6. Добавим к счетчику Т единицу и перейдем на шаг 2. 

Шаг 7. с:=С.. 

Замечание 5. Модифицированный алгоритм в целом выполняет меньше итераций, чем ВЕ-декодер, т. к. 
на шаге 4 порог выбирается динамически. Следовательно, декодер не выполняет бесполезные итерации, на 
которых не изменяются биты вектора С’. Значение порога в модифицированном декодере зависит от числа 
ошибок, повредивших кодовое слово, и сразу устанавливается таким, что в зашумленное кодовое слово С’ 
гарантированно вносятся изменения. 

Оценим сверху количество операций сложения, сравнения и присваивания в обоих декодерах. В 
оригинальном ВЕ-декодере [М№, К, 4]-кода С производится р(КК + (М — КЮ)М + 1) операций сложения, р(3ЗМ — 
2К +2) операций сравнения и р((М — КЕ +2М+ 3) +1 операций присваивания. В ВЕ-декодере с 
динамическим порогом — р(КК + (М — К)М +3) операций сложения, р(5М№М — 2К + 3) операций сравнения и 
р((М —К)(К +1) +2М + 4) + 1 операций присваивания. Здесь р — параметр декодера, устанавливающий 
максимальное количество операций, К — вес строк проверочной матрицы кода. Заметим, что при реализации 
алгоритма фактически не используются операции умножения и деления, т. к. на втором шаге для вычисления 
синдрома $ удобно использовать операции сложения вместо умножения. Напомним, что матрица Н обладает 
разреженной структурой, и ее строки удобно хранить в виде списка номеров ненулевых элементов. 
Следовательно, вместо умножения вектора С’ на матрицу Н необходимо суммировать координаты вектора С’, 
номера которых совпадают с номерами ненулевых элементов в соответствующей строке матрицы Н. 

В сравнении с оригинальным алгоритмом модифицированный ВЕ-декодер выполняет больше 
операций, но не значительно: число операций сравнения увеличилось на р(2М + 1), присваивания — на р, 
сложения — на 2р. 

Для сравнительного исследования корректирующей способности оригинального и модифицированного 
алгоритмов декодирования создано программное средство, реализующее имитационную модель двоичного 
симметричного идеально синхронизированного помехоустойчивого канала связи согласно [14—16]. Для 
обеспечения помехоустойчивости в модели использованы ГОРС-коды и ВЕ-декодеры (оригинальный и с 
динамическим порогом). Ошибки в канале моделировались независимыми и равномерно распределенными. 

В экспериментах использованы специально найденные проверочные матрицы, задающие ГОРС-коды. 

Опишем основные параметры этих матриц, используя стандартные обозначения основных параметров кода, а 
также: ] и К — вес каждого столбца и вес каждой строки проверочной матрицы соответственно; 4, 6 —4иб 
циклов в графе Таннера, соответствующем проверочной матрице. 
Регулярная матрица Ну: М = 20, К = 5, =3, К = 4, 4 = 6, в. = 0, «в = 41. 
Регулярная матрица Но: М = 28, К = 7, ] =3, К = 4, а = 6, в. = 0, 5 = 42. 
Регулярная матрица Нз: М = 28, К = 7, ] =3, К =4, а = 6, в. = 0, вв = 29. 
Нерегулярная матрица Н»: М = 32, К = 5, ] = 3, а = 12, в. = 0, ®в = 0. 

С использованием этих матриц были построены кодеки ГОРС-кодов и проведены имитационные 
эксперименты. На рис. 2—5 представлены графики зависимости корректирующей способности построенных 
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кодеков ГОРС-кодов от вероятности ошибки в канале. Обоснование выбора значений порога Т =ТиТ=2в 


ВЕ-декодере см. в замечаниях 3, 4 и лемме. 
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Рис. 2. График корректирующей способности декодеровдля [ГОРС-кодов, заданных матрицей Н> 
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Рис. 3. График корректирующей способности декодеров для 1/)РС-кодов, заданных матрицей Нз 
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Рис. 4. График корректирующей способности декодеров для ГОРС-кодов, заданных матрицей Н. 
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Рис. 5. График корректирующей способности декодеров для ГОРС-кодов, заданных матрицей Н1 


В диапазоне реального уровня ошибок [8, 13, 14] на рис. 2—4 можно наблюдать, что ВЕ-декодер при 
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значении порога Т = 2 показывает лучшие результаты, чем при Т = 1, а модифицированный ВЕ-декодер 
обладает лучшей корректирующей способностью по сравнению с оригинальным. 

Декодеры демонстрируют схожую эффективность при малых значениях длины кода, однако при ее 
102 увеличении модифицированный декодер показывает лучшие результаты. Так, при вероятности ошибки в 
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непомехоустойчивом канале 0,05 различие в вероятности ошибки в помехоустойчивом канале при 
использовании ВЕ-декодера с порогом Т = 2 и Т=1 составляет от 0,005 до 0,03 в пользу использования 
большего значения порога. Если же задействовать ВЕ-декодер с порогом Т = 2 и модифицированный декодер, 
это различие колеблется в зависимости от [ГОРС-кода в интервале от 0,001 до 0,003. При вероятности ошибки в 
непомехоустойчивом канале 0,1 вероятность ошибки в помехоустойчивом канале при использовании ВЕ- 
декодера с порогом Т = 2 меньше, чем с порогом Т = 1, на величину от 0,001 до 0,02. При использовании ВЕ- 
декодера с порогом Т = 2 и модифицированного декодера это различие колеблется в зависимости от ГОРС- 
кода в интервале от 0,002 до 0,01. 

Оба декодера чувствительны к количеству циклов в графе Таннера, соответствующему проверочной 
матрице ГОРС-кода. Чем больше отношение количества циклов к общему количеству элементов в матрице, тем 
хуже исправляет ошибки любой ВЕ-декодер. При проведении экспериментов было интересно выяснить, можно 
ли так увеличить количество циклов в матрице, что модифицированный декодер покажет худшие результаты по 
сравнению с ВЕ-декодером. Опытным путем была найдена такая матрица — Н:, содержащая 41 цикл длины 6. 
Результаты исследования корректирующей способности декодеров для этой матрицы показаны на рис. 5. 
Заметим, однако, что матрица Н> содержит еще больше циклов длины 6, а именно 42. Принципиальное отличие 
матриц Н1 и Н> — в плотности единиц: 

— в Н1 — 60 единичных элементов на 300 элементов матрицы, 

— в Н> — 84 единицы на 588 элементов матрицы. 

Напомним, что особенностью ГОРС-кодов является разреженная структура проверочной матрицы, 
поэтому Н> более типична для ГОРС-кодов. 

Обсуждение и заключения. В работе рассмотрен декодер 6-Е ррше двоичных ГОРС-кодов. Даны 
рекомендации о выборе таких входных параметров декодера, как порог и количество итераций алгоритма. 
Сформулирована и доказана лемма о максимальном значении порога декодера. Разработана модификация ВЕ- 
декодера двоичных [ОРС-кодов, в которой предлагается задавать порог динамически во время выполнения 
алгоритма в зависимости от полученного синдрома. Для оригинального и модифицированного декодеров 
найдены верхние оценки количества операций. Эти оценки показывают, что модификация усложняет декодер 
незначительно. Проведенные имитационные эксперименты демонстрируют лучшую корректирующую 
способность модифицированного декодера по отношению к оригинальному. Эксперименты также показали 
зависимость качества декодирования от степени разреженности матрицы и количества циклов длины 6 в графе 
Таннера, соответствующего проверочной матрице ГОРС-кода. Таким образом, возникает задача построения 
проверочных матриц с малым количеством коротких циклов, что является предметом дальнейших 
исследований. 
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итодисйоп. Т№$ рарег 15 4еуме т Фе 4еуеортепе оЁР а ше#о юг езитайпе Фе сштеп ите-Кедаепсу 
рагатеег$ оЁ еасВ оЁа $61 оЁ зиииапеоч$[у ап ш4дерепдепйу ореганпе оепегаюг$ ш Фе гафо @есготс зузет. А 
сепега| сазе 15 сопз14егеа, ш \мсВ Фе деулайоп оЁ ие уамез оЁ 1е поплтпа| сепегаюг рагатеег$ Нот Фе аззитеа 
уа[аез ог Фезе рагатеегз 1$ д&егтише пой ошу Бу гапдот, БЕ а150 Бу 10п2-{егт Недиепсу демайоп. ТЬе \огК 
оБуесйуе 15 №0 оепегаПте Фе тефо4 Юг езитайп» Фе ите-Недиепсу рагатеег$ оф $12па1$ (юпе-егиа потита! 
Кеадиепсу ап сштепе Недиепсу деумаНоп Кот Фе поттта| уаме) Базе4 оп Фе зппиЦапеоч$ теазигетепЕе оЁ фе 
рВазез о $1епа1[5 сепеге ш фе зузет оГ ш4дереп4депЙу Рлиснопше сепегаогз. ТВе гезеагсВ фазК 1$ №0 сопзл4ег а 
зузет оЁ зиииНапеоч$[у ап ш4ерепдеп Йу орегайпе сепегаюгз. ЕасВ оГ Фет сепегаез$ Вагипотис 51епа1$, \позе 
ите-Недиепсу рагатеег$, зас аз фе ауегазе Недаепсу, аге сопзапё фигие а сейаш пищегуа| оЁ обзегуайоп. Ви 
Беге\/И, безе ите-Недиепсу рагатаег$ аге Клпо\п УИ тпзиРАслеп{ ассигасу 4ие ю Фе шЙчепсе оГежегпа1 асюг$ 
(сВапоез ш фетрегайхге, ргеззиге, зирр!у уо[асе, ес.). [Ц 1$ гедапей №0 оБаш езИтаез$ оЁ Фе Ите-Кедаепсу 
рагатеег$ ог 51епа[5 (фе Фигайоп оЁГ Фе теазиние ищегуа|, уаез ог Фе 1оп2-еги Недиепсу ап Фе запдага 
деу1айоп) Нотт Ще гези {$ оЁ теазигетеп{$ ог Фе рВазез ог $12па15 огте4 Ъу Ше сепегаюг$ аё теазигие ищегуа|5 
Беопеотас фю Ше обзегуаноп ищегуа|, уф уЛасЬ Фе ауегасе Недаепсу гетатз сопзапй. 

Маептай апа Мейфосб. А зузет о? зипаЦапеоч$[у ап ш4ерепдепЙу РапсНопио сепегаюг$ 15 сопз14еге4. ТВе 1опз- 
{егт уае о? $1епа| Недиепсу Юг еасП о{ Фе сепегаюг$ оуег Фе обзегуаНоп ищегуа| гета1п$ сопзбапь Би (1$ Кпо\/п 
\И зоте таготи ОР еггог. Оигие фе обзегуаноп пщегуа|, зеуега| плеазигетеп{ оЁ Ше 51епа| рВазе оф еасВ ог е 
сепегаогз аге ре огтед. А{ Пе зате ите, Ве сиггеп{ уа[ез оЁ е з1епа| Недаепсу ап Фе 4иганоп о Фе теазиги? 
и\егуа| Бауе гапдот Чемайноп$ Нот Ше 1опэ-еми уаез, ап ЮПо\ Фе погта| а151биНоп |а\у у тего 
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Метод оценивания длительности временного интервала с использованием 
одновременных измерений фазы в системе одновременно и независимо работающих 
генераторов 


Д. Д. Габриэльян, А. А. Костоглотов', О. А. Сафарьян’, Л. В. Черкесова”, О. В. Дворников 

' ФГУП «Ростовский-на-Дону научно-исследовательский институт радиосвязи» (г. Ростов-на-Дону, Российская Федерация) 
? ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 
ОАО «Минский научно-исследовательский приборостроительный институт» (г. Минск, Республика Беларусь) 


Введение. Статья посвящена разработке метода оценивания текущих частотно-временных параметров каждого 
из совокупности одновременно и независимо функционирующих генераторов в составе радиоэлектронной 
системы. Рассмотрен общий случай, при котором отклонение истинных значений номинальных параметров 
генераторов от предполагаемых значений указанных параметров определяется не только случайным, но и 
долговременным отклонением частоты. Целью работы является обобщение метода оценивания частотно- 
временных параметров сигналов (долговременной номинальной частоты и текущего отклонения частоты от 
номинального значения) на основе одновременного измерения фаз сигналов, формируемых в системе 
независимо функционирующих генераторов. Постановка задачи — рассматривается система из одновременно и 
независимо функционирующих генераторов, каждый из которых формирует гармонические сигналы, частотно- 
временные параметры которых, такие как средняя частота, являются в течение некоторого интервала 
наблюдения постоянными, но известными с недостаточной точностью из-за влияния внешних факторов 
(изменение температуры, давления, напряжения питания и т.д.). Требуется по результатам измерений фаз 
сигналов, формируемых генераторами на измерительных интервалах, принадлежащих интервалу наблюдения, в 
пределах которых значения средней частоты остаются постоянными, получить оценки частотно-временных 
параметров сигналов — длительности измерительного интервала, значений долговременной частоты и 
среднеквадратического отклонения частоты. 

Математическая модель и методы решения. Рассматривается система одновременно и независимо 
функционирующих генераторов. Долговременное значение частоты сигнала для каждого из генераторов на 
интервале наблюдения остается постоянным, но известным с некоторой погрешностью. В течение интервала 
наблюдения проводится несколько измерений фазы сигналов каждого из генераторов. При этом текущие 
значения частоты сигнала и длительности измерительного интервала имеют случайные отклонения от 
долговременных значений и подчиняются нормальному закону распределения с нулевым математическим 
ожиданием и известной величиной дисперсии. 

Оценивание частотно-временных параметров по результатам измерения фаз сигналов проводится с 
использованием многомерной функции правдоподобия. В качестве оценок случайного отклонения частоты 
сигналов и длительностей интервалов измерения, выбираются значения, доставляющие максимум функции 
правдоподобия. Нахождение максимума проводится на основе решения переопределенной системы линейных 
алгебраических уравнений. 

Результаты исследования. Разработана новая математическая модель и численно-аналитический метод 
определения частотно-временных параметров сигналов, учитывающие как долговременное постоянное 
отклонение частоты, так и кратковременные отклонение, носящие случайный характер. 

Обсуждение и заключения. Полученные результаты могут быть использованы при разработке и создании 
информационно-измерительных и —’информационно-телекоммуникационных систем, в том числе 
территориально распределенных систем. Получаемые оценки частотно-временных параметров позволяют 
повысить стабильность частоты сигналов и соответственно повысить точность проводимых измерений и 
качество передачи информации. 


Ключевые слова: генераторы высокочастотных колебаний, высокочастотные радиосигналы статистический 
метод стабилизации частоты, стабильность частоты, метод наименьших квадратов. 
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Шеодисйоп. ТБе соп$фапё ого\/Ь оЁ гедатетеп{5 Юг аа41о ап у14ео шЮюгтайоп {гап$111$510п, Юг даа 
ГКогтайоп ап4 1гап$1115$10п зу$еил$, песеззцайез ап шсгеазе ш забИу оЁ Фе ите-Неаиепсу рагашаег$ оЁ $1°1а15 
сепегаййе ш гад1о-@есготис зу$етз (КЕЗ) [1]. Тве шЮюгпайоп-@есоттитсаной зу$бетз ап ажа-теазитие 
опез аге {ур1са| ехатр[ез оЁ засВ зузепа$. Опе оЁ Фе сотропеп5 ш ргоу1Аше Б1>В еЁЯслепсу оЁ засВ зуз{еит$ 15 


Габриэльян Д. Д. и др. Метод оценивания длительности временного интервала с использованием одновременных измерений фазы 





аззостайед \уив Фе Тогтайоп оР $1епа1$ \уив Б1еШу за е ите-Недиепсу рагатаегз. ТВе 1авег 1$ гееуап ш 
рагис\аг, Гог Фе зуетз оЁ огтайоп 1гап$ип11$$101 УИ сотр] ех $1епа[5, гадаг ап гад1о паутеайоп зу${еттз, 
ап аи410о ап4 у14ео шРогтаноп ап$111$5101 зу$ета$ аз мей. 

Сштепйу, фе таш арргоасВ ш Фе сгеайоп оЁ ВЕ$ 15 а тод\аг сопзбгасйоп рипстре, уБлсВ даегилитез Ве 
и5е ог КисНопаПу сошр!&е4 ЫосКз$ ап 4еу1сез соппеце4 о еасВ оег. ТВе аррПсайоп оР 1$ арргоасЬ 1еа4$ ю 
ргезепсе оЁРа 1агое питбег оЁ №121 -Кедиепсу 5$1епа| (НЕ $1епа1) сепегабог$ ш апу гад1о-ейесгоплс зубет. А{ Ве зате 
ите, РисНопа|! сотр!е{епез$ оЁ еасВ оЁ Фе 4е\мсез шса4де4 ш Фе зубет епаез © сопз14ег Фе НЕ-$1епа1 
сепегабог$ ог езе 4еу1се$ аз ип 5 орегайи® зипиЦапеочц$[у ап4 ш4дерепдепйу оп а сецаш те ищегуа1. 

Резрие Ше Тас{ Фа 51епа| сепегафогз ш уапоиз$ 4е\мсез Вауе Ч1егепе рагатеегз, Ше ргезепсе оРа 1Тагое 
питьег оЁ зппиКапеочц$[у ап шдереп4еп у РлисНоп1ие оепегабог$ ш фе КЕЗ ргоу14ез$ даегилите Фе сиггеп{ уа[ез 
ог фе ите-Неаиепсу рагатеег$ оЁ Везе 51епа1$ @гоиеВ теазигие апа за6зеааепа ргосеззше оРе оепегае4 5$1°1па15 
рВазез' [2—6]. Те 1айег такез И роз е ю езитае Ше сие уаез оЁ 1е ите-йеаиепсу рагатлейегз \ИВ Меег 
ассигасу. ш фаги, Фе оМаштед уашез оЁ Ме ите-Недиепсу 51епа| рагтате{егз аге аззослайеа у е рагатеег$ оЁ Фе 
сепегайог$ Ша аге Кюгише Фет. ТБа{ аПо\уз ейпег ю уаБШте Ше сепегаютз Неаиепсу, ог ю сопз14ег $ сиггеп 
уае дигио Ше забзеацеп: $1опа1 ргосеззште [2—6]. 

\/е|-Кпомт зо[иНопз ю Фе рго ет оЁ езитайпе фе сиггепЕ ите-Ееацепсу рагатаегз {аКе шю ассоцпь аз 
а пе, оту гапдот сотропеп оЁ Ще Недиаепсу демайоп 4ие ю Фе шйНоиепсе оЁ уатом$ с@югз. ш 15 сазе, фе 
уа[ез ог зисВ рагатщег аз бе попа! (теап) Кедаепсу оЁ Фе сепегаюг аге сопз14еге Кпо\тп. АЕ Фе зате ите, ш 
тапу сазез, 1оп2-{егт Недиепсу 4емлайоп$ аззослаеа уу Бой Фе шЙчепсе оГ ежегпа1 Гасбог$ ап Фе {1есбпо]о2у оЁ 
сепегабог ргодисйоп аге по{ апа[уте ш Фе рарег$ ип4дег ехаплтаноп. 

Тра$, фе зоийоп ®ю Фе рго ет о{ езитайи> Фе ите-Недиепсу рататаегз оф 512па15, засВ аз фе потитта| 
(теап) Недиепсу ап4 гапдотп 4еулайоп$ Нот Ве ауегасе гедчепсу оЁ сепегафогз, 1$ а сваПепое. 

Тре рарег об]есйуе 1$ Фе сепега|таНоп ог фе шефоа оГезитайпе фе ите-Неачиепсу з1опа| рагатаетз (е 
1опе-{егт попа! Недиепсу апа Фе сиггет{ Недаепсу деуланоп Кот Фе потлша| уаше) Базе оп е зппиКапеоц$ 
теазигетеп о Ве 51епа| рВазез Гогилие ш фе зубет оЁ шдерепаеп Иу РлисНопше оепегаюгс. 

Аз ш Ще рарег [3], 1е{ из сопзлаег 1е зубет оЁ М№М+1 зепегаюог$, еасВ о? \мсВ Югтз Фе Баглпогтис 51епа[5, 
уВо$е ите-едаепсу рагатщег$, засВ аз ауегазе Недиепсу ©, , аге сопуапь Баё Кпоуп ул шзаЁРйсет ассигасу 


ие ю Ше шНиаепсе оЁ ежегпа! асбюгз (сВапоез ш фетрегаехге, ргеззиге, заррЙу уоЦасе, ес.) диние а сецаш 
обзегуайоп ищегуа|. 

Е 15 гедише4 №0 обаш езитаез ог Ше Ите-Недаепсу рагатеег$ оЁ 51°1па15 (фе аиганоп оЁ Фе теазагие 
и\егуа|, Неадиепсу уаез ап КМЗ (гоо-теап-зацаге, ог зап4ага, деманоп) оЁР 51е1па15) Нот Ше гези5 оЁ 
теазигетеп5 оЁ Фе рБазез оЁ е 51епа|$ Югшеа Бу Ше оепегабюг$ а Фе теазитие пищегуа15 Беопеоше © Фе 
обзегуайоп ищегуа[, ут ууРсВ Ше ауегасе Кедиепсу уа[аез гетали соп$апй. 

Ма "етайса| поде апд зом@оп тео. Еог Фе зует ог М№М+1 сепегаюгз ип4аег сопз14еганогп, [е{ и$ 
ргезеп+ Пе $12па[5 Гогте4 Бу еасВ оЁ е сепегафотс, аз 

5, (Е) = А, с0з(®, -1+ф,), п=1....,М-+1. (1) 

Озше опе оЁ Фе сепегаюогз (Бегешайвег Фе №М+1 сепегаог), ме уШ Юпи шеазипие ищегуа|$ оЁ поплта| 
Чиганоп #°” (т=1.....М), засб Фа 4итпе а Фе зрестИеЧ пиегуа!$, е ауегасе Недиепсу апа ВМ$ оё Ше 
Кеадиепсу оРеасВ ое сепегаог$ сап Бе соп$14еге4 сопзапй. 

\У!е гергезепе Ме пеазигеЧ уа[иез оЁ 5151а1$ рвазе Ф, „оЁ еасН оЁ М зепегаюгз (п=1..... №) оп фе т-@ 


теазигие ищегуа1 (т =1,..., М) аз ЮПо\з: 
Ф, „= (©, +4, (1 +4Е,), п=Ь,..,М, т=Ь,...М, (2) 
уреге До, „-— 15 гап4от Чдеулайоп оЁ Ше п-@ репегабюг Недиепсу оп фе т-@ теазигие ищегуа! Нот Фе 


ауегасе уаше ®,; #” ап ДЕ, аге, гезресиуе1у, потлта] 4игайоп ап де\ланоп Нот Фе попа] уаше оЁ Пе т-® 
теазигие> пфегуа1 дигайоп. 
Еог фе теап уаме НКедаепсу, фаКие шю ассоци Ше тодеги ргодисНоп {есбпо]огу ог НЕ-сепегаог$, \е 
сап утие 
<, = ©” +60, бы <”, п=|,..., М-+Т, (3) 


) 


\’Веге пе уаиез с” аге Кпо\уп, ап до, , \ИВ ассоип Юг 1е шНиепсе оГежегпа]| асюг$, аге ипкпо\уп. 


Такше шю ассоци{ Ше ехргезз1оп$ (3), фе г@аНоп (2) фаКез Фе юг 


" Сабмейуаю ОО, Ргусипоу АС, КвиюизеуУУ, её а|. Сепегаюг Неаиепсу заб Штаноп тево4. ВЕ Райепи по. 2219654, 2003. (ш Виз.) 
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Ф, „= (+6, +Ло, (+4, ), п=Ь,...М, т=Ь...М. (4) 


АЁег Ппеаг1тайоп (Ще {егтз$ 0%,.ДЕ ап До, „-ДЕ, , аге 415саг4е4), Фе ехргеззлоп (4) епа ез {0 езйтайе 


Фе гапдот деуланоп ое Кедиепсу оР Ве и-@ сепегабог оп Фе т-@ теазигетепе ищегуа| аз оПо\5: 


Ди ии и=|,..М, ШЕЬ...М, (5) 
| [А 


0 0 0 
упеге Ф'") Е. 


Тре Недаепсу демайоп$ Ао’ (п=1....,№М,т=1.....М) ап4 демаяоп$ оЁ теазигетепе пфегуа| дигайоп 


п,т 


ДЕ (т=1..... М) ЮПо\уи Фе погта] 415 1бийоп Тау [7]: 


| (ло, 1 (ль, 
р(4в.,)- 2по“ вы 25°. 6 ) РМ.) 2По, ыы 25 ое о | _ 


у\Пеге о” аге Кпо\уп уашез оЁ Ве ге!айуе ша БИиу ое Нгедибйсу оРе и- зепегабюг (и =1...., №). 


Тре олуеп г@аНоп$ аезсге ше таетайса| под] ог Ше зузет оЁ оепегаюг$ фа{ псйоп зипиЦапеоч$[у 


ап шдереп4епйу. 
зошНоп те#фод. ТаК1их шю ассочпЕ Фе геаНопз (5) апа (6), \е сап \тце Фе шии4итеп$1опа1 1огатИис 


ПКейвоо4 Рапсйоп: 


мм Ф -Ф® бо 1-ю. Е) 
[ (бо, №) =-У' У 1 ш 2 +, ‚ ©. иво, М.) (7) 


т= п= 2(о” . 00 . #0 ) 
\расЬ шс[а4ез Фе уесюг$ 0% апа Де, позе еетепт аге ипКпо\и уашез, гезресиуу, б®, (и=1...., №) 
апа ДР (т=|,...,М ). 
ТБе езйтаез$ до, ап@ ЛЕ, аге юипа Нот Фе тахипит соп@ оп (7) ап4 сотгезроп@ ю Фе зойоп оЁ Фе 


зузет оЁ Ппеаг а[2ебга1с едцаноп$ (ЭГАЕ). 


ОГ, (60, +) 
о =О, И=[...., М, 
000, (8 
ОГ (66, +) 
о =0 т=|.....М. 
ОДЕ, 
Тре едцаноп$ ш (8) Вауе фе огт 
(Ф =Ф®—д0. 0 ==О № о 
й,т пт ПИ т = т т =0, арм, 
(0 ‚0®. о 
(9) 


—_ @Ф® — (0) со (0) 
(Ф,„-Фо 60, 1-0”. в 
о =0, т=1,....М. 


| р . | | | | | | ) 
А›, ыы А В, ? 


—2 . . 
уПеге А’ 1$ Фе Ыоск эй 4птепзюп$ №М.МхМ \миН еетеш а” =(0°.®°”) ; А. 15 Фе Ыоск ми 
11 1,2 


п-тьт п п 


_2 —1 
Читепз10$ М.М хМ И е!етеп @`”), = (с”) (© и 


А 


Т 
‚, 18 Фе Моск у Фитепзюп$ М. М х М ми @етеп а”) = в. 


‚т 


= 
А, 13 Ше Моск уп 4ипепзюпт$ М.М хМ ми @етеш а.) = (с | ©) 


п-т,т 
Ф _Ф<® 
(0 (0)\? +00)” 
о, ‘®, | -Ё 


(0) 
пт 


В, 15 Фе Моск \иИВ 4итепзюопт$ №. М х1 \иВ @етет В” = 


пт 


В, 15 Ше Моск ий 4ипепзюоп$ №М- М х1 \и еетет 6 = (о 
о: = 


Г 1$ 1е $11 оЁе таблх шапзроз1 оп. 


Габриэльян Д. Д. и др. Метод оценивания длительности временного интервала с использованием одновременных измерений фазы 





Тре зубет оГ едиайоп$ (9) сощашз № оЁ Фе ипКпо\п сотропеп5 до апа М оЁ Фе ипКпо\п 
е]етеп ДЕ’. ТВе {1юа питбег оЁ теазигед уамез оР $12 па1 рВазе 15 №. М. Те езитаноп оЁ Феш уа[чез, {аК1пе шо 


ассоип ше теазигетепе еггог$, зо Ъе сагле ош Бу Фе 1еа${-зааге тефоа (Г$М) [8-10]. ш 1$ сазе, Фе 
сопЧоп №М.М > М+М таз Бе шеф ап Фе патфег оЁ теазигетепе пегуа[5 плиз зайз$Ёу фе пехё сопа1оп: фе 
пипабег оЁ теазигетет М позе ВИН] бе сопаюп М > №М/(М-П. 


Тре гергезещаноп о ЗЭГАЕ чзше Ще ехрге$101$ (9) ап4 (10) 4ейпез Фе питепса|-апа]уйса! тефоа г Ше 
ргоет зо[уш2. АП еетеп5 ог Фе таблх Вауе ап апа!уйса1 гергезетщайоп. А Фе зате ите, уБеп раззте тю Ше 
зузет оЁ погта| едаайопз$, аз а гше, изе4 ш фе ГМ, оМашше Фе апа[уйса| ехргез$101$ 15 а150о поЁ ЧИсий. 
Нозеуег, фе шуегз1оп оЁ]агое-дитепт$1опа] тайлх сап Бе регюгте4 ошу изше Фе питепса| тефо4$ [11-13]. 


Тре езитаез (бе, ) апа (АЕ, ) ораше Нот йе зоийоп ю Ше зузет оЁ едцайоп$, гергеземе4 Бу Ше 


ехргез51оп$ (9), аПо\м из ю даегише слигепе ите-Недиепсу рагате{ег$ оЁР Фе сепегабог$ ап4 $12па15 югише Бу 
ет. 

Везеагсв Вези!65. У\е ргорозе а тафетайса| плоде! дезст1те а зипаКапеоц$ ап ш4ерепдет РлисНопе 
ог Ше оепегаюг зузет, ап а питегса|-апа]уйса| тебфо4 г деегишше Фе Ите-Кедаепсу рагате{ег$ оЁ Фе 
$121а1$ м ассоциЕ юг Бо Фе 1опэ-егт сопзап{ Недаепсу аемлайоп ап Фе звой-егт аеулайоп оЁ тапдот па@хе. 
ТЬ1$ арргоасв епа]ез: 

— © еуамжме, ассог4ше фю Фе гези5 оЁ теазитие фе рВазез оф 51епа1$ Ююгте4д Бу зипаЙапеон$ Ту апа 
шДерепдеп у Рпсйопшо сепегаотз, по{ оу гапаот 4деулайоп$ ш Фе Неачепсу оЁ Фе сепегме4 510па[5, ап ю 
оБбаш Ше езИтаез$ оЁ Ше ауегасе |опо-{егт Недиепсу оЁеасВ ое сепегаюогз; 

— ®ю ехс4е ап шуаБИиу шИчепсе ог Фе ите пщегуа! дагайоп оР теазигетеп{5 оп Фе тез те 
езитаез оЁ е 1епа| {ите-Кедаепсу рагатаегз [14-16]. 

015си5510п апа Соп$10п$. ТБе гези{$ оашеД сап Ъе изе4 ип4ег Фе 4еуеортепЕ ап сгеайоп о{ да- 
теазигие ап шЮгтаноп-{еесотититсайоп зузетл$, шсаатше Фе оеоэтаршсаПу 41516 ще4 зу%етаз. ТБе гез по 
езитаез$ ог фе Ите-Кедиепсу рагатаег$ ргоу14е шсгеазш» Фе 51епа| Недиепсу зЗаБИИу ап4, ассогат®у, 
пиргоуе Фе ассигасу оР Фе теазигетеп5 ап Ше дааЩу оЁ шРгтайоп гап$111$$101. 
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